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ABSTRACT 


The  theory  and  practice  of  far-infrared,  two-beam 
interferometry  is  summarized.  A  description  of  the  instru¬ 
mentation  and  performance  of  the  Aerospace  Corporation  two- 
beam,  double-beam,  far-infrared  lame  liar -grating  inter¬ 
ferometer  is  given.  Detailed  experimental  measurements 
and  analyses  of  the  pure  rotational  spectra  of  the  atmospheric 
molecules  nitric  oxide  (NO)  and  water  (H^O)  are  presented. 
Also  discussed  in  depth  are  measurements  of  the  optical 
constants  of  some  important  far-infrared  window  materials 
and  the  fabrication  and  performance  of  narrow  bandpass 
filters  for  the  far-infrared.  Included  is  a  discussion  of 
noise,  experimental  errors,  and  problems  in  the  practice 
of  two-beam  interferometry.  Finally,  a  summary  is  given 
of  the  computational  programs  necessary  and  convenient  for 
reducing  the  raw  interferometric  data  to  their  more  readily 
used  spectral  form. 
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1.  INTRODUCTION 


It  is  the  purpose  of  this  report  to  summarize  the  experience  and 
results  obtained  with  the  Aerospace  Corporation  two-beam,  double-beam 
lamellar- grating  interferometer. 

The  history  and  status  of  two-beam  interferometry,  sometimes  called 
Fourier  transform  spectroscopy  or  autocorrelation  spectroscopy,  up  to 
October  1965  has  been  summarized  by  Loewenstein  (Ref.  1),  and  the  inter¬ 
ested  reader  is  referred  to  his  article  for  a  comprehensive  survey. 

Briefly,  Fourier  transform  spectroscopy,  which  had  its  inception  in  the 
works  of  Michelson  (Ref.  2),  has  proved  to  be  a  valuable  tool  in  the  in¬ 
vestigation  of  the  far-infrared  region  (FIR)  of  the  electromagnetic  spectrum. 
Fellgett  (Refs.  3,  4)  has  shown  that,  when  detector  noise  is  the  limiting 
factor,  the  gain  in  signal  to  noise  (S/N)  of  an  interferometer  over  a  conven¬ 
tional  single-slit  dispersive  spectrometer  is  n/N/2,  where  N  is  the  number 
of  spectral  elements  in  the  region  studied.  This  is  called  the  multiplex  or 
Fellgett  advantage. 

There  are  two  types  of  two-beam  interferometers  used  in  the  infrared 
(IR).  The  classic  Michelson  interferometer,  which  utilizes  division  of 
amplitude  to  produce  the  two  interfering  beams,  has  been  developed  by 
Gebbie  (Ref.  5)  for  use  in  the  FIR.  J.  Connes  (Ref.  6)  has  also  used  a 
Michelson  for  low-level  intensity  studies  in  the  middle-  and  near-IR. 
Recently,  J.  Connes  and  P.  Connes  (Ref.  7)  have  reported  on  truly  spectac¬ 
ular  spectra  of  Mars  and  Venus  in  the  near-IR  using  a  Michelson  interfero¬ 
meter.  In  the  FIR,  Bell  (Refs.  8,  9)  and  Chamberlin,  et  ai.  (Ref.  10)  have 
exploited  the  fact  that  in  a  Michelson  interferometer  the  interfering  beams 
are  well  separated  in  space  before  recombination.  This  allows  placing 
reflective  or  absorbing  samples  in  one  arm  of  the  Michelson  interferometer 
and  thus  obtaining  phase  information. 


A  second  type  of  interferometer,  the  lamellar  grating  interferometer, 
produces  the  two  interfering  beams  by  dividing  the  wave  front  (thus  avoiding 
experimental  difficulties  typical  of  a  beamsplitter).  This  interferometer 
has  been  developed  by  Strong  and  Vanasse  (Refs.  11-13),  A  comparison 
between  these  two  types  of  interferometers  and  a  conventional  grating 
spectrometer  has  been  reported  by  Richards  (Ref.  14).  Recently  a  high 
resolution  lamellar  grating  has  been  described  by  Hall,  et  al.  (Ref.  15). 

It  is  an  FIR  instrument  which  has  a  double-beam  differencing  mode  of  opera¬ 
tion  that  helps  alleviate  the  dynamic  range  problem. 

The  present  report  consists  of  a  brief  survey  of  the  theory  of  two- 
beam  interferometry  with  emphasis  on  the  experimental  techniques 
necessary  for  obtaining  the  interferometric  data.  A  brief  discussion  of 
some  idealized  cases  is  presented  to  familiarize  the  reader  with  interfero- 
grams  (the  raw  data  from  interferometers)  and  their  direct  interpretation 
in  some  simple  cases.  Next,  the  instrumentation  of  the  two-beam,  double¬ 
beam  lamellar  grating  is  described  and  the  advantages  and  disadvantages  of 
its  double-beam  differencing  mode  of  operation  discussed.  The  next  section 
consists  of  experimental  results  obtained  with  this  instrument.  Included 
are  the  measurement  and  analysis  of  pure  rotational  spectra,  linewidth 
measurements  and  pressure  shifts,  index  of  refraction  measurements  of 
gases,  optical  constants  for  solids  in  the  FIR,  and  the  fabrication  and  per¬ 
formance  of  narrow  bandpass  filters  for  the  FIR.  Some  of  these  results 
have  been  published,  others  are  being  prepared  for  publication,  and  some 
are  preliminary  in  nature  (and  in  these  cases  the  conclusions  drawn  must 
be  considered  somewhat  speculatory).  The  next  section  deals  with  some  of 
the  more  important  experimental  errors  encountered  in  two-beam  inter¬ 
ferometry  and  how  these  errors  are  reflected  in  the  spectrum.  Also  in  this 
section  are  discussed  some  of  the  problems  we  have  encountered  which  are 
peculiar  to  two-beam  interferometry  (and  some  peculiar  to  our  instrument) 
and  how  their  eventual  solution  will  undoubtedly  lead  to  a  higher  S/N  ratio 
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for  FIR  interferometry.  The  last  section  deals  with  the  computer  programs 
and  options  that  are  necessary  and  useful  to  convert  the  interferogram  into 
the  more  familiar  spectrum.  Included  is  a  discussion  and  listing  (FORTRAN) 
of  the  conventional  "direct  sum"  Fourier  cosine  transform  program  and  the 
less  flexible,  but  much  faster,  Cooley- Tukey  program. 


II.  THEORY 


There  are  many  excellent  treatments  (Refs.  6,  12,  16-20)  that  cover 
with  varying  degrees  of  rigor  the  theory  and  practice  of  two-beam  interfer¬ 
ometry.  In  this  section  we  shall  only  summarize  the  theory  and  then  discuss 
a  few  simple  theoretical  examples  which  will  be  useful  in  understanding  the 
experimental  results  in  Section  IV. 

Consider  a  beam  of  monochromatic  electromagnetic  energy  which  is 
passed  through  a  two-beam  interferometer.  The  intensity  F'(x)  detected 
would  have  the  following  form  (Ref.  21). 


F'(x)  =  1+1  cos  Zimv 


(1) 


where  I  is  the  average  intensity  (i.  e.  ,  the  average  of  F'(x)),  x  is  the  optical 
path  difference  generated  by  the  interferometer  of  one  of  the  interfering 
beams  with  respect  to  the  other,  and  v(=  X  )  is  the  wavenumber  of  the  radi¬ 
ation.  If  instead  of  monochromatic  energy  we  have  a  band  of  wavenumbers, 
then  Eq.  (1)  may  be  generalized 

v)dv  + 


I(  v)  cos  2wxv'dv 


(2) 


The  first  integral  is  a  constant  and  is  the  total  intensity  of  the  spectrum  I(v) 
(which  also  equals  F'feo)).  In  terms  of  F'(x)  it  is  seen  to  be  (1/2)F'(0)  by 
evaluating  Eq  (2)  at  x  =  0.  Thus 


F'(x)  -  j  F’(0) 


/ 


00 

I(v)  cos  2trxydv«F(x) 


(3) 
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In  the  literature  F(x)  is  called  the  interferogram  or  interferogram 
function,  as  is  F'(x).  To  avoid  confusion  we  shall  call  F(x)  the  interferogram 
function  and  F'(x)  the  interferogram. 


F(x)  is  the  Fourier  Cosine  Transform  (FCT)  of  I(v).  It  then  follows 
that  if  F(x)  is  Fourier  cosine  transformed,  the  result  will  be  I(^).  That  is, 
aside  from  a  multiplicative  constant  and  recognizing  that  F(x)  is  an  even 
function  of  x. 


I(v) 


/ 


00 

F(x)  cos  2irxvdx 


o 


(4) 


In  this  form  Eq.  (4)  is  of  limited  practical  importance.  In  general, 
if  F(x)  is  known,  a  rather  large  computer  is  needed  to  calculate  I(v),  and 
Eq.  (4)  must  be  put  into  a  form  for  which  a  computer  program  can  be  written. 
Analog  computers  have  been  used  but  usually  with  digital  input  data.  It  is 
usually  digital  computers  that  are  used  to  solve  Eq.  (4)  for  I(v).  We  there¬ 
fore  consider  only  the  requirements  for  F(x)  in  digital  form  that  will  lead 
to  an  accurate  representation  of  I{v). 

Let  F(x)  be  measured  at  points  x  ,  x . ,  x7 ,  •  ■  ■  ,  then 

O  I  « 


qo 

I(v)  ar  \  "  F(x  )cos(2rrx  y)Ax 
/  -j  m  m  m 

m=0 


which  reduces  to. 


if  all  the  Ax  are  equal  and  x  is  at  x 
m  n  o 


0/ 


T - - — • — 

A  number  of  papers  describing  such  computers  were  given  at  the 
"Colloque  Sur  Les  Methods s  Nouvelles  de  Spectroscopic  Instrumentals,  " 
April  25-29,  1966,  Bellevue,  France.  Many  of  these  papers  are  scheduled 
for  publication  in  the  Journal  de  Physique  early  in  1967. 

2 

This  is  very  convenient  for  computer  programming  and  for  speed  in  the 
computer  calculation.  See  Section  VI, 
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(5) 


N 

I(v)  =  Ax[F(0)/2  + 

m=  1 


F(mAx)cos(2irmAxv)] 


Now  it  is  necessary  to  see  how  closely  this  summation  will  approximate 
the  true  I(v).  This  has  been  discussed  by  Strong  and  Vanasse  (Ref.  12). 

Making  use  of  sampling  theorems  from  information  theory,  they  obtain  if 

Ax  ~  TT 
c 

where  v  is  the  highest  wavenumber  having  non-zero  energy  in  the  spectrum, 
c  3 

then  l(v)  is  completely  determined. 

This  relation  is  quite  important  as  Ax  should  be  as  large  as  possible  in 
order  to  allow  more  time  to  measure  each  F(mAx).  We  anticipate  the  dis¬ 
cussion  of  Section  V  by  pointing  out  it  is  better  to  spend  more  time  measuring 
the  F(mAx)  with  the  largest  allowed  Ax  than  it  is  to  measure  a  larger  number 
of  the  F(mAx)  for  a  shorter  time  each  using  a  smaller  Ax.  What  happens 
when  Ax  is  chosen  too  large  is  discussed  in  Section  V.  It  ’S  also  apparent 
that  if  the  spectral  range  of  interest  extends  to  v  ,  then  vc  should  be  chosen 

(by  filtering,  etc.  )  as  close  to  v  as  is  compatible  with  good  S/N. 

ft 

Now,  what  is  the  effect  on  I(v)  when  the  sum  (or  integral)  is  not  carried 
out  to  infinity?  An  infinitely  long  (analog  or  digital)  interferogram  is  a  physical 


This  is  a  sufficient,  but  not  necessary,  condition.  Conr.es  (Ref.  6) 
has  shown  that  if  the  spectral  range  of  v  being  passed  by  the  interferometer 
lies  between  v.  and  v,,  it  is  possible  to  choose  Ax  =  [  2(  v,  -  v  j)]  or  some  - 
what  smaller,  out  certainly  >  (2  v  )“  .  Also  if  the  double-^beam  differencing 
technique  is  used  (see  Section  III)  and  if  the  highest  observed  absorption  is 
at  v  ,  then  one  can  choose  Ax  =  (2  v  )**  (really  a  bit  smaller  since  the  feature 
will  “In  general  have  finite  width)  eveti  though  there  is  energy  passed  by  the 
instrument  at  higher  wavenumber. 
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impossibility.  There  are  two  important  questions  to  consider:  First, 
in  the  absence  of  noise,  what  is  the  effect  of  the  truncation  of  the  interfero- 
gram,  i.  e.  ,  taking  the  integration  out  to  xmax  =  or  the  sum  out  to 
m  =  l/^x?  This  iB  answered  later  in  this  section  where  the  scanning  function 
of  the  two-beam  interferometer  is  discussed.  The  second  question  concerns 
the  effect  of  noise.  This  is  discussed  in  Section  V. 

Except  to  the  practitioners  of  the  art  of  two-beam  interferometry,  raw 
interferometric  data,  i.  e.  ,  the  interferograms,  appear  strange  and  unsug- 
gestive.  We  shall  now  discuss  some  simple  (idealired)  examples  of  interfero¬ 
grams  and  their  interpretation.  This  material  provides  the  back®*  s';nd  for 
interpretation  of  the  experimentally  observed  interferograms  discussed  in 
Section  IV. 

The  simplest  broadband  l(v),  such  as  a  source  would  present,  is 
I(v)  =  ,  a  constant  up  to  v  ,  and  zero  beyond.  Substitution  into  Eq.  (2) 

yields 

Fi<x)  =  I0vc[l  +  (siny/y)] 


where  y  =  2trxv  . 

c 

A  more  realistic  source  function  would  be  one  with  zero  or  little 
energy  at  low  wavenumber  (because  of  theblackbody  radiation  law),  rising 
toward  larger  v,  then  falling  to  zero  at  (due  to  the  filtering  optics). 

A  function  meeting  this  requirement  is 


The  factor  6  is  included  so  that  both  1( v)  functions  considered  here  have  the 

same  total  intensity  between  0  cm  *  and  v  cm*'. 

’  c 


Substituting  into  Eq.(2),  we  ootain 


f;(x)  =  i 


V 

o  c 


s  in  y 


cos  y  +  r 


with  y  defined  as  before. 

For  small  values  of  x.  F'^(x)  and  Fi,(x)  fluctuate  rapidly,  but  the  ampli¬ 
tudes  are  damped  such  that  both  F'(x)  curves  approach  asymptotically  the 
horizontal  straight  line,  F'(x)  =  -0^c  •  The  F'^(x)  fails  off  as  (2-rrxv^)  whereas 
F^(x)  falls  off  more  sharply,  approximately  as  (2ttxwc)  ^ .  In  each  case, 
however,  after  only  moderate  x  has  been  reached,  the  F'(x)  curves  become 
essentially  straight  lines.  This  is  what  as  observed  experimentally.  The 
F'j(x)  and  Fi,(x)  are  shown  in  Fig.  1,  with  ^  =  100  cm  Thus  when  an 
absorbing  medium,  such  as  a  gas 
whose  molecules  have  a  dipole  mo¬ 
ment,  is  placed  in  the  infrared  beam, 
any  variation  of  the  F(x)  curve  from 
a  straight  line  (for  intermediate  and 
large  values  of  x)  may  be  attributed 
to  the  absorbing  medium 


Now  let  us  consider  what  sort 
of  ante rferogram  a  single  absorption 
line  with  Lorentz  line  shape  would 
look  like  Let 


I(v)  -  IQMe 


•Y{v) 


whe  re 


Y(v> 


A  it 


o 


Fig.  1.  Interferograms  for 

"source -like "  I(v).  Line  curve 
Uv)  ~  i>lQ  [(‘'/v'c)  - 
dotted  curve  I(v)  -  IQ,  a  constant 
F’ (x)  is  plotted  in  units  of 
lo  vc(vc  =  100  cm"1). 
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where  A  is  a  constant  (independent  of  v),  L  the  absorbing  path  length,  e 
the  linewidth  parameter  (which  for  weak  absorption  is  the  half-width  at  half¬ 
maximum  intensity),  w  the  wavenumber,  and  the  resonant  wavenumber. 
Then 

F'(a>  =  ^l^di'  +  cos  2wxvdw 


Expanding  the  exponential  in  a  power  series  we  obtain 


F'(x)  =  I IQ(v)dw  +  I Iq(v)  cos  Z-nxvdv  + 


r  oo 

z 


4-  dv 


r  00 

M  z 


cos  2*txvc.v 


The  first  and  second  integrals  are  merely  the  source  interferogram. 
As  discussed  previously,  this  yields  essentially  a  constant  contribution  to 
F'(x)  for  intermediate  and  large  values  of  x.  The  third  integral  is  a  con¬ 
stant  and  represents  the  intensity  removed  from  the  beam  as  a  result  of  ab¬ 
sorption.  The  fourth  integral  will  give  the  deviation  of  the  interferogram 
from  a  straight  Fne  at  intermediate  and  large  values  of  x  and  thus  infor¬ 
mation  about  the  absorbing  line. 

Now  in  order  to  avoid  overly  cumbersome  mathematical  expressions 
let  us  simplify  the  above  expression  by  neglecting  all  except  the  first  terms 
in  the  summation.  A  more  general  treatment  may  be  found  in  the  literature 
(Ref  22).  This  will  be  a  good  approximation  for  weakly  absorbing  lines  and 
will  yield  a  qualitative  picture  for  absorption  lines  of  similar  shape.  There- 
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cos  2-rrxy dv 


F'(x) 


where  F'  (x)  is  the  inte  rfe  rogram  of  the  source, 
s 

Now  [{v-  vq)^  +  f4"]  *  is  a  sharply  peaked  function  (as  long  as  v  -■-.() 
and  ha.;  larg®  values  in  the  region  v  =  v  .  Thus  small  error  will  be  made 
in  the  following  approximations;  (a)  set  Jq(v)  r  Iq ( w'0) »  aod  (b)  extend  the 
limits  of  integration  to  t_oc,  so 


F  (x)  -  F'  (x)  -  I  (v  )A d 
s  o  o 


F (Jx)  -  I  (v  /  AJ?tt  ( 1  +  e  cos  2ttx v  ) 


The  effect  of  the  single  absorption  line  is  given  by  the  second  term 
and  its  dependence  on  x  by  the  term  in  brackets,  as  shown  in  Fig  2.  The 
amplitude  of  this  factor,  i  a)A is  the  intensity  removed  from  the  beam 
by  the  absorbing  line,  and  in  general  it  is  quite  small  in  comparison  with 
the  total  incident  intensity.  Thus.  F’  (x)  will  be  oniy  slightly  modified  by  a 
single  absorption  line.  As  an  illustration, for  /-  15  cm.  poe400  Torr,  the 
FIR  absorption  spectrum  of  CO  between  15  and  1  15  cm  *  (—20  lines,  not 
just  one)  removes  less  than  5°^o  of  the  incident  intensity  This  means  that 
the  variation  of  F'(x)  about  its  average  value  (1/2)F'(P)  will  be  quite  small  (less 
than  5%  at  best)  and  according  to  the  above  treatment,  if  the  lines  have 
approximately  Lorentzian  shape,  will  decay  exponentially  Making  c  smaller, 
which  can  be  done  for  CO  bv  lowering  the  pressure,  will  make  the  argument 
of  the  damping  exponential  smaller,  but  of  course  the  integrated  intensity 
will  also  become  smaller  (A  or  pressure).  This  is  a  characteristic  oroblem 
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of  absorption  interferometry  and  is  called  the  "dynamic  range"  problem. 
How  it  has  been  solved  is  discussed  in  Section  III. 

Another  simple  example  is  the  interferogram  generated  by  a  spectrum 
which  consists  of  evenly  spaced  (in  cm  *)  lines.  The  theoretical  interfere- 
grams  can  be  obtained  by  summing  over  the  absorption  line  contribution  in 
Eq.  (7).  Let  the  v  ’s  be  designated  v.  (lowest),  v, ,  i'~,  •  •  •  ,  v  (highest)  and 
be  equally  spaced  such  that 


vl  =  *2 


u  =  v  -  v  - 

1  3  2 


=  v  -  u  ,  = 
m  rn- 1 


or  in  terms  of  X. 


X  ~  1  X  1  *  - 

\  ~  ‘I  X3  " 


m 


The  varying  part  of  the  interferogram  will  be  a  sum  of  contributions  such  as 
that  shown  in  Fig.  2.  When  x  =  X^  (=l/6v),  then  x  =  2X2  =  3X^  =  •  •  •  =  m^m» 
so  all  the  cosine  contributions  will  have  a  maximum.  Likewise,  when 
x  =  2Xj  (=2/6v),  then  also  x  -  4X2  =  6X3  =  •  •  *  =  2m\m,  and  again  all  the  co¬ 
sines  will  have  their  maximum  value.  The  same  situation  occurs  at 
x  -  3Xj,  4Xj,  •  •  •  .  Thus,  significant  variations  of  F‘(x)  about  its  average 

value  (1/2)  F’  (0)  are  expected  when  x  =  mX^  =  m/6v.  Such  features, 
called  signatures  in  interferometry,  have  been  observed  in  the  FIR  inter- 
ferograms  of  linear  (Ref.  23)  and  diatomic  (Ref.  24)  molecules.  This  will 
be  further  amplified  and  discussed  in  Section  IV. 

We  now  come  to  the  question  of  finite  or  truncated  interferograms. 

For  this  topic  it  i3  natural  to  inquire  what  is  the  scanning  function  of  a  two- 
beam  interferometer.  The  scanning  function  of  a  spectrometer  or  inter¬ 
ferometer  may  be  thought  of  as  the  output  spectrum  of  the  instrument  when 
a  monochromatic  line  is  the  input.  Thus 


and 


I(v)  =  6(vo) 


F(x)  = 


cos  2iTXvdv  =  cos  2trx^ 

o 


Now  suppose  we  FCT  F{x)  out  to  x  =  L,  rather  than  x--*oo.  Then 


IL(v)  =  L 


sin  2ttL(  v  ~  v  ) 
o' 


■2VL7w-Voy 


a  L  sine  2irL (v-  v  ) 


(8) 


and  this  is  the  scanning  function  of  a  two-beam  interferometer.  This  is 

shown  in  Fig.  3.  The  IT  (v)  does  have  its  largest  maxima  at  v  ,  but  its 

resemblence  to  6(vq)  is  not  too  good  for  finite  L.  Perhaps  most  annoying 

(to  the  spectres copist  at  any  rate)are  the  extraneous  sidelobes,  or  feet, 

whose  amplitudes  decay  as  ~(v  -  v  )”*.  The  two  sidelobes  adjacent  to  the 

central  maxima  are  not  small  and,  as  seen  in  Fig.  3,  are  about  (first  minima 

to  second  maximum)  34%  of  the  peak  value  at  /  =  .  The  width  of  the 

scanning  function  is  determined  by  L  and  the  first  minima  are  at  v  - 

»i.O.  7  15(L)”\  The  relative  amplitudes  of  the  sidelobes  do  not  depend  on  L. 

In  contrast,  the  scanning  function  for  a  conventional  grating  spectrometer  is 
2 

of  the  form  (sine)  which  has  much  smaller  sidelobes. 

We  are  now  in  a  position  to  discuss  resolution.  Resolution  is  a  some¬ 
what  arbitrary  criterion,  usually  stated  in  terms  of  how  well  two  equally 
intense  lines  appear  separated  in  an  output  spectrum.  The  Rayleigh  criterion 
for  a  diffraction  grating  spectrometer  is  that  two  lines  of  equal  intensity 

3 - 

To  within  a  multiplicative  factor  I^(v)-»6(  v  )  as  L-*®.  See  W.  Heitler, 
Quantum  Theory  of  Radiation  (Clarendon  Press,  Oxford,  1954),  Chapter  II, 
page  66. 
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Fig.  3 


.  The  scanning  function  of  a  two-beam  interferometer. 
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should  be  considered  to  be  just  resolved  when  the  principal  intensity  maxi- 
mum  of  one  coincides  with  first  intensity  minimum  of  the  other.  On  this 
basis  there  will  be  a  dip  of  intensity  of  18.  9%  between  the  two  lines.  Figure 
4  shows  two  such  lines  separated  byAvsO.  715(L)’*,  i.  e.  ,  applying  the  first 
part  of  the  Rayleigh  criterion  to  two  lines  as  would  be  seen  by  a  two-beam 
interferometer.  The  dip  in  intensity  between  the  two  lines  is  considerably 
less  than  18.  9%.  Strong  and  Vanasse  (Ref.  12)  have  introduced  the  modular 
resolution  criterion  and  it  is 


'^modular  =  L'“  I’1 

This,  of  course,  is  more  conservative  thanAv  =  0.715  (L)~\  and  two  lines 
separated  by  thisAv  are  shown  in  Fig.  5.  The  dip  in  intensity  between  the 
two  lines  is  greater  than  18.  9%.  Additional  reasons  for  preferring  the 
modular  resolution  criterion  are  given  later  in  Section  IV. 

The  sidelobe  problem  appears  to  be  serious,  especially  as  compared 
with  its  almost  negligible  importance  in  a  grating  spectrometer.  It  is  pos¬ 
sible  to  suppress  these  sidelobes  by  apodization.  Essentially,  apodization 
consists  of  multiplying  the  interferogram  function  by  a  monotonically  de¬ 
creasing  function  of  x  such  that  F(L)  =  0.  The  theory  and  practice  of  apodi¬ 
zation  has  been  discussed  by  Filler  (Ref.  25).  One  of  the  disadvantages  of 
apodization  is  lower  resolution. 

Indiscriminate  apodization  can  be,  and  usually  is,  useless.  It  is 
instructive  to  consider  what  effect  the  scanning  function  has  upon  a  line  of 
finite  width  (lines  with  zero  width  are  quite  unphysical).  One  can  show  by 
convolution  theory  that  a  line  whose  shape  is  given  by  f(v)  will  appear  in  the 
output  spectrum  as  f(v)*S(v),  i.e.  ,  the  convolution  of  f(v)  with  S(v)  ar  the 
scanning  function.  We  shall  again  consider  a  line  with  a  Lorentz  shape. 


I(v)  = 


A€ 

- 2 — 7 

(w  -  vQ)  +* 
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Fig.  4.  Two  6  function  line#  separated  by  0.  715  (L)"1  as  would  be 

seen  by  a  two-beam  interferometer  which  has  been 

scanned  out  to  x  =  L. 

max 
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F(x)  =  irAe  cos  2rrx0 

Now  if  F(x)  is  FCT  (or  I(v)  is  convolved  with  S(v))  out  to  x  =  L,  we  obtain 


I  H  =  IM  !  1  -  e-2"Le 

L  I 

When  L-*oo,  it  is  obvious  I — * I(v')  as  it  should.  Also  to  within  a  mul¬ 
tiplicative  constant,  l(v)—  6(v  )  as  €  -•  0  and  I ,  (v )  — •  sine  2ttL{  v  -  v  )  as  €  -*  0. 
For  finite  L,  the  exponential  term  and  the  term  it  multiplies  describes  the 
distortion  of  the  line  from  its  true  shape.  Figures  6,  7,  and  8  are  plots  of 
Eq.  (10)  for  L£  =  oo,  3/4,  and  1/2.  On  the  scale  of  these  plots,  the  case 
L<  -  1  would  appear  identical  to  LC  =  oo.  In  general,  the  maximum  intensity 

of  IT  (v)  =  I.  (v  )  differs  from  I(v  )  as 
L  L  o  '  o 

IL(vo)  =  I(»0)  (1  -  «-2,rL<) 

(On  the  graphs,  A  is  chosen  such  that  I(vo)  =  1.  )  For  L<  =  3/4,  the  feet  are 

just  becoming  discernible  and  the  maximum  intensity,  IjJvq)  differs  from 

I(vo)  by~0.  8%.  For  L<  =  1/2,  the  eidelobes  are  now  quite  apparent  and 

I.  (v  )  is  ~4.  6%  weaker  than  I(v  ). 

L  o  o 

Let  us  define 


v  - 


cos  2ttL(  v  - 


«o>  - 


sin  2trL(w 


*o> 


(10) 


dl(v)  =  I(v)  -  IL(v) 


and  furthermore  let 

Vj  =  the  value  of  v  where  4I(v)  has  its  first  minima,  and 

=  the  value  of  v  where  4l(v)  has  its  second  maxima  (the  first 

maxima  is  of  course  at  v  =  v  ). 

o 
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Now  we  define 


IF  *AIfv2)  -  4l(Vl) 

so  that  Ip  measures  the  maximum  sidelobe  intensity.  A  plot  of  IjJv0)/Ip 
versus  L«  is  shown  in  Fig.  9.  This  curve  rises  quite  steeply  as  a  function 
of  L>(  ,  i.  e.  ,  as  the  width  of  the  scanning  function  (f"werned  by  L)  approaches 
€,  the  sidelobe  intensity  decreases  very  rapidly.  As  an  example  of  its  use¬ 
fulness,  suppose  the  S/N  in  the  spectrum  is  100/1;  then  there  is  no  point 
in  making  Ip{vQ) /Ip  much  greater  than  100/ 1,  which  means  L£-»0.S1.  For 
L€~0.  81,  the  interferogram  e  ^Trx*  COs  2irxvQ  will  be  down  to  e of 
its  value  at  x  =  0,  which  is  down  by  a  factor  of  162.  The  use  of  this  curve 
when  noise  in  the  interferogram  is  considered  is  discussed  in  Section  V. 


H)-4«  »9-Z*  fg  -«  >•  rgl«  »|  +  Z«  H+J« 


Fig.  6.  The  Lorentt  line  shape. 
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Fig.  9.  Plot  of  Ij^  ( vQ)/Ip  versus  L*  . 


III.  DESCRIPTION  OF  THE  INSTRUMENT 


A.  OPTICS 

The  optical  diagram  of  the  lamellar  grating  interferometer  is  shown 
in  Fig.  10.  The  source  of  radiation  is  a  high-pressure,  85-W  fused  quartz. 


Fig.  10,  Lamellar  grating  interferometer  optical  diagram. 

The  foci  of  the  ellipsoids  and  M2  are  between 
M5  and  M^  (sample  beam)  and  at  M3  (reference 
beam),  respectively. 


mercury  arc  lamp  (G.  E.  H85A3)  with  the  outer  glass  envelope  removed. 

The  arc  is  mounted  in  a  water-cooled  housing  with  two  cutouts  through  » 

which  the  reference  and  sample  beams  emerge.  The  two  beams  are  caught 

and  focused  by  two  ellipsoidal  mirrors  M^  and  M.,.  The  upper  or  sample 

beam  is  folded  by  two  flats  M^  and  Mg  and  then  passed  through  the  sample 

chamber.  Mirror  Mg  is  a  polished  aluminum  flat  to  match  the  reflectance 

of  the  aluminum  chopper  wheel  Mg.  The  beam  comes  to  a  focus  in  the 

middle  of  the  sample  chamber  allowing  the  use  of  minimum  size  samples. 

The  focus  is  a  circle  of  0.  95  cm  diam.  The  lower  or  reference  beam  passes 
under  the  sample  chamber  and  comes  to  a  focus  on  the  flat  Mg.  Flat  Mg 
is  located  on  a  six-position  slide  assembly  which  allows  changing  between 
the  mirror  and  any  of  five  2.  5 -cm  diam  samples  for  reflection  studies 
from  outside  the  tank  without  breaking  the  vacuum  in  the  chamber.  The 
flat  Mg  folds  up  the  reference  beam  perpendicular  to  the  sample  beam.  The 
two  beams  are  then  recombined  by  the  polished  aluminum  semicircular 
chopper  mirror  Mg  which  rotates  at  19  cps.  To  match  precisely  the  intensi¬ 
ties  of  the  two  beams,  motorized  iris  diaphragms  have  been  placed  in 
equivalent  positions  in  both  beams  approximately  13  cm  before  the  first 
focus. 

Mirror  M^  is  a  33- cm  diam  spherical  mirror  which  refocuses  the 
combined  beams.  A  limiting  aperture,  which  reduces  the  angular  spread 
of  the  beam  by  a  factor  of  2,  can  be  moved  into  place  in  front  of  M^  from 
outside  the  tank.  This  feature  is  of  great  use  in  studies  of  solid  samples. 

After  striking  Mg  and  M^,  the  beam  is  refocused  on  the  circular  entrance 
aperture  (not  shown  in  the  figure)  of  the  interferometer  proper.  The  size  of 
the  aperture  is  determined  by  the  high-frequency  radiation  cutoff  of  the 
system  (Ref.  12). 

The  interferometer  optics  consist  of  two  45.  7 -cm  diam  spherical 
mirrors  and  the  lamellar  grating.  The  two  spherical  mirrors  are  arrayed  , 

in  a  Czerny -Turner  configuration  which  allov/s  the  grating  to  be  used  in 
parallel  light  in  the  zeroth  order.  1 
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The  lamellar  grating  consists  of  two  sets  of  interleaving  plane  mirrors, 
as  shown  in  Fig.  11.  There  are  24  facets  in  each  set,  each  facet  being 
0.  635  cm  wide  and  30.  5  cm  long,  giving  a  total  area  of  30.  5  X  30.  5  cm  for 

5 

the  grating.  The  grating  is  made  of  Pyrex.  The  individual  facets  of  each 
grating  are  flat  and  coplanar  to  two  fringes  of  the  5461  X  line  of  Hg,  except 
for  the  outside  facets  of  each  grating,  which  deviat.  by  about  four  fringes. 
These  four  facets  have  been  blanked  off  to  eliminate  distortions  in  the  inter- 
ferograms  caused  by  their  optical  imperfections. 

After  the  optical  path  difference  has  been  introduced  by  the  grating, 
the  beam  is  reflected  to  the  second  Czerny-Turner  mirror,  which  after 
being  folded  by  the  diagonal  flat  (scatter  grating),  refocuses  it  on  the 
exit  aperture  mounted  on  the  end  of  a  light  pipe.  Flat  is  adjustable 
from  outside  the  tank  to  permit  the  final  alignment  to  be  made  in  vacuum. 

This  is  desirable  because  nearly  all  the  FIR  energy  is  absorbed  by  water 
vapor  in  the  air.  The  size  of  the  exit  aperture  is  the  same  as  the  entrance 
aperture.  The  circular  light  pipe,  65  cm  in  length,  takes  the  beam  to 
the  wall  of  the  vacuum  tank  where  a  gallium-doped  germanium  bolometer 
(Texas  Instruments)  is  mounted.  The  detector  operates  at  4.  2°K  and  has 
a  sensitive  area  of  5  X  5  mm. 

One  of  the  severe  problems  of  the  FIR  is  radiation  purity.  The 
mercury  arc  lamp  emits  orders  of  magnitude  more  energy  in  the  ultraviolet 
(UV),  visible,  and  near-IR.  Eliminating  this  high-frequency  radiation  re¬ 
quires  an  elaborate  array  of  transmission  and  reflection  filters.  However, 
in  interferometry  the  problem  is  not  nearly  as  difficult  as  in  conventional 
dispersive  spectroscopy.  Filtering  is  accomplished  in  this  instrument  by 
the  following  means:  At  positions  and  M^,  coarse  (<20  lines/mm) 
gratings  are  placed,  which  provide  the  high-frequency  cutoff  filtering  of  the 
source  radiation.  Wavelengths  shorter  than  approximately  1.5  times  the 

Manufactured  by  Davidson  Optronics  of  West  Covina,  Calif. 
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Fig.  11,  Front  view  of  the  lamellar  grating.  At  the  right  is  the  chopper 
assembly,  the  back  side  of  Mg  and  the  six-position  slide 
assembly  on  which  Mg  is  mounted.  Mirror  Mjq  is  just  below 
the  grating;  Mjg,  which  is  replaced  by  a  scatter  grating  after 
alignment,  is  above  the  grating.  Also  shown  (in  front  of  Mg) 
are  the  filter  wheel  and  variable  diaphragm. 
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blaze  wavelength  of  the  grating  are  thrown  out  of  the  beam  at  the  blaze 
angle  of  the  grating.  To  longer  wavelengths,  the  gratings  appear  as  plane 
mirrors.  By  changing  these  gratings,  the  high-frequency  cutoff  can  be 
changed.  Two  identical  gratings  in  series  sharpen  the  cutoff.  Various 
sets  of  gratings  are  available  to  vary  the  high-frequency  cutoff.  In  addition 
to  filtering  by  the  two  scatter  plate  gratings,  there  is  filtering  by  black 
polyethylene  and  crystalline  quartz.  A  sheet  of  black  polyethylene  is 
located  between  the  source  and  chopper  in  front  of  the  sample  area  on 
the  grating  side  of  the  source  housing  after  mirrors  Mg  and  M.,.  This 
cuts  down  the  UV  and  visible  light  entering  the  sample  chamber  and  helps 
to  eliminate  the  sometimes  deleterious  effects  of  UV  light  on  samples 
under  study.  If  gas  samples  are  to  be  studied,  they  must  be  contained  in 
a  conventional  gas  cell  that  has  windows  offering  two  more  filtering  positions. 
Wedged  crystalline  quartz  windows  are  used.  To  compensate  for  the  ab¬ 
sorption  in  the  sample  beam,  two  nearly  identical  pieces  of  crystalline 
quartz  are  placed  in  the  reference  beam.  These  are  located  between  Mg 
and  Mg  on  a  four-position  filter  wheel  controlled  from  outside  the  tank. 

In  the  other  three  positions  of  the  filter  wheel  are  an  opaque  shutter  for 
blanking  off  the  reference  beam,  an  open  hole  for  studies  where  the  filtering 
of  crystalline  quartz  is  not  wanted,  and  one  additional  space  where  any  other 
compensation  filter  can  be  placed  (e.  g.  ,  polyethylene).  The  crystalline 
quartz  absorbs  the  intermediate  IR  radiation  of  the  source.  The  quartz 
pieces  in  the  reference  beam  are  slightly  thinner  --  by  a  few  hundredths  of 
a  millimeter  --  than  those  on  the  sample  cells.  The  four-position  filter 
wheel  can  be  rotated  about  an  axis  perpendicular  to  the  optic  axis.  This 
allows  the  effective  thickness  of  the  filters  to  be  varied  to  obtain  a  precise 
equivalence  of  both  beams,  which  is  necessary  to  ensure  complete  balance 
of  the  spectral  transmission  of  both  beams.  Slight  deviations  in  balance 
lead  to  anomalies  in  intensity  for  the  region  of  zero  optical  path  difference. 

Figure  11  shows  the  lamellar  grating  and  mirrors  Mg,  M^,  Mg,  M10 
and  M  j g.  Figure  12  shows  a  side  view  of  the  interferometer.  Figure  13 
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Fig.  12.  Side  view  of  the  interferometer.  At  the  left  is  the  source 
housing  and  optics.  The  side  view  of  the  chopper  is  in  the 
center  of  the  figure.  A  "swing-away"  limiting  diaphragm 
can  be  seen  in  front  of  mirror  M7.  The  interferometer  is 
on  a  dolly  with  rails  matching  those  in  the  vacuum  tank 
(extreme  left). 


gives  an  overall  view  of  the  interferometer.  The  sample  chamber  slides 
into  position  between  the  source  housing  and  the  chopper  wheel.  The  sample 
chamber  has  been  put  in  place  on  the  left  side  and  the  light  pipe  attached  to 
the  bolometer  on  the  right  side.  The  bolometer  is  located  in  an  elbow  just 
outside  the  hole  from  which  the  light  pipe  is  protruding.  All  electrical  and 
cooling  water  connections  are  made  by  vacuum-tight  connectors  in  the  bulk¬ 
head  at  the  lower  left. 


Fig.  13.  The  interferometer  in  place  in  the  vacuum  tank.  The  light 
pipe  can  be  seen  protruding  from  the  right -hand  side  of 
the  chamber  wall.  Ine  sample  chamber  is  in  place  on  the 
left  side  (appears  black).  The  rails  upon  which  the  instru¬ 
ment  rides  are  at  the  bottom  and  the  electrical  cooling 
water  connections  at  the  bottom  left.  The  diameter  of  the 
vacuum  tank  is  1.68  m. 
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B.  MECHANICS 


The  interferometer  is  mounted  on  a  heavy  steel  box  frame  which  is 
kinematically  supported  in  a  vacuum  chamber.  Thus  deflections  in  the 
vacuum  chamber  walls  will  not  alter  the  alignment  of  the  optics.  ^  Evacua¬ 
tion  is  necessary  to  eliminate  the  strong  water  vapor  absorption  in  the 
FIR.  A  10-in.  vacuum  pump  system  can  evacuate  the  chamber  with  all  the 

_5 

apparatus  inside  to  a  pressure  of  10  Torr  in  a  few  hours.  Liquid- nitrogen 
trapping  helps  to  remove  the  last  traces  of  water  vapor.  Figure  14  shows 
the  vacuum  tank  and  the  electronics  console. 


F*g.  14,  Electronic  console  and  vacuum  tank 


°The  interferometer  frame  and  mirror  mounting  were  manufactured 
by  American  Astrophysics,  Monrovia,  Calif.  The  vacuum  chamber  was 
manufactured  by  Advance  Tank  and  Manufacturing  Company,  Gardena,  Calif. 
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The  sample  chamber  can  be  isolated  from  the  rest  of  the  vacuum 
chamber  by  gate  valves,  which  also  function  as  ports  through  which  the 
sample  beam  passes.  With  the  gate  valves  closed,  the  sample  chamber 
can  be  brought  up  to  atmospheric  pressure  without  impairing  the  vacuum 
in  the  main  chamber  for  rapid  change  of  samples.  By  closing  one  of  the 
gate  valves,  single-beam  interferograms  of  the  reference  beam  can  be  ob¬ 
tained. 

Access  ports  have  been  provided  in  the  chamber  walls  to  allow  minor 
adjustments  to  be  made  from  outside  the  chamber  without  rolling  out  the 
entire  interferometer.  The  entire  interferometer  can  be  lowered  down  on 
wheels  which  ride  on  rails  for  servicing  outside  the  chamber.  A  matching 
set  of  rails  on  a  dolly  mate  with  these  rails,  ard  the  interferometer  can 
be  rolled  out  onto  the  dolly. 

The  entire  vacuum  chamber  is  vibrationally  isolated  from  the  floor 
by  air  suspension  mounts  under  each  of  the  six  legs.  These  mounts  isolate 
the  chamber  from  vibration  and  shock  down  to  a  few  cycles  per  second, 
eliminating  spurious  modulation  due  to  external  vibration. 

The  heart  of  the  interferometer  is  the  grating  and  its  drive  system. 

The  lower  set  of  facets  of  the  grating  is  mounted  on  a  movable  platform 
suspended  from  a  stationary  platform  holding  the  upper  set  of  facets  by 
three  flexure  hinges  similar  to  those  used  by  Strong  (Ref.  13).  These  three 
hinges  theoretically  allow  the  lower  platform  to  move  in  a  perfectly  parallel 
fashion  with  respect  to  the  upper  platform.  Unfortunately,  torsional  oscil¬ 
lations  and  pitch,  roll,  and  yaw  errors  of  the  order  of  1  arc-deg  necessi¬ 
tated  further  constraint.  A  flat,  parallel-aided  steel  bar  was  installed 
below  the  lower  platform.  Two  sets  of  ball  race  bearings  were  attached 
to  the  lower  platform  contacting  the  bar  on  its  two  vertical  surfaces  in 
order  to  eliminate  torsional  oscillations  and  reduce  the  pitch  and  yaw  errors 
to  less  than  i.  1  arc  sec  ov-ec  83  mm  of  travel.  The  errors  in  pitch  and  yaw 
are  smoothly  varying  functions  of  the  travel  and  are  the  worst  at  the  extreme 
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end  of  the  travel.  The  roll  error,  net  important  optically,  v.-a *  reduced 
to  below  the  point  where  scraping  of  the  facets  was  a  danger.  Sines  the 
platform  rises  through  an  arc  about  0.  i  in.  during  the  travel,  no  restraint 
could  be  applied  to  this  motion 

The  grating  platform  it  driven  by  a  stepping  motor  which  has  200 
steps  per  revolution.  The  use  of  a  stepping  motor  eliminates  coasting  and 
allows  fine  control  of  the  drive  by  controlling  the  drive  pulse  rate  to  the 
motor.  The  stepping  motor  is  geared  down  by  very  low  back?-  - h  gears. 

The  drive  is  transmitted  to  the  grating  platform  by  mean*  of  a  prtoiaio:. 
lead  screw  and  a  pin  riding  on  the  driver  nut.  The  entire  piatxorm  nre 
loaded  to  ensure  that  the  driver  nut  always  rides  on  the  *ame  side  of  the 
lead  screw.  The  drive  system  is  capable  of  moving  the  grating  platform  in 
increments  less  than  0.  5p. 

The  relative  position  of  the  two  sets  of  facets  is  me  asured  by  two 
separate  systems,  coarse  and  fine.  The  coarse  system  is  capable  of 
measuring  the  absolute  position  of  the  facet*  to  within  1  mm  over  the  entire 
path.  The  fine  system  can  measure  the  absolute  position  of  the  facets  to 
±1.  5p  over  a  2  mm  cycle.  The  two  complementary  systems  together  can 
measure  the  absolute  position  of  the  facets  to  ±  1. 5p  over  the  entire  path. 

The  position  of  zero  optical  path  difference  can  be  determined  u,  an  accuracy 
of  better  than  £  l.  Sp.  Relative  to  this  position  there  is  a  constant  precision 
of  better  than  £l.  5p  for  the  location  of  any  point  over  the  entire  range  of 
travel.  Unfortunately,  the  error  varies  its  a  periodic  manner,  which  leads 
to  apparent  "ghosts"  in  the  spectra  obtained.  This  is  discussed  in  Section 
V.  A  change  in  the  method  of  measuring  grating  position  is  being  planned, 
employing  •  laser  interferometer. 

The  coarse  system  uses  a  synchro  to  determine  the  position  of  the 
grating  platform.  Approxi  tely  one  revolution  of  the  synchro  corresponds 
to  the  full  travel  of  the  grating  platform  . 


The  fine  system  uses  a  linear  metric  Inductosyn  scale  (Farrand 
Controls,  Inc.  ,  Valhalla,  New  York)  to  determine  the  position  of  the  grating 
platform.  This  system  similar  to  thos'  used  on  automatic  milling  machines, 
measures  position  by  means  of  the  inductile  coupling  between  a  stationary 
set  and  moving  set  of  linear  pole  patterns.  The  moving  set,  supported  on  a 
glass  blank,  is  mounted  on  the  side  of  the  grating  platform.  The  stationary 
set  is  mounted  on  the  grating  platform  immediately  opposite  the  moving  set 
on  the  side  support  member.  The  stationary  set  is  250  mm  long,  allowing 
its  use  over  the  full  travel  of  the  grating  platform.  Figure  15  shows  a  back 
view  of  the  lamellar  grating.  The  Inductosyn  scales  can  be  se  n  at  the  right 
and  the  drive  system  underneath  the  grating  platform. 

The  grating  platform  has  a  maximum  travel  of  10.  0  err.,  of  which 
5.  3  cm  are  useful,  0.  3  cm  on  one  side  of  the  coplanar  position  of  the  grating 
and  8,  0  cm  on  the  other  side,  A  travel  of  8.  0  cm  is  equivalent  to  an  optical 
path  difference  of  16.  0  on,  giving  a  modular  resolution  of  0.  063  cm-1.  The 
intensity  passed  by  the  interferometer  as  a  function  of  the  optical  path  dif¬ 
ference  generated  by  the  interferometer,  known  as  the  interferogram,  is, 
in  theory,  symmetric  about  the  coplanar  position.  In  practice  we  have 
found  this  will  not  be  true  unless  the  incident  parallel  beam  is  strictly  per¬ 
pendicular  *o  the  grating  facets  in  the  a  *;*,.«***  perpendicular  to  the  grating 
grooves.  The  incident  beam  does,  of  course,  deviate  from  perpendicularity 
in  the  a&imuth  parallel  to  the  grating  grooves  because  of  the  Cxerny-Tui net 
optical  System.  This  means  that  the  optical  path  difference  generated  is 
not  2d  (d  =  geometrical  separation  of  two  sets  of  facets),  but  rather  2d  cos  a, 
where  a  is  the  angle  between  incident  beam  and  the  normal  to  the  facet  faces 
The  angle  is  approximately  9  deg  and  it  is  measured  with  a  theodolite  to  an 
<ccuracy  of  $.15  arc  sec,  which  is  sufficiently  accurate  to  eliminate  it  as 
s  source  of  systematic  error  in  the  measurement  of  optical  path  difference 
The  normals  of  the  two  set*  of  gratings  are  made  parallel  by  adjusting  the 
gratings  in  their  mounts  while  monitoring  the  adjustment  with  a  comparison 
autocollimator  capable  of  detecting  deviation  from  parallelism  of  leas  than 


Fig.  15.  Back  view  of  lamellar  grating.  The  flexure  hinge  supporting 
the  movable  platform  is  visible  in  the  middle  of  the  figure. 

The  e  are  two  more  flexure  hinges  in  the  front  of  the  platform. 
On  he  right  are  the  linear  inductosyn  plates.  The  drive  system 
can  be  seen  be'ow  the  movable  platform. 
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1  arc  sec.  The  entire  usable  travel  of  the  grating  (8.  3  cm)  has  been 
monitored  and  is  rechecked  after  any  adjustment  with  this  comparison  auto- 
collimatcr.  With  these  precautions,  very  symmetric  interferograms  are 
obtained. 

C.  ELECTRONICS  AND  DATA  RECORDING  SYSTEM 

The  interferogram  is  not  recorded  continuously,  but  is  sampled  at 
equally  spaced  increments.  The  system  is  designed  to  sample  at  multiples 
of  1 0 jj.  in  optical  path  difference.  The  actual  sampling  interval  is  dictated 
by  the  highest  frequency  present  and  has  been  discussed  previously. 

The  grating  is  allowed  to  sample  each  point  for  2  to  8  sec  to  im¬ 
prove  the  S/N.  During  this  time  the  output  of  the  detector,  suitably  ampli¬ 
fied  and  synchronously  detected,  is  fed  into  a  voltage-to-frequency  converter 
that  converts  the  voltage  output  linearly  into  a  frequency  which  is  then 
counted  for  the  gate  or  sampling  time.  At  the  end  of  this  time,  the  number 
of  counts  (0  to  9999)  accumulated  in  the  counter  and  the  number  of  the  point 
are  punched  on  paper  tape  and  printed  on  a  digital  printer.  The  grating 
then  advances  to  the  next  point,  and  the  process  is  repeated.  This  process 
is  continued  until  the  entire  interferogram  has  been  punched  out  on  the 
paper  tape.  This  tape  output  is  then  processed  by  the  data  reduction  center 
and  the  Fourier  cosine  transformation  is  computed  yielding  the  frequency 
spectrum.  The  digital  printer  output  can  be  used  to  check  for  errors  in  the 
input  data  to  the  computer.  Unfortunately,  the  use  of  punched  paper  tape 
has  not  been  entirely  satisfactory.  We  are  in  the  process  of  converting  to 
punched  cards. 

The  grating  drive  system  has  been  designed  to  operate  in  several 
different  modes.  These  range  from  a  completely  manual  mode  to  a  com¬ 
pletely  automatic  sampling  mode.  Other  modes  of  the  grating  drive  system 
arc:  the  automatic  mode,  in  which  the  grating  will  proceed  to  and  remain 
at  any  position  dialed  into  the  controls;  and  the  manual  mode,  in  which  the 
grating  can  be  moved  under  direct  control  of  the  operator  and  thus  is  useful 
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for  survey  work 

Figure  16  is  a  block  diagram  of  the  grating  drive  and  data  recording 
systems.  The  strip  chart  record  is  used  primarily  for  visual  indication  of 
performance  during  a  run.  The  oscilloscope  is  used  in  a  similar  fashion. 

The  interferogram  can  be  recorded  either  single-beam  or  double¬ 
beam.  In  double-beam  operation,  the  outputs  of  the  reference  and  sample 
beams  are  differenced  automatically  by  the  broadband  ac  preamplifier 
following  the  detector.  This  difference  signal  is  recorded.  In  single-beam 
operation,  one  of  the  two  beams  is  blanked  off.  As  has  been  discussed 
previously, there  is  a  problem  of  "dynamic  range"  in  absorption  interfero¬ 
metry.  Double -beaming  circumvents  this  problem.  In  double -beaming 
two  beams  are  used,  one  of  which  passes  through  the  sample  cell  (sample 
beam)  and  the  other  through  a  matched  reference  cell  (reference  beam); 
then  the  sample  beam  intensity  (or  signal)  is  subtracted  from  the  reference 
beam  intensity  (or  signal).  When  this  is  done,  most  of  the  information  about 
the  spectral  distribution  of  energy  from  the  source  is  lost,  but  such  infor¬ 
mation  may  be  recovered  by  a  prior  or  subsequent  recording  of  single-beam 
interferograms  of  the  source  alone.  In  practice,  this  is  not  difficult  or 
inconvenient  as  the  single -beam  source  interferogram  need  run  out  only  to 
a  3-  to  5-mm  optical  path  difference  and  can  be  obtained  in  a  time  interval 
short  compared  with  the  higher  resolution  double -beam  interferograms  where 
the  optical  path  difference  typically  runs  out  to  100  to  160  mm. 

We  write  the  appropriate  equations  for  the  intensity  passed  by  the 
interferometer  for  the  reference  and  sample  beams: 


1  ,  c 

reference  beam:  F*(x)  =  *  j  Ir(v')cos2irxv'dv 
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Fig.  16.  Block  diagram  of  lamellar  grating  interferometric  system. 
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"C 

sample  beam:  F'(x)  =  iF'(0)  +  I  I  (v)cos2Trxvdv 

s  c  a  j  s 

o 

where  all  the  symbols  have  their  previous  meaning,  and  the  subscripts  r 
and  s  refer  to  the  reference  and  sample  beam,  respectively.  Differencing 
these  equations  and  transposing  the  constant  term  to  the  left  side,  we  ob¬ 
tain  the  double -beam  interferogram  function,  F  (x): 

r-s 

vc 

Fr  g(x)  s  F'  (x)  -  F;(x)  -  \  [F^(0)  -  F'(0)]  =  j  [Ir(v)  -  I#(v)]  cos2Trxvdv  (11) 


With  matched  beams  and  no  sample,  the  detector  will  see  essentially  a  dc 
signal  that  will  not  be  amplified  by  the  broadband  ac  preamplifier.  In 
practice,  because  of  the  finite  dimensions  of  the  chopper  (chopping  is  not 
done  at  a  focus,  see  Fig.  10),  there  is  ripple  in  the  signal  at  the  chopping 
frequency;  however,  the  ripple  is  nearly  90  deg  out  of  phase  with  any  true 
signal  and  is  largely  filtered  out  by  the  phase  lock  amplifier. 

To  balance  the  beam,  it  is  necessary  to  have  matched  spectral  dis¬ 
tribution  of  energy  in  the  two  beams.  The  total  intensity  is  matched  by  a 
motorized  iris  in  each  beam  (set  back  ~18  cm  from  the  focus).  Making  the 
beams  spectrally  equivalent  is  more  difficult  to  achieve  since  there  must 
be  equal  filtering  in  each  beam.  For  gas  spectroscopy,  a  cell  is  required 
in  the  sample  beam  (usually  quartz  windows  are  used)  so  that  an  equivalent 
set  of  windows  must  be  inserted  in  the  reference  beam.  As  mentioned  pre¬ 
viously,  they  may  be  rotated  through  small  angles  to  provide  a  variable 
thickness  of  quartz.  This  is  necessary  because  of  wedging  of  windows  (to 
prevent  a  troublesome  channel  spectrum)  and  optical  inhomogeneity  of  the 
windows.  In  practice,  perfect  balance  is  not  obtained,  but  with  the  above- 
mentioned  modifications,  near  perfect  balance  can  be  obtained  everywhere 
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Fig.  17.  Single-beam  (lower)  and 
double- beam  (upper)  inter¬ 
ferometric  features  of  the 
pure  rotational  spectrum  of 
carbon  monoxide. 


Fig.  18.  Single-beam  (lower)  and 
double -beam  (upper)  inter¬ 
ferometric  features  of  the 
channel  spectrum  of  a  plane 
parallel  window  of  crystalline 
quartz. 


except  in  the  region  of  x  =  0.  The  unbalance  there  is~20  times  the  noise 
level  of  the  instrument. 

The  Fourier  cosine  transform  of  the  left  side  of  Eq.  (11)  will  yield 
[lr(v)  -  IgM]  .  the  difference  between  the  reference  beam  spectrum  and 
the  sample  beam  spectrum.  Thus  an  absorption  spectrum  appears  as 
though  it  were  an  emission  spectrum.  As  discussed  above,  the  source 
(reference  beam)  spectrum  can  be  easily  obtained  from  single-beam  runs. 

The  advantages  of  double-beam  differencing  over  single-beam 
operation  can  b»  discussed  in  terms  of  Figs.  17  and  18,  which  compare 
single-beam  and  double-beam  interferograms  of  carbon  monoxide  and 
quartz,  all  obtained  with  the  central  maximum  at  full  ccale.  As  can  be 
seen,  the  "dynamic  range"  problem  has  been  greatly  alleviated  by  double¬ 
beam  differencing, 
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IV.  RESULTS 


A  HIGH  RESOLUTION  SPECTROSCOPY  OF  GASES 

1.  CARBON  MONOXIDE 

A  portion  of  the  double -beam  interferogram  of  the  pure  rotational 
spectrum  of  carbon  monoxide  (CO)  is  shown  in  Fig.  19.  The  absorption 
path  length  was  15  cm,  and  the  pressure  in  the  gas  cell  was  7  32  Torr.  Had 


this  interferogram  been  obtained 
single  beam,  keeping  the  central 
maxima  full  scale,  the  features  due 
to  the  pure  rotational  spectra  of  CO 
would  be  down  in  intensity  by  a  L  c- 
tor  of  ~15. 

Carbon  monoxide  is  diatomic 
and  therefore  its  pure  rotational  g 
spectrum  consists  of  almost  evenly  ^ 
spaced  lines  (in  cm  *).  As  dis-  1 

cussed  in  Section  II,  its  pure  rota-  5 

tional  interferogram  should  have 
signatures  spaced  at  equal  intervals 
in  x,  the  optical  path  difference. 

These  signatures  are  apparent  in 
Fig.  19,  and  the  reciprocal  of  their 
spacing  gives  an  average  value  of 
2B0  (the  separation  of  the  lines  in 
cm'*).  It  is  alsc  seen  that  the  sig¬ 
natures  are  distorting  with  larger 
optical  path  diffe  rence .  Qualitative  - 
ly,  the  first  large  minimum  is  losing 
amplitude  and  the  second  minimum 


OPTICAL  PATH  WFTtAWCt  (■«) 


Fig.  19.  The  pure  rotational 
double-beam  interferogram  of 
CO.  The  evenly  spaced  features 
are  characteristic  of  diatomic 
and  linear  molecules. 
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is  gaining,  which  is  due  to  the  fact  that  the  pure  rotational  spectrum  of  CO 
is  not  made  up  of  uniformly  spaced  lines,  but  rather  lines  that  converge  at 
higher  wavenumber  because  of  centrifugal  distortion.  This  leads  to  the 
distortion  of  the  signatures.  If  the  lines  had  been  diverging,  the  appropri¬ 
ate  lobes  would  have  lost  and  gained  amplitude  in  the  reverse  order.  Also, 
it  is  seen  that  signatures  decay  very  rapidly  in  amplitude,  over  and  above 
that  caused  by  the  distortion  of  the  signatures,  and  this  is  due  to  the  finite 
width  of  the  absorption  lines.  In  this  interferogram  ~25  signatures  were 
observable  before  they  were  lost  in  the  noise.  At  lower  pressures  one 
would  of  course  expect  to  see  more  signatures  (lines  are  narrower),  and  in 
fact  at  p  =  202  Torr  we  have  seen  the  6 1st  signature  at  x  *  16  cm.  Finally, 
there  is  an  extra  signature  at  x®  3.  6  mm.  That  is  caused  by  a  channel 
spectrum  generated  in  the  quartz  envelope  of  our  source.  It  produces  a 
cosine  modulation  in  the  spectrum.  We  have  found  that  by  proper  editing  of 
the  interferogram,  th'*  channel  spectrum  can.  be  eliminated.  This  problem 
is  discussed  further  in  Section  V. 

Figure  20  shows  the  result  of  performing  the  Fourier  cosine  trans¬ 
formation  on  the  CO  interferogram.  that  is, the  pure  rotational  spectrum  of 
CO  between  19  and  100  cm’*,  The  dashed  line  is  the  approximate  envelope 
of  the  source  spectrum  (absorption  goes  up)  and  would  be  the  place  to  which 
a  100%  absorbing  line  would  go.  The  variation  in  the  base  line  is  the  channel 
spectrum  mentioned  previously.  These  features  are  quite  broad  and  are  not 
mistaken  for  the  much  sharper  pure  rotational  lines  of  CO.  The  wave- 
numbers  for  the  pure  rotational  spectrum  of  CO  can  be  determined  quite 
accurately  from  microwave  (Ref.  26)  and  near-IR  (Ref.  27)  data;  therefore, 
the  measurement  of  this  spectrum  will  not  be  discussed  further,  although 
it  was  very  useful  for  calibration  purposes. 

It  is,  however,  important  to  look  at  the  line  shape.  As  is  well 
known  in  interferometry,  if  even  a  small  error  is  made  in  the  measurement 
of  optical  path  differences  (see  Section  V),  shifts  of  the  lines  and  serious 
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Fig.  20.  The  pure  rotational  spectrum  of  CO.  The  intensity 
scale  is  arbitrary,  but  the  dashed  line  indicates  the  approxi¬ 
mate  intensity  contour  of  the  source.  The  cosinusoidal  vari¬ 
ation  of  the  base  line,  especially  noticeable  at  low  wave- 
number,  Arises  from  the  channel  spectrum  due  to  the  lamp 
envelope. 


distortion  of  line  shapes  occur,  Figure  21  shows  a  typical  result  from  our 
interferometer.  Shown  is  the  J  =  16 -*”17  transition  of  CO.  The  open  circles 
are  the  intensity  points  as  given  by  the  computer  The  line  looks  and  is 
very  symmetric.  In  order  to  examine  the  line  shape  more  closely,  the  first 
derivative  of  this  line  shape  was  computed  (open  squares).  Ii  has  the  ex¬ 
pected  shape,  but  there  is  about  a  1%  difference  between  the  maxima  and 
minima.  If  the  line  had  been  perfectly  symmetric  these  would  have  been 
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equal  but  of  opposite  sign.  This 
excellent  symmetry  of  line  shape 
and  reasonably  good  first  deriva¬ 
tive  shape  have  encouraged  us  to 
pursue  line  shape  measurements, 
and  this  will  be  discussed  later 
in  this  section.  Also  if  the  line 
were  perfectly  symmetric,  the 
extrema  of  the  first  derivative 
curve  would  be  equidistant  (in  cm  *) 
from  vQ  .  These  wavenumber  inter¬ 
vals  were  measured  for  a  number 
of  CO  spectra,  and  a  typical  set 
of  results  is  shown  in  Table  1. 

The  lines  measured  spanned  the 
region  26  to  81  cm  ‘.  In  the  first  column,  is  the  lower  state  rotational 
quantum  number,  (AvD)  is  the  wavenumber  interval  between  vq  and  the  ex¬ 
trema  of  the  first  derivative  curve  on  the  low  wavenumber  side  of  vQ ,  and 
(Avp)+  is  the  corresponding  interval  on  the  high  wavenumber  side  The 
last  column  gives  the  differences  (6)  between  (Av^)  and  (Av^)  +  The  agree¬ 
ment  between  these  two  quantities  is  quite  good;  the  average  value  of  6, 

-  3  -  1 

without  regard  to  sign, is  2.  3  x  10  cm  and  in  our  opinion  is  a  good 
indication  of  the  accuracy  in  measurement  attainable  with  our  instrument 
(at  least  for  isolated  lines  with  reasonably  good  S/N).  The  differences,  6  . 
are  of  random  sign,  thus  indicating  that  any  error  in  locating  x  =  0  in 
the  interferograrn  is  negligible  (see  Section  V). 

At  this  point  it  is  worthwhile  examining  interferograrn*  from  another 
point  of  view  (Ref  28)  In  a  recent  series  of  articles  Gordon  (Ref  29)  has 
developed  the  theory  concerning  the  information  obtained  when  the  band 
shapes  of  IR  and  Raman  bands  are  Fourier-transformed  Of  particular 
interest  is  his  treatment  for  pure  rotational  absorption  His  development 


Fig.  21.  The  J  =  16  —  17  transition 
of  CO.  The  circles  are  intensity 
points  computed  from  the  inter- 
ferogram.  The  squares  are  the 
first  derivative  points  computed 
frcm  the  intensity  points. 
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Table  1.  Fir»t  Derivative  Extremes  (Referenced  to  vQ)  for  the  Pare 
Rotational  Lines  of  Carbon  Monoxide4 


JL 

(AVq)  _ 

-1 

cm 

,4',dU  . 

-1 

cm 

6  =  (AvD)_  -  (Avd)+  > 
on.'1 

6 

0.  0477 

0.  0441 

+  0.  0036 

7 

0.  0498 

0. 0562 

-0. 0064 

8 

0.  0570 

0. 0539 

+  0.  0031 

9 

0.  0649 

0.  0618 

+0. 0031 

10 

0  0590 

0.  0593 

-0. 0003 

11 

0.  0649 

0.  0669 

-0. 0020 

12 

0.0593 

0.  0567 

+0  0026 

13 

0  0611 

0. 0625 

-0  0015 

14 

0  .  0615 

0. 0589 

+0  0026 

15 

0. 0599 

0.  0603 

-0. 0004 

16 

0  0531 

0.  0550 

-0  0019 

17 

0  0502 

0  0526 

-0  0024 

18 

0.  0539 

0.  0525 

+0  0014 

19 

0.  0496 

0.  0471 

40.  0025 

20 

0. 0516 

/ 

0. 0527 

-0. 001 1 

a 


From  an  interferogram  with  absorbing  path  -  15  cm,  p  =  608  Torr 


show,  that  a  a  pure  rotational  absorption  band,  suitably  normalized  to 
unit  area,  is  Fourier-transformed  (frequency  to  time)  then  the  dipole  mo¬ 
ment  correlation  function  <  u(o)  •  u(t)  >  is  obtained.  Gordon  interprets  this 
correlation  function  (in  the  classical  limit)  as  follows: 

"Imagine  that  one  could  observe  the  microscopic  re¬ 
orientation  of  a  single  molecule  in  a  system  of  many 
molecules  in  thermal  eouilibrium.  Suppose  that  at 
a  time  "o"  the  dipole  of  this  molecule  points  along  a 
direction  u(o).  Then  we  follow  the  thermal  motion? 
of  the  molecule  and  at  a  time  later  we  measure  the  pro¬ 
jection  of  u(t)  on  the  original  direction  :  u(o)  u(t). 

Now  we  make  this  measurement  again  and  agaTii,  picking 
out  different  reference  times  "o  '  Finally,  we  average 
all  of  these  trajectories  u\o)  ■  u(i)  to  obtain  the  correla¬ 
tion  function  <  ufo)  u(t)  >  .  " 

An  interferogram  is  the  FCT  of  the  spectrum  and  is  the  intensity 
passed  by  the  rometer  measured  as  a  function  of  *he  optical  path 

difference  generated  by  the  interferometer  The  optical  path  difference 
is  linearly  related  to  the  time  of  retardation  by  the  velocity  of  light;  thus, 
irttrrie vagravns  may  be  cosily  converted  from  the  optical  path  differ*  .ce 
scale  to  the  time  retardation  scale  by  dividing  by  the  velocity  of  light.  Early 
interferograms  (Refa.  23,  24)  are  complicated  by  the  tact  that  features  of 
the  inte rfe rogram  corresponding  to  the  absorption  of  the  sample  are  super¬ 
imposed  on  the  source  inte  rfe  rogram  By  using  a  "double -beam"  (reference 
and  sample  beam)  two-beam  interferometer  we  can  to  a  large  extent  sub¬ 
tract  the  source  inte  rfe  rog  ram .  and  suitably  normalise  the  resulting  inter- 
ferogram  so  that  it  corresponds  to  the  dipole  moment  correlation  turret  urn 
described  by  Gordon  ; the  correspondence  is  not  exact,  see  Ref  30)  A 
portion  of  such  inte  rfe  rog  ram  is  shown  »n  F’-g  22  tor  gaseous  CO 
(p  n78  dorr,  and  absorption  path  lengih  lb  cm) 

The  interpretation  of  the  inte  rfe  rog  ram  on  the  basis  of  wavelength 
and  optical  path  difference  has  been  given  previously  (Ret  23  3  1  )  and 

the  interpretation  on  the  basis  of  frequency  and  time  is  easily  «■  een  from 
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velocity  of  light. 


Gordon's  development.  In  the  semiclassical  limit  at  a  time  t  =  0,  a 

molecule  in  the  Jth  angular  momentum  state  (  =  -ft  [J(J+1)]  J/I, 

where  I  is  the  moment  of  inertia  of  the  molecule)  is  perfectly  correlated 

(with  the  system  and  itself)  since  t  =  0  is  chosen  as  the  reference  point. 

After  the  time  of  one  revolution  has  passed  the  correlation  will  now  be  less 

(or  even  negative),  since  molecules  in  other  different  angular  momentum 

states  (0,  1,  2 . J-l,  J+ 1 ,  J+2,  ...)  will  have  made  more  or  less 

revolutions  After  times  longer  than  the  time  of  average  rotation,  the 

correlation  function  will  oscillate  (with  decreasing  amplitude)  about  the 

zero  value.  However,  in  the  case  of  diatomic  and  linear  molecules  we  can 

expect  the  correlation  function  to  become  significantly  positive  again  for 

times  longer  than  the  average  rotational  period.  This  will  happen  when 

the  molecule  in  the  Jth  state  has  made  J  rotations,  then  the  molecules  in 

any  other  J1  state  will  have  made  J'  rotations,  and  if  the  system  has  not 

been  disturbed  the  state  of  the  system  at  this  time  will  duplicate  (ne- 

glecting  translation)  the  state  of  the  system  at  t  =  0.  Also,  on  the  basis  of 

this  interpretation,  smaller  but  quite  significant  negative  correlations  are 

expected  just  before  and  after  the  large  (-1)  positive  correlation.  This 

behavior  is  shown  quite  clearly  in  Fig.  22.  For  CO  a  large  positive- cor  - 

- 12 

relation  occurs  at  t  -  8.7  x  10  sec  ,  the  dipole  correlation  function 
attaining  the  value  -0,  90  at  this  time.  Extending  the  argument  to  longer 

times  shows  that  the  correlation  function  should  show  the  same  behavior 

/ 

at  integer  multiples  of  this  time.  Indeed,  the  interferogram  from  which 
Fig.  22  was  taken  showed  24  such  features  (out  to  t~2  x  10’  U  sec).  Two 
general  observations  mi;  he  mMe  abua*.  Viiwwt  k  *  (veiled  C 

in  interferometry). 

First,  the  correlation  peak  for  any  given  signature  is  always  less 
than  the  one  preceding  it.  This  is  understandable,  using  a  semiclassical 
approach.  Comparing  the  state  of  the  system  at  t  =  0  with  some  later  time 
t',  we  see  that  some  of  the  molecules  will  have  had  their  rotational  motion 
interrupted  (by  inelastic  collisions,  or  even  elastic  collisions,  that  d:  slurb 
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the  phase  of  the  rotation),  and  thus  these  molecules  will  be  uncorrelated 
(on  the  average).  At  moderately  low  pressures  and  above,  this  collisionai 
type  interruption  is  the  main  contribution  to  the  breadth  of  the  spectral 
line.  In  fact,  it  is  possible  to  analyze  the  fall-off  of  intensity  (or  corre¬ 
lation  factor)  of  these  signatures  for  line-width  parameters. 

Second,  the  signatures  are  not  symmetric  about  the  maximum  posi¬ 
tive  correlation.  The  negative  correlation  just  before  the  large  positive 
correlation  is  less  negative  than  the  one  just  after.  This  is  because  .150 
molecule  is  perfectly  rigid  but  rather  distorts  when  rotating  (centrifugal 
distortion);  therefore,  even  in  the  semiclassical  limit  the  angular  velocities 
of  the  various  states  are  not  integer  multiples  of  the  lowest  rotating  state 
but  rather  converge  slightly,  since  the  moment  of  inertia  increases  slightly 
for  each  higher  state  of  angular  momentum. 

Quantum  mechanically  the  above  arguments  would  have  to  be  modi¬ 
fied  since  it  is  not  the  rotational  frequencies  that  are  observed,  but  rather 
the  differences  in  rotational  frequencies  (usually  AJ  =  ±1).  It  is  expected 
that  the  more  rigorous  argument  will  lead  to  much  the  same  interpretation. 

2.  DEUTERIUM  CHLORIDE 

A  portion  of  the  pure  rotational  interferogram  of  deuterium  chloride 
(DC1)  is  shown  in  Fig.  23.  The  expected  signatures  characteristic  of  dia¬ 
tomic  and  linear  molecules  are  seen  out  to  x*2.  5  cm.  Between  2.  5  and 

3.  0  cm  there  are  still  variations  in  the  interferogram,  but  no  signatures 

are  apparent.  After  about  3.  2  cm  the  signatures  reappear,  although  they 

are  broader  and  down  considerably  ;n  intensity.  This  is  interpreted  as  due 

to  the  isotopic  splitting  of  the  pure  rotational  spectrum  of  DC1,  due  to  the 

35  37 

two  chlorine  isotopes  Cl  and  Cl  .  DCl  is  essentially  a  mixture  of  two 

35  37 

molecular  species  DCl  and  DCl  The  rotational  constants  for  the  two 
molecules  are  slightly  different;  thus,  the  interferogram  consists  of  two 
sets  of  signatures  that,  at  small  x  (low  resolution)  are  in  step  (superpose 
constructively),  at  larger  x  are  out  of  step  (superpose  destructively),  and  at 
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Fig.  23.  The  pure  rotational,  double-beam  interferogram  of  DC1, 


still  larger  x  get  back  in  step  one  lobe  removed.  The  fact  that  this  phenom¬ 
enon  is  observed  in  the  interferogram  indicates  that  at  least  some  of  the 

3^  37 

isotopic  lines  in  the  DC1  '  and  DCi  pure  rotational  spectrum  are  resolved. 

Two  other  comments  are  in  order.  The  last  part  of  the  interfero¬ 
gram  (x  *  13.  5  cm)  is  shown,  and  the  estimated  signal  (true  variation  about 
average  value  of  I(x))  to  oise  is  about  one.  Also,  there  is  a  drift  down¬ 
ward  (in  intensity)  of  the  interferogram.  This  is  instrumental  in  nature 
and  is  unimportant,  since  a  long  term  drift  (in  x)  corresponds  io  extremely 
low  wavenumber  (frequency)  components  in  the  spectrum  (0  to  ~2  cm"*). 

In  Fig.  24  are  shown  selected  regions  the  pure  rotational  spectrum 
of  DCI.  The  J  =  i-*2  (21.  5  cm'1),  J  =  2 -*3  (32.  3  cm'1),  and  J  =  3~*4 
(43.  C  cm’1)  regions  are  shown.  The  J  =  2— *3  transition  is  clearly  split, 
and  it  is  evident  that  the  resolution  displayed  .5  better  than  the  separation 
of  the  lines,  which  is  0  094  cm.  The  J  =  1  —►2  is  not  split  but  is  definitely 
asymmetric.  The  calculated  splitting  is  0.063  cm’1  (Refs.  32,  33).  The 
theoretical  resolution  of  this  run  was  0.  074  cm’1.  It  is  clear  from  the 
mora  th«m  adequately  split  J  =  2 -►3  transition  and  the  unresolved  but 
definitely  asymmetric  J  =  1—+2  transition  that  the  theoretical  resolution 
is  closely  approached. 

3.  NITRIC  OXIDE 


Nitric  oxide  (NO)  has  been  of  great  scientific  interest  because  of  its 
unique  property  of  being  a  stable  diatomic  molecule  having  an  odd  number 
of  electrons.  (The  pure  rotational  spectrum  of  NO  is  reported  .'n  Ref  34. ) 

The  ground  state  of  NO  ®  st-'te  v'h’rh  i«  split  into^n^ and  ^3/2 
components  by  spin-orbit  interaction.  The  splitting  is  small  ( =  122  cm"1), 
resulting  in  appreciable  population  of  both  components  at  room  temperature. 
The  consequence  of  this  splitting  in  the  pure  rotational  spectrum  of  NO  is 
the  appearance  of  two  series  of  lines  that  can  be  assigned  to  the  two  compo¬ 
nents  of  the  ir  state.  Rotational  constants  can  be  de  r iv ed TorT> o t h  sene*  , 
but  because  of  the  perturbation  by  the  iow-lying  electronic  state  the  constant? 
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ao  derived  are  only  effective  constant*  which  are  combinations  of  the  usual 
rotational  constants  and  additional  terms  due  to  the  spin-orbit  perturbation. 

The  microwave  spectrum  cf  NO  was  reported  by  Gallagher  and 
Johnson  {Ref.  35),  who  determined  the  effective  molecular  constants  for 
the  8tate  and  calculated  the  constants  for  the  state.  Favero, 

Mirri,  and  Gordy  (Ref.  36)  more  recently  have  reported  the  millimeter  wave 
spectrum  of  the  state  and  derived  constants  for  it.  Recently  there  has 

been  a  great  deal  of  work  in  the  near  -  and  intermediate -IR  spectral  regions. 
Areas,  et  al.  (Ref.  37)  studied  the  3-0  band  of  NO  with  a  precision  of  ~0.  03  cm"!* 
Olman,  et  al.  (Ref.  38)  studied  this  band  and  also  the  2-0  band  with  a  pre¬ 
cision  of  about  0.  0 1  cm-1.  Meyer,  at  al.  (Ref.  39)  studied  the  2-0,  3-0,  and 
4-0  bands  with  a  precision  of  better  than  0.  01  cm’*.  James  arxl  Thibault 
(Ref.  40)  studied  the  fundamental  region  with  a  precision  comparable  to  that 
of  Olman,  et  al.  Palik  and  Rao  (Ref.  41)  reported  results  obtained  in  the 
FIR  spectral  region  with  a  precision  of  about  0.  05  cm'*. 

7het«  i»  •igitiiicaui  ulvagi In  values  of  the  effective  constants 
as  determined  by  the  numerous  experimenters,  and,  in  particular,  the 
effective  rotational  constants  determined  in  previous  IR  investigations 
(Refs.  37-41)  differ  significantly  from  those  determined  in  microwave  in¬ 
vestigations  (Refs.  35,  36). 

Therefore,  in  order  to  resolve  the  discrepancies  betv/een  microwave 
and  IR  results,  it  was  decided  to  re -investigate  the  pure  rotational  spectrum 
ot  NO  at  high  resolution.  Ti.-  ’*f  -  -ilim  -  were  determined  once  with 
high  resolution  (  *  0.  06  cm**)  and  were  corroborated  by  a  lower  resolution 
(  *  C.  1  cm  )  run.  The  precision  of  the  line  positions  is  estimated  to  be 
about  0.  003  cm  except  at  the  high  and  low  extremes  of  the  spectrum 
where  the  S/N  was  poorer. 

The  parameters  for  the  high  resolution  run  are  as  follows:  The 

-mteHerogyfrm  wa#  digitally  sampl«d~ai— ' 40.-. ^optical  path,  dilferejic.? _ 

increments.  Each  point  was  sampled  for  7  seconds  to  improve  the  S/N.  The 
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maximum  optical  path  difference  generated  was  160  mm  A  total  of  4000 
points  were  recorded  in  the  interferogram.  The  transf  jrmed  spectmm  was 
computed  at  0.  02  cm  *  interval* 

Air  Products  and  Chemicals  NO  of  99.  0%  purity  was  used.  The 
gas  was  additionally  purified  by  passing  through  a  dry  ice-acetone  trap. 

The  sample  was  slightly  contaminated  by  NO^  as  was  evidenced  by  its  pale 
yellow  color  when  viewed  through  the  cell.  Some  of  the  weak  unidentified 
features  in  the  spectrum  may  be  due  to  NO^  absorption.  The  sample  was 
contained  in  a  15-cm  path  length  cell  with  crystalline  quartz  windows  at 
a  pressure  of  399  Torr  at  room  temperature  (298°K). 

Figure  25  shows  selected  regions  of  the  interferogram.  Many 
qualitative  features  of  the  pure  rotational  spectrum  of  NO  can  be  deduced 
directly  from  the  interferogram.  First,  it  can  be  seen  that  there  are 
actually  two  sets  of  signatures  (denoted  by  the  "a"  and  "b"  sequence)  that 
are  well  separated  beyond  about  30-mm  optical  path  difference.  It  follows 
from  this  that  the  spectrum  will  consist  of  two  sets  of  almost  uniformly 
spaced  lines,  the  spacing  for  the  "a"  sequence  being  larger  than  that  for  the 
"b"  sequence.  Second,  the  "a"  sequence  is  distorting  (smearing  out)  more 
rapidly  than  the  "b".  It  follows  then  that  lines  corresponding  to  the  "a"  set 
will  show  a  greater  convergence  (because  of  centrifugal  distortion  and  higher 
order  terms)  than  those  of  the  ’b"  set.  Finally,  adjacent  signatures  in  each 
set  are  nearly  inverted  copies  o(  each  other,  which  indicates  that  the  quantum 
number  7  -V  r<?rgy  ’*'*'?!  exprcs  !on  is  nonintegr*!  Intact,  since  adjacent 
signatures  are  completely  inverted  with  respect  to  each  othei ,  J  is  half-integral. 

Figure  26  shews  the  spectrum  computed  from  the  interferogram 
shown  in  Fig.  25.  The  relative  intensities  of  the  lines  are  influenced  by 
the  envelope  of  the  source  output  spectrum  and  are  not  true  indications 
of  the  band  contour.  The  source  output  spectrum  has  been  discussed  pre¬ 
viously.  This  run  was  not  calibrated  for  intensity,  but  the  shape  of  the 
source  spectrum  is  very  similar  to  that  shown  in  Fig.  20. 
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Table  2  lists  the  frequencies  of  the  lines  observed  and  their  assign¬ 
ments,  together  with  deviations  from  frequencies  calculated  using  constants 
derived  from  analyses  of  these  data.  Frequencies  of  peak  absorption  were 
computed  by  the  method  of  zero  slope  from  the  0.  02  cm”  *  interval  computer 
output.  As  mentioned  above,  the  accuracy  of  the  observed  frequencies  is 
+  0.003  cm"*.  The  frequencies  are  quoted  to  0.0001  cm"*.  Little  confidence 
is  placed  in  the  last  figure,  but  it  is  included  to  offset  roundoff  errors  in  the 
comparison  with  calculated  values.  Transitions  not  used  in  the  least  squares 
analyses  are  enclosed  in  parentheses. 


Table  2.  Observed  Frequencies  and  Deviations  of  Calculated 
Frequencies  of  Nitric  Oxide 


J  +  1/2 

Obsarvad, 

cm'1 

2 

'1/2  , 
(obs-cilc) 

*10^,  cm  ^ 

(obs-calc)^ 
XIO4,  cm-1 

Obterv«d, 

cm*1 

2 

'3/2 

(obs- calc) 
X104,  cm’1 

(oba-calc)1* 
xlO4,  cm"1 

5 

(18.392,) 

+25 

+24 

(18  907,) 

-73 

-78 

6 

21.731, 

-19 

-20 

22.  350, 

-12 

-18 

7 

25.  075, 

-16 

-17 

25.784, 

-  9 

-16 

8 

28.  422 ^ 

♦21 

+21 

29.220, 

+22 

+  15 

9 

31.7640 

+  8 

+  7 

32.  648, 

0 

-  6 

10 

35.  107, 

+  8 

+  8 

36.078. 

+  16 

+n 

11 

38. 450, 

+14 

+14 

39. 502, 

-  3 

.  i 

12 

41.793, 

+21 

+21 

42. 926, 

+  2 

0 

13 

45.  132, 

-  7 

-  7 

46  346, 

-  4 

-  5 

14 

48.  475, 

♦  11 

4i: 

49.763, 

-  2 

-  2 

IS 

51.814, 

-  7 

-  7 

53.  178, 

+  5 

+  7 

16 

55. 155z 

+  6 

+  6 

56  5882 

-  9 

+  6 

17 

58.  492, 

-12 

-12 

59  996, 

-  6 

-  i 

18 

61.829, 

-17 

-16 

63.  399, 

-10 

-  4 

19 

65.  1*7^ 

-  5 

-  5 

66.800, 

+  8 

+  15 

20 

68. 501, 

-15 

•  14 

70. 198, 

+25 

+31 

21 

71.  837  j 

+  3 

+  4 

73.  586, 

-11 

-  5 

22 

75. 166, 

-17 

-17 

76.  973, 

-25 

-21 

23 

78.  499, 

+  7 

+  7 

80.  357, 

-25 

-25 

24 

81. 828, 

+27 

+26 

83.740, 

+  14 

+  9 

2S 

85.  150, 

+  1 

0 

87.  1 14, 

+  2 

-10 

26 

<88.  478,) 

+71 

♦69 

(90.  482,) 

-26 

-47 

27 

(91784,) 

-48 

-51 

(93.856,) 

+  47 

+  12 

‘Ui,l  equarss  analysis  dstermining  Bq,  Dq,  and  Hq 

^Lsast  squarss  analysis  holding  B  at  microwave  value  (L. R*f.  J5; 
2  o  mc 

Raf.  36)  and  dstsrmining  Dq  and  H^. 
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The  data  were  analyzed  by  computing  a  least  squares  fit  to  the  fol¬ 
lowing  power  series  in  (J  +  1): 

v  =  2B JJ  +  1)  -  (4D  -  2H  }{J  +  l)3  +  6H  (J  +  l)5  (12) 

o  o  o  o 

where  vis  the  frequency  of  the  transition,  J  the  quantum  number  of  the  lower 
state,  Bq  the  effective  rotational  constant,  Do  the  effective  centrifugal  distor¬ 
tion  constant,  and  Hq  the  effective  higher  order  distortion  constant,  which  is 
not  negligible  for  NO,  Least  squares  analyses  were  made  on  each  series 
and  on  the  difference  frequencies  between  the  two  states.  It  ia  felt  that  the 
difference  frequencies  are  more  accurate  than  the  frequencies  themselves  , 
since  the  differences  tend  to  minimize  any  errors  in  the  wavenumber  scale 
such  as  might  be  generated  by  a  phase  erro-  :  ie  interferogram. 

The  first  line  in  Table  3  show  a  the  result  of  these  analyses.  The 
errors  given  are  95%  confidence  inervale.  The  diifcrences  between  the  ob¬ 
served  frequencies  and  those  caicilateJ  uainj  the  constants  determined  by 


Table  3.  Effective  Rotational  Con  *  tents  of  Nitric  Oxide  (in  cm'  ) 


Reference 

VK,*r!/Z) 

Vff<S/2> 

ABceff(2+3/2-  2»1/2) 

D„eff(2.i/2) 

jj 

:  i 

-- 

This  work 

1.671854*0.  000081 

1.720178*0.000062 

0.  048323*0.  000057 

(3.4*3.  0)xl0'7 

i 

-- 

Thie  work* 

[1.6718614] 

[1.7196531] 

[0.04779172] 

(3.6+1  4)xl0*7 

f 

-- 

This  work** 

[1.6718614] 

[1.7202435] 

[0.04838213] 

(3.6+1.  4)xl0'7 

35 

Microwave 

1.6718614 

1. 7196531 

1. 13xl0'6 

36 

Microweve 

1.7202435 

V 

40 

Nea.'-IR 

1. 67229*0.  00025 

1.72014+  0.00025 

(7.728+  3.  51  10"7 

37 

N<--r  IR 

1.67199 

1.71958 

1.77xl0*6 

t 

f 

38 

Ne»r-'.R 

1.67233*0.00019 

1.72012+0.  00018 

(1.6*0.  3)xl0’6 

39 

Near-IR 

1.67223 

1.72002 

1  5xI0‘6 

4  1 

F1T3 

1. 6’20 

1.7198 

1.8x10  6 

!  > 
i 

r 


these  analyses  are  shown  in  the  third  and  sixth  columns  of  Table  2.  The 

average  deviations  from  frequencies  used  in  the  analyses  for  either  state 

ar<»  about  £0.  0012  cm  *.  Since  the  analyses  reproduced  the  microwave 

rotational  constants  so  well,  it  was  decided  to  use  the  microwave  values  of 

the  effective  BQ's  and  perform  least  squares  analyses  for  the  effective  Do's 

and  Ho's,  The  results  of  these  analyses  are  shown  as  the  second  and  third 

lines  in  Table  3.  The  second  line  is  the  result  of  using  Gallagher  and 

2  2 

Johnson's  {Ref.  35)  values  for  both  Bq(  rr^^)  an<*  Bq(  tt^^).  The  assumed 

values  are  given  in  brackets.  {The  constants  of  the  8tate  derived 

by  the  analysis  reproduced  the  observed  frequencies  with  significantly 

greater  deviations  of  a  nonrandom  character  than  the  constants  given  in 

2 

the  first  line,  indicating  that  the  value  of  Bq(  calculated  by  Gallagher 

and  Johnson  is  in  error.  The  third  line  is  the  result  of  using  Gallagher 

2 

and  Johnson's  (Ref.  35)  value  of  Bq(  ^ )  and  Favero,  Mirri,  and  Gordy's 

(Ref.  36)  value  of  B^it^^).  The  constants  derived  from  these  \nalyses 


Table  3,  Continued 


Doeff(2r3/2)  ADoeff(2w3/2-2*1/2) 

Vo*“'\lZ>  UHoem\n-\/2) 

<10.24±0.  23)xl0"6 
(8.  51±0.40)xl0'6 
(10.  46t0.  I2)xl0*6 
10.  64xl0’fc 

(9.  90+0.  21)xl0'6 
(8.  14+0.  40)xl0-6 
(10. 10±0.i0)xl0‘6 

(- 1. 42±.C.  29/xlO*” 
(- 1 . 40±0.  18)xI0'9 
(-1.  40+0.  »8)xU>'9 

(+  1 . 00+0.  23)xl0*9 
(-0.  51  +  0.  53)xl0'9 
(+1.  19±0.  151xl0'9 

(+2.42±0.  21)xl0'9 
(+0.  89±0.  54)xl0'9 
(+2. 59±0.  14)xl0’9 

(10.  70*0.  19)xl0'6 

9.  77x10 

(9.  5+0.  3)xl0’6 

9. 3xl0'6 

10.  OxtO*6 

-2. 091xi0'9 (calc) 

+  2.  093xl0'9(ca>c) 

^Constants  derived  by  ••turning  the  microwave  valuei  o f  B  from  Ref.  35 
and  performing  least  inmm  cilci’lat;on»  *o  determine  Dq  and  Hq. 

^Constants  derived  by  assuming  the  miciowav;  value  of  *3  rff(2ir,  ,,)  fr' ~n 
Ref.  35  and  B  eff(2ff  .  )  from  Ref.  36  and  performing  leSst  sqdarer  cal¬ 
culations  to  dAermine  TD  and  H  . 

o  o 
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have  significantly  smaller  confidence  intervals  and  reproduce  the  observed 
frequencies  as  well  as  the  complete  analysis.  The  differences  between  the 
observed  frequencies  and  those  calculated  using  the  constants  given  in  the 
third  line  arc  given  in  the  fourth  and  seventh  columns  of  Table  2.  The  aver¬ 
age  deviations  for  either  state  are  about  0.  0012  cm"  ,  the  same  as  for  the 
complete  analysis. 

The  values  of  the  difference  constants  are  also  included  in  Table  3. 

It  should  be  noted  that  the  confidence  limits  are  smaller  than  for  either 
state  separately.  This  corroborates  the  belief  that  the  difference  frequencies 
are  known  better  than  the  frequencies  themselves.  Finally,  xable  3  includes 
a  tabulation  of  constants  derived  by  other  experimenters  given  the  spectral 
region  in  which  data  were  obtained. 

Once  the  values  of  the  effective  constants  have  been  determined,  the 
unperturbed  ground  state  rotational  constants  can  be  obtained  from  them. 
Using  the  formulation  developed  by  James  and  Thibault  (Ref.  40)  the  true 
constants  are  given  by 

Bo  =  1/2  (Bol  +  Bo2  -  3Do> 

Table  4.  Rotational  Constants  of  ^o  ^ol  +  ^o2^ 

Nitric  Oxide 

Ho  =  1/2  (Hol  +  Ho2) 

where  the  subscripts  1  and  2 
2  2 

denote  the  ir^2  and  '"'3/2 
states,  respectively.  The  con¬ 
stants  obtained  from  these  re¬ 
lations  are  given  in  Table  4. 

The  errors  quoted  are  approxi¬ 
mately  95%  confidence  limits. 


B  =  1.696008  1  0.  000072  cm 
o 

Dq  =  (5.29  1  0.27)  x  10"6  cm"1 

H  =  (-2.  1  12.6)  x  10"10  cm"1 
o 
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One  constant  of  NO  which  is  of  great  interest  is  the  spin-orbit 
coupling  constant,  Aq.  The  conventional  relation  (Ref.  40)  used  for  calcu¬ 
lating  Aq  is 


Using  the  values  of  these  constants  given  in  Tables  2  and  3,  one  obtains  a 

„  1 

value  of  Aq  of  122.  49  cm  *.  But,  as  James  and  Thibault  (Ref.  40)  have 
pointed  out,  this  formulation  is  incorrect.  They  derived  a  rather  extensive 
series  of  equations  from  which  the  correct  value  of  Aq  can  be  obtained.  Un¬ 
fortunately  these  equations  depend  strongly  on  the  values  of  Dqj,  D^,  and 
D  ,  constants  which  have  not  been  well  determined.  Using  these  equations 
for  Aq>  a  value  of  128.  6  cm  *  is  obtained.  Little  confidence  can  be  placed 
in  this  value,  however,  for  reasons  stated  above.  Information  obtained  in 

the  vibrational  spectrum  of  NO,  namely  the  difference  in  band  origins  of 

2  2 

the  satellite  band  and  the  band,  gives  a  direct  measure  of  6q. 

The  value  of  Aq  is  very  sensitive  to  this  parameter.  Aq  is  related  to  6q  by 
the  following  equation  (Ref.  39): 

f  2  3  2l 

A  =  64  +  3(3  +y  D  )4  +  2B 

o  1.0  o  2  o  J  o 

Using  James  and  Thibault's  (Ref.  40)  value  of  6q,  119.732  cm'1,  and  the 

values  of  B  and  D  from  Table  3,  one  obtains  a  value  of  A  of  123.  1600  cm-1, 
o  o  o 

The  accuracy  of  this  value  is  limited  by  the  precision  with  which  6Q  is  known. 
The  inaccuracies  of  the  other  constants  are  in  the  fifth  and  greater  decimal 
places. 

Included  in  Table  3  are  values  of  the  effective  rotational  constants 
obtained  in  the  other  investigations.  It  should  be  noted  that  our  value  of 
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2 

B^(  *1/2)  rePr°duceB  the  microwave  value  of  Gallagher  and  Johnson  (Ref.  35) 

well  within  our  experimental  uncertainty.  Our  uncertainty  does  not  bracket 

any  of  the  other  IR  values  of  this  constant.  It  should  be  pointed  out  that  in 

two  of  the  near-IR  investigations  (Refs.  38,  40)  the  data  were  not  analyzed 

7 

using  the  well-known  combination  difference  relations  ,  but  rather  were 
directly  analyzed  in  terms  of  a  power  series  in  J  whose  coefficients  are 
combinations  of  upper  and  lower  state  constants.  This  type  of  analysis 
suffers  the  usual  drawbacks  of  data  fitted  to  a  high  order  polynomial,  many 
coefficients  of  which  are  small.  The  other  near-IR  investigations  (Refs,  37, 
39)  did  use  combination  relations. 

Our  value  of  Bq(^tt is  bracketed  by  the  two  microwave  values 
and  does  not  overlap  either  within  our  uncertainty.  No  uncertainty  was 
quoted  with  either  microwave  value.  We  are  very  close  to  the  value  of 
Favero,  et  al.  (Ref.  35).  However,  when  the  microwave  values  of  the 
effective  Bq's  were  used  in  the  least  squares  analyses  of  our  data,  only  the 

effective  B  of  Favero,  et  al.  (Ref.  35)  gave  good  agreement  with  observed 

0  2 
frequencies  and,  importantly,  the  agreement  using  their  value  of  Bq( 

was  no  different  than  that  using  our  value. 

The  values  of  the  centrifugal  distortion  constants  are  much  less 

accurately  determined.  In  general,  the  value  of  Dq(  ir^^)  agree*  with  other 

values,  «  10  X  10  ^  cm  The  uncertainty  of  this  value  is  of  the  same  order 

2 

as  that  given  for  other  values.  However,  the  value  of  Dq{  it^)  is  con¬ 
siderably  lower  than  any  other.  This  value  is  the  least  well  determined  of 
all  the  effective  constants.  One  overriding  factor  with  respect  to  both  Dq 
and  Hq  is  that  our  data  are  the  only  data  accurate  enough  to  allow  significant 

7 - 

Because  many  of  the  lines  in  the  vibration -rotation  bands  were  not  resolved 
or  were  blended  with  overlapping  water  vapor  lines,  the  authors  of  Refs.  38 
and  40  felt  there  were  insufficient  combination  relations  involving 
unblended  lines  to  permit  an  accurate  analysis. 
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determinations  o£  the  effective  higher  order  distortion  constants,  Hq.  Our 
values  of  the  effective  H^'s  agree  in  sign  with  the  calculated  values  of  James 
and  Thibault  (Ref.  40).  Their  data  were  not  accurate  enough  to  permit  a 
significant  experimental  determination  of  these  constants. 

4.  WATER  VAPOR 

K.  M.  Randall's  1938  review  article  (Ref.  42)  on  FIR  spectroscopy 
aptly  described  the  FIR  absorption  spectrum  of  water  vapor  as  "a  complex 
and  apparently  chaotic  arrangement  of  lines.  "  Since  the  first  observation 
by  Wright  and  Randall  (Ref.  43)  in  1933  of  a  water  vapor  spectrum  in  which 
rotational  structure  beyond  100p  was  resolved,  the  pure  rotational  spectrum 
of  water  vapor  has  served  as  a  test  of  instrumental  performance  and  a 
background  with  which  nonvacuum  spectrometers  have  had  to  contend. 

Early  workers  (Refs.  43-45)  in  the  FIR  spectral  region  were  able  to 
achieve  a  resolution  of  up  to  0.  5  cm*1  in  their  observations  of  the  spectrum 
of  water  vapor.  However,  using  those  results  the  energy  levels,  and  hence 
rotational  constants,  could  not  be  determined  as  accurately  as  possible 
from  the  near-IR.  vibration -rotation  spectrum.  Therefore,  observations  of 
the  pure  rotational  spectrum  of  water  vapor  were  used  only  for  illustrations 
of  instrumental  performance  and  secondary  frequency  calibrations  for  the 
FIR.  No  further  analysis  of  the  spectra  was  attempted.  Genzel  and  Eckhardt 
(Ref.  46)  first  reported  the  very  far  (beyond  200p.)  IR  spectrum  of  water 
vapor  in  1954.  Since  then  there  have  been  many  papers  (Refs.  47-55)  on  the 
FIR  spectrum  of  water  vapor;  however,  none  with  resolution  comparable  to 
that  obtainable  in  the  near-IR  vibration- rotation  spectrum.  Four  separate 
measurements  of  the  spectrum  were  made  using  a  cell  with  15-cm  absorbing 
path  length  and  with  pressures  of  7.  5,  10.  0  (2  runs),  and  16.  5  Torr  at  room 
temperature  (298°K).  The  modular  resolution  of  the  four  runs  was  between 
0.  06  and  0.  09  cm"  1 . 

The  parameters  for  the  runs  are  as  follows.  The  interferograms 
were  digitally  sampled  at  40-p  optical  path  difference  increments.  Each 
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point  was  sampled  for  5  to  7  seconds  (constant  for  a  given  run)  to  improve  the 
S/N.  The  maximum  optical  path  difference  generated  was  between  114  and 
160  mm.  Between  2850  and  4001  points  were  recorded  in  each  interferogram. 
The  transformed  spectrum  was  computed  at  approximately  0.02  cm’*  intervals. 

Figure  27  shows  the  spectrum  from  5  to  125  cm  *  computed  from  one 
of  the  interferograms.  This  spectrum  was  obtained  with  10.  0  Torr  pressure 
(298°K)  in  a  15-cm  cell  and  has  a  modular  resolution  of  0.063  cm  *.  The 
relative  intensities  of  the  lines  are  influenced  by  the  envelope  of  the  source 
output  spectrum  and  are  not  true  indications  of  actual  intensities.  These 
runs  were  not  calibrated  for  intensity,  but  the  shape  of  the  source  spectrum 
is  very  similar  to  that  shown  in  Fig.  20. 

Below  the  spectrum  are  shown  the  calculated  (Ref.  56)  frequencies  and 
relative  intensities  for  pure  rotational  transitions  of  water  vapor  using  the 
energy  level  scheme  of  Gates,  et  al.  (Ref.  57).  These  relative  intensities 
have  not  been  normalized  to  the  source  spectrum. 


Figure  28  shows  details  of  two  portions  of  the  spectrum  illustrating 
the  resolution  attained.  Figure  28a  shows  the  triplet  at  38  cm  The 
calculated  separation  of  the  541  —  6^  (38.  642 ^  cm"1)  and  321  791^ 

These  results 


cm  )  transitions  is  0.  149  cm’*.  Figure  28b  shows  the  doublet  at  59  cm 
The  calculated  separation  of  two  components  is  0.0678  cm 


are  the  first  published  to  resolve  this  doublet. 


*1, 


Table  5  lists  the  frequencies  of  the  lines  (in  cm  )  observed  and 
their  assignment,  together  with  calculated  frequencies,  calculated  relative 
intensities  (Ref.  56),  and  observed  relative  strength.  The  calculated  fre¬ 
quencies  were  computed  from  energy  levels  determined  in  this  study.  The 
calculated  relative  intensities  are  normalized  to  100  for  the  ^q3~*^12 
tions  at  36.6055  cm"1,  the  line  of  greatest  apparent  intensity  in  our  spectra. 
The  relative  strengths  are  rough  guides  of  apparent  intensity,  and  therefore 
precision,  of  the  lines  as  observed  in  our  spectra.  Whenever  possible, 
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Table  5.  Pure  Rotational  Transitions  of  Water  Vapor 


Observed 

Frequency, 

-1 

cm 

Calculated 

Frequency, 

-1 

cm 

Relative 

Strength 

Calculated 

Relative 

Intensity 

Assignm< 

18.5706  (47)* 

(18.5706)** 

M 

13.7 

101“*  1 10 

25.085.  (19) 

->  - 

(25.0855) 

S 

15.  1 

202”~211 

30.  556g  (45r 

30. 560^ 

w 

0.8 

331~422 

32.3655  (42) 

32.364- 

9 

w 

1.2 

43l"~524 

32.9512  (10) 

(32.  9 5 1 2) 

s 

6.6 

1 1 1~*  202 

36.6055  (15) 

(36.6055) 

s 

100.0 

3  —  3 

J03  12 

37. 1 37 j  (7) 

(37. 137j) 

s 

4.4 

i 

38.465Q(3) 

38.4682 

s 

10.9 

22I~*312 

38.  642^  (79) 

38.  643 

w 

2.2 

541~634 

38.7916  (18) 

(38. 7916) 

s 

108.6 

3 12~*  321 

40.2830  (10) 

(40. 2830) 

s 

43.6 

413~*422 

40. 987 j  (6) 

(40. 987 3) 

s 

20.9 

211~220 

42.633.  (34)2 

42.622 

w 

0.7 

652~743 

43.250q  (34r 

43.248 

w 

0.8 

734~*827 

44. 098^  (20) 

(44.099) 

M 

5.2 

532^625 

47.0529  (29) 

(47.053) 

s 

135.7 

514~523 

48.0602  (58) 

(48.0602) 

w 

1.0 

633^726 

51.4376  (57) 

51. 433g 

w 

1.4 

542~633 

53. 444 6  (20) 

(53.4446) 

s 

54.3 

404"*413 

55.4063  (16) 

(55.406j) 

s 

79.1 

212~221 

55.7046  (78) 

(55.7046) 

s 

113.0 

101**212 

57.269q  (31) 

(57.2690) 

s 

175.7 

212~*303 

58.772,  (32) 

(58.7722) 

8 

34.5 

624-633 

58.9030  (79 

58.  W0 

w 

1.5 

642~735 

68- 


Table  5.  (cort) 


Observed 
F  requency, 

_  i 

cm 


Calculated 
Frequency,  Relati 

ve 

cm  Strength 


Calculated 

Relative 

Intensity 


Assignment 


59.871^  (48) 
59.9393  (44) 
62,  3008  (19) 
62.8713  (-)1 
64.0104  (31) 

67,  240 8  (26) 

68.  059g  H2) 
69. I91g  (24) 

72.  186}  (16) 

73.  259^  (32) 
74. 1080  (30) 

74.  875j  (73) 

75.  5232  (24) 
77.3138  (38) 
78.  1936  (30) 
78.914?  (26) 
79.7754  (19) 
80.  997 ^  (37) 
81.6175  (43) 
82.  152g  (28) 
85.6295  (37) 
85.  7834  (-)1 
87,7594  (32) 


(59.  871 ?) 
<59.  939)' 
(62.  30  ^ ) 
62. 84, 


?64.  028  . 

64.  040"J 

(67,  ?.40) 

(68.059  \ 
8' 

(69.  1 9 1  g ) 
(72.  186  j ) 
(73.259^ 
(74. 10g) 
(74.  87^) 


(75.  523L,) 
[77.  32]*** 
(78.  193) 
(78.914?) 
(79.7774) 
[81.  GOi 

(81.  6lg) 


(82. 152g) 
<85.  629) 
{85,  80) 
(87.7594) 


s 

41.8 

63.  3 

6l5-624 
?25  ?34 

s 

149.  2 

523  ~532 

w 

s 

1.3 

55.  3 

22.  2 

752  ~845 

313-322 

4  —  5 

*32  523 

M 

12.  0 

8_  /  -  8,c 
2o  35 

S 

57.  1 

4  —4 

22  31 

s 

30.  1 

322~413 

s 

116.  5 

202~313 

s 

127.  4 

3  —  3 

3Z1  J30 

s 

206.  0 

505  ~514 

M 

4.  1 

743  "*836 

3 

233.  4 

414  ^423 

M 

12.  8 

936~945 

S 

101. 6 

7l6-725 

S 

47.  4 

322  ^331 

S 

180.  8 

313  404 

M 

18.  0 

927  -936 

M 

11.0 

835  844 

S 

222.  1 

423  -432 

M 

11.  1 

643  -?34 

W 

1.8 

U38  U47 

S 

78.  1 

524-533 

>69- 


Table  5.  (cont) 


Observed 

Calculated 

Calculated 

F  requency, 

-1 

cm 

Frequency 

-1 

cm 

Relative  Relative 

Strength  Intensity 

Assignment 

88. 077 

4 

(26l  f 

88.  079  j 

S 

706.  8 

3  —  4 

03  14 

88. 877  j 

(51) 

(88, 87g) 

M 

70.  0 

7 34  '  43 

89. 58 1 g 

(51) 

(89.  581g) 

M 

86.  2 

315~~524 

92.  5 1 5 5 

(245)  t 

92.  540^ 

M 

262.  0 

1 10  ~~  221 

96. 067, 

(99) 

(96.067) 

M 

198.  7 

625  -634 

96.2i28 

(32) 

96.  186g 

M 

69.  7 

606~'6i5 

96.  24fe 

40.  6 

633  ~642 

98. 8006 

(40) 

(98.  8006) 

21. 8 

533  “*  624 

99. 066g 

(179) 

99.0240 

M 

'  83.9 

ln  ~~  Z20 

|99.  14 

275.  8 

23  3 1 4 

10G.0210 

(54) 

[99.96] 

W 

1.4 

U47-n56 

[99.98] 

23.  7 

817  ^S26 

100. 52lg 

(256) 

[100.  2^ 

M 

3.  1 

r- 

O 

T 

CO 

<M 

O 

* 

100. 55r 

b 

975.  5 

414  *  ^05 

101. 5234 

(84) 

(101. 823) 

M 

160.  5 

532  ~^41 

104. 2884 

(61) 

(104.  2884) 

W 

44.  4 

4  —4 

31  40 

104. 5705 

(39) 

(-)1 

(104.  570s) 

M 

359.  0 

4  —  5 

04  D15 

105.  1 1 3  j 

(105.  12) 

W 

2.  2 

854  -945 

105. 637, 

(79)  j 

105,  59e 

M 

136.  5 

432  ~441 

D 

105.  658 

237. 2 

6l6 

106.  139q 

(65) 

(106.  1390) 

W 

58.  0 

533  -542 

107.075, 

o 

(96' 

(107.  075^) 

W 

47.  7 

?26  “~735 

107. 7382 

(54) 

(107. 73g) 

w 

148.  2 

634  -643 

Ill. 120, 

(139) 

Ill. 10j 

M 

33.  5 

?35  “*  '44 

111. 128g 

172.  9 

2ll-322 

Observed 

Frequency, 

-1 

ern 

Calculated 

Frequency, 

-1 

cm 

Relative 

Strength 

Calculated 

Relative 

Intensity 

Assignment 

116.5959 

(70) 

(Uv.  i96) 

W 

57.  7 

836  “~845 

117.06?8 

(154) 

[117.07] 

w 

13.  1 

945  ~"954 

117. 966^ 

(51) 

(117.  96?) 

w 

173.  6 

?07  ""?16 

120, G83g 

(134) 

(120.0838) 

w 

418.9 

515  ~60o 

120. 5049 

(142) 

O 

(120.  50 g) 

w 

89.2 

827  "*836 

121. 913g 

d<r 

121. 903 

w 

1303. 3 

505  ~616 

I 

r  (  )  average  deviation  for  four  observations  in  0,  0001  cm" 

%u%*> 


t 

1 

? 


(  )  transition  used  in  the  determination  of  the  energy  levels 
[]  calculated  from  the  energy  levels  of  Ref,  5? 

line  width  very  pressure  sensitive,  frequency  quoted  is 
from  the  lowest  p’  assure  observation 

only  observed  in  one  spectrum 
only  observed  in  two  spectra 
only  observed  in  three  spectra 
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38.20  0.40  0.60  0.80  39.00  S9.60  0.80  60.00  0.20 

v  icjir1} 


Fig,  <.8.  The  38.  fe  and  59.9  cm”  region  of  the  pure  rotational  spectrum 
of  water  vapor.  The  circles  are  the  intensity  points  computed 
from  the  inter  ferogram. 

strong  and  medium  lines  were  used  in  the  analysis.  The  frequencies  of  peak 
absorption  were  computed  by  the  method  of  zero  slope  from  the  computer 
output.  The  accuracy  of  the  observed  frequencies  for  strong  and  medium 
lines  is  estimated  to  be  ±0.  003  cm"  and  slightly  more  for  the  weak  lines. 
Frequencies  are  quoted  to  0.  0001  cm"  and  slightly  more  for  the  weak  lines. 
Frequencies  are  quoted  to  0.  0001  cm  *  only  to  offset  roundoff  errors  in  the 
analysis.  Little  confidence  is  placed  in  the  last  figure. 

The  transitions  listed  in  Table  5  were  assigned  from  calculations  based 
on  the  energy  levels  of  Gates,  et  al.  (Ref.  57),  which  were  the  most  recent 
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complete  set  of  energy  levels  published.  Once  the  transitions  were  assigned, 
new  values  of  the  energy  levels  could  be  determined  by  adding  up  transitions 
through  the  levels  from  the  ground  state. 

Since  the  water  molecule  is  a  b-type  asymmetric  rotor,  selection  rules 
allow  only  ee—oo  and  eo-«oe  transitions,  where  a  and  o  refer  to  the  evenness 
or  oddness  of  the  limiting  prolate  and  oblate  symmetric  rotor  quantum  num¬ 
bers.  This  is  reflected  in  the  spectrum  by  two  independent  sets  of  transi¬ 
tions.  The  eo  and  oe  sets  cannot  combine  with  the  oo  and  ee  sets.  Hence, 
the  eo  and  oe  sets  cannot  combine  with  the  ground  state,  which  is  ee. 

This  results  in  the  1  ^  level  being  the  lowest  state  of  the  eo  and  oe  levels. 

The  energy  difference  between  the  ground  state,  0qo,  energy  level  and  the  1  j 
energy  level  cannot  be  determined  directly.  Fortunately,  the  sum  rules  of 
Fraley  and  Rao  (Ref.  58)  make  possible  the  calculation  of  this  parameter. 

Fraley  and  Rao  (Ref.  58)  have  developed  a  set  of  sum  rules  for  ortho¬ 
rhombic  asymmetric  molecules  that  are  very  useful  in  the  analysis  of  the 
pure  rotational  spectrum  of  water  vapor.  These  sum  rules  are: 


,  ,.J  (ee  +  eo  -  oo  -  oe)  .  1  /r_o  .  _o  .  ,  t 

(-!)  - - JUTIJ - -  =  Ao  -  1  <Taaaa  +  Taa>  +  3<aaaaaa> 


+  J(J  +  1)  [£  T°aaa  -  3(aaaaaa)] 


+  J^(J  +  l)2(aaaaaa) 


(13) 


(-1) 


J  (ee 


+  oo  -  eo 

Jirriy 


oe) 


Bo-5!Tbbbb+Tbb>+?TbbbbJ(J+l>  <14> 


(-i  r 


(ee  +  oe  -  eo  -  oo) 

- JfJT  I'J -  " 


=  C  -  i 
o  3 


(T  + 

cccc 


)  +  4-  Tu  J(J+  1) 
cc  3  cccc 


(15) 
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(C  T/T  °0\U71TTT)  =  T  <A  +  B  +  C  )  -  X  (T°  +  T°  +  T°  +  T?, 

J(J  +  1)(2J  +  1)  3  o  o  o'  60  aaaa  aa  bbbb  bb 


+  Ttccc  +  Tcc>  +  IT  +  W  +  »  [w  <3Taaaa 


+  3T°  ,  +  3T°  -  2T°  -  2T°  -  2T°  ) 

bbbb  cccc  aa  bb  cc' 


-  ~  (aaaaaa)  j  +  J2(J  +  l)2  j  (aaaaaa) 


(16) 


The  notation  ee,  for  example,  stands  for  the  sum  of  all  energy  levels  of  a 

given  J  whose  limiting  prolate  and  oblate  symmetric  rotor  quantum  numbers 

are  both  of  even  parity.  A  ,  B  ,  and  C  are  the  rotational  constants  associ- 

o  o  o 

ated  with  each  of  the  principal  axes  of  the  molecule.  T  ,  T  ,  etc.  ,  are 

aaaa  aa 

the  quartic  distortion  constants,  (aaaaaa)  is  the  higher  order  distortion 
constant.  Similar  terms  for  (bbbbbb)  and  (cccccc)  would  be  added  to  Eqs.  (13) 
through  (15)  if  analysis  showed  them  to  be  necessary. 

Combinations  of  these  sum  rules  give  four  sets  of  four  simultaneous 
equations,  with  the  energy  difference  A  between  the  0qo  and  1qj  energy  levels 
as  one  of  the  unknowns.  This  was  done  for  the  transitions  presented  here 
giving  the  result  A  =  23.7913  +,  0.  0016  cm’  as  the  average  for  the  four 
sets  of  equations.  This  agrees  well  with  the  value  quoted  by  Gates,  et  al. 

(Ref.  57)  of  23.79  cm'*,  Having  this  value,  we  can  calculate  the  energy  of 
all  the  eo  and  oe  levels  relative  to  the  true  ground  state. 

Because  some  of  the  transitions  observed  here  were  analyzed  as  blends 
of  very  closely  spaced  lines  and,  thus,  are  unsuitable  in  the  determination 
of  energy  levels  and  since  it  was  desirable  to  extend  the  calculations  to 
J  =  7,  it  was  necessary  to  use  transitions  observed  by  other  experimenters. 
Thet  i  are  listed  in  Table  6.  Two  of  these  are  microwave  transitions  (Ref. 

59)  and  are  more  accurate  than  our  measurements;  seventeen  others  were 


-74- 


Table  6.  Other  Transitions  Used  to  Derive  Energy  Levels 


Frequency, 

-1 

cm 


Calculated 

Frequency,*  Assignment  Ref. 

cm"  * 


0.7416830 

0.74 

523"~616 

59 

12.68200 

12.  68 

3  —4 

J2 1  14 

59 

314.74 

314.74 

312  ~*441 

52 

340. 55 

340. 55 

726~~835 

5 

2 

349.77 

349.79 

661  ~~770 

349.79 

660~~77lf 

384.88 

384. 86 

615  ~~?44 

385.54 

385. 52 

414~*541 

396.44 

396. 43 

7  17  ~*826 

426. 33 

426. 31 

422  ~*551 

431.  16 

431.  14 

515~642 

441.75 

441.75 

423~~550 

442.09 

442.  10 

523~~652 

456.  87 

456. 88 

624~~753 

472  .  39 

472.41 

524~*651 

506.  94 

506.93 

625~*752 

536.25 

536.26 

532~661 

541.07 

541.07 

533~*660 

554.63 

554.64 

633~*762 

51 

567.23 

567.23 

634  "*761 

51 

♦Energy  levels.  Ref.  57 


taken  from  the  spectra  of  Rao,  et  al.  (Ref.  52)  and  Blaine,  et  al.  (Ref.  51) 
in  the  300  to  600  cm’*  region.  These  transitions  were  unassi^ned  in  the 
original  papers.  They  were  chosen  and  assigned  on  the  basis  of  frequency 
and  relative  intensity  calculations  (Ref.  56)  using  the  energy  levels  of 
Gates,  et  al.  (Ref.  57).  The  transitions  were  chosen  when  calculations  indi¬ 
cated  that  a  single,  relatively  intense  transition  should  occur  we  1  separated 
from  its  nearest  intense  neighbors.  As  can  be  seen  in  Table  6  the  lines 
chosen  agree  with  calculated  values  to  within  0.  02  err  the  quoted  accuracy 
of  the  other  spectra.  However,  since  the  accuracy  of  these  frequencies  Is 
lower  than  our  measurements,  the  energy  levels  derived  using  these  transi¬ 
tions  are  quoted  to  one  less  significant  figure. 

The  energy  levels  derived  from  these  calculations  are  given  in  Table  7, 
together  with  values  obtained  from  near-IR  vibrational -rotational  spectra 
(Refs.  57,  60,  61).  The  accuracy  of  levels  determined  solely  from  transitions 
measured  in  this  work  is  estimated  to  be  1.0.  003  cm"  .  These  energy  levels 
are  indicated  by  four  figures  after  the  decimal  point.  Energy  levels  which 
were  determined  using  transitions  from  higher  frequency  spectra  are  indi¬ 
cated  by  three  places  after  the  decimal  point  and  are  estimated  to  be  accurate 
to  1.0.  02  cm’*,  The  energies  of  the  6^-,,  6^,  6,  j,  and  6^q  appear  anomalous. 
Perhaps  there  was  an  error  accumulating  in  one  of  the  two  series  that  be¬ 
came  obvious  only  in  these  almost  degenerate  levels;  or  perhaps  one  of  the 
transitions  used  in  the  determination  of  these  levels  was  misassigned.  In 
general,  though,  the  agreement  in  the  values  of  the  energy  levels  between 
this  work  and  the  near-IR  data  is  good. 

Having  determined  a  complete  set  of  energy  levels  up  to  J  --  7,  we 
can  calculate  frequencies  of  other  transitions  which  were  not  used  in  the 
calculations.  Table  8  shows  a  comparison  of  observed  and  calculated  fre¬ 
quencies  for  three  microwave  transitions  (Ref.  59).  The  agreement  is 
within  the  quoted  accuracy  of  the  energy  levels.  Also,  in  Table  5,  a  similar 
comparison  can  be  made  for  FIR  transitions.  Again,  the  agreement  is  good. 
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Table  7.  (cont) 


Jk-1-K+1 

This  Work 

Fraley 
(Ref.  60) 

Benedict 
(Ref.  61) 

Gates 
(Ref.  57) 

431 

383.8372 

.821 

.83 

.85 

441 

488. 110 

.110 

.  16 

.  10 

4 

’40 

488. 1256 

.  136 

.  17 

.  13 

505 

325. 399 

.  328 

.  36 

.35 

515 

326.6203 

.601 

.  64 

.64 

514 

399.  50? 

.  440 

.45 

.46 

524 

416.  202j 

,174 

.  18 

.22 

523 

446.  56q 

.487 

.  50 

.50 

533 

503.  96lg 

.938 

.98 

.00 

532 

508. 86j 

.782 

.80 

.81 

542 

610. 1005 

.  108 

.  12 

.  12 

541 

610. 384 

.334 

.34 

.35 

551 

742. 10? 

.071 

.08 

.  10 

Wl  i 
O  1 

742. 11? 

.074 

.08 

.  10 

606 

446. 704j 

.671 

.66 

.71 

6.6 

447. 302 

.231 

.24 

.24 

615 

542.8909 

.873 

.87 

.91 

*25 

552. 96Q 

.880 

.92 

.92 

*24 

602. 762j 

.728 

.67 

.77 

*34 

649. 02? 

.952 

.97 

.97 

Sj 

661. 5343 

.514 

.  56 

.  56 

<-43 

756.  765 

.704 

.72 

.75 

*42 

757. 78Q 

.  77 

.  78 

*« 

888. 65q 

.602 

.61 

.  60 
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Table  7.  (cont) 


Jk-pk+i 

This  Wcrk 

Fraley 
(Ref.  60) 

Benedict 
(Ref.  61  ) 

Gates 
(Ref.  57) 

*51 

888. 592 

.635 

.65 

.  63 

661 

1045. 11 j 

.045 

.09 

.07 

660 

1045. 032 

.045 

.09 

.07 

707 

586. 295 

.213 

.28 

.26 

7 17 

586. 464 

.449 

.43 

.48 

716 

704. 262 

.194 

.  18 

.22 

726 

709.  594g 

.558 

.58 

.60 

725 

782. 45^ 

.383 

.41 

.40 

735 

816.670J 

.671 

.65 

.72 

734 

842. 394 

.334 

.36 

.38 

744 

927. 771 

.723 

.76 

.77 

743 

931. 272 

.228 

.23 

.22 

7  53 

1059. 632 

.62 

.65 

752 

1059. 900 

.824 

.83 

.85 

762 

1216. 164 

.  199 

.  19 

.20 

761 

1216. 25? 

.199 

.  19 

.20 

771 

1394, 802 

.815 

.91 

.86 

770 

1394. 88^ 

.815 

.91 

,  86 

827 

885. 642 

.569 

.62 

.  62 

826 

982. 904 

.92 

.91 

836 

1006. 14? 

.090 

.  14 

.  12 

835 

1050. 144 

.  20 

.  15 

®45 

1122.  743 

.694 

.78 

.72 

8« 

1131. 762 

.78 

.  76 
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Fig.  29. 


The  tom  rule 
The  quantity 
traded  from  iumi  to  n 


Filot  for  the  "a 
27. 70  -  r 


axis  of  the  water  m 
72  J  {J  +  1)  }  cm*  has  b< 
te  the  curvature  more  a 


Fi,.  32.  The  lum  rule  plot  for  the  combination  of  ell  uet  of  the 
water  molecule. 
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Table  8.  Comparison  of  Observed  and  Calculated  Microwave 
Transitions  of  Water  Vapor 


Observed 

(Ref._J9), 

cm 

Calculated  f 
cm"  1 

Assignment 

6. 114549 

6. 1133 

2  —  3 

*20  13 

10. 84596 

10. 8429 

4  —  5 

22  15 

14.  94370 

14.  946 

3  —  4 

■*30  23 

Also,  with  a  complete  set  of  energy  levels  the  sum  rules  of  Fraley 
and  Rao  can  used  to  determine  the  rotational  constants.  Plots  of  the  sum 
rules  are  shown  in  Figs.  29  to  32.  Only  the  plots  in  Figs.  29  and  32,  which 
are  associated  with  the  a-axis  of  rotation,  show  significant  curvature.  This 
indicates  that  the  data  can  give  a  significant  value  for  the  higher  order  dis¬ 
tortion  constants  only  for  the  a-axis  (aaaaaa). 

Least  squares  fits  to  the  sum  rules  shown  in  Fqs.  (13)  through  (16) 
resulted  in  the  values  shown  in  Table  9.  The  errors  quoted  are  95%  confidence 
intervals.  The  iollowing  formulas  (Ref.  58)  applicable  to  planar  orthorhombic 
molecules  are  used  to  obtain  the  remaining  quartic  distortion  constants  and 
the  rotational  constants: 
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Table  9-  Parameters  Derived  from  Least  Squares  Fits 
to  the  Sum  Rules  of  Fraley  and  Rao  (cm*') 


A  -iT°  -  27.8801  +  0.0034 
o  4  aa 


B  -ir®  *  14.  513410.  0090 
o  4  bb 


C  ~l‘£°  a  9.2813  1  0.0021 
O  4  cc 


T  =  -0.0084  ±0.0012 
aaaa 


T”  a  -0.  0083  ±0  0012 
bbbb 


T  =  -0.00107  ±0.00027 
cccc 


(aaaaaa)  =  (6.  89  ±  0.  26)  x  10 


(bbbbbb)  *»  (cccccc) »  0 


A  tB  K  =  51.6742  ±0.0027 
o  o  o 


T°  +  T°  +  T°  s  -0.  00829 
aa  bb  cc 


(aaaaaa)  +  (bbbbbb)  +  (cccccc)  =  (6.  02  ±  0.  98)  X  10 
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The  subscript  e  refers  to  equilibrium  values  of  the  rotational  constants  that 
are  related  to  the  ground  and  excited  state  constants  by  relations  of  the 
following  type: 


Av"  Ae  ‘  £  Qi<v  +  I >  (2D 

where  v  is  the  vibrational  quantum  number  and  sum  is  over  the  fv  ndamental 
vibrations  of  the  molecule.  The  a's  arid  equilibrium  rotational  constants 
were  evaluated  by  weighted  least  squares  fits  to  our  values  of  ground  state 
constants  plus  approximately  40  excited- state  constants  given  by  Benedict 
(Ref.  62).  The  results  of  these  analyses  are  given  in  Table  10.  With  these 
values  of  the  equilibrium  rotational  constants,  the  remaining  quartic  dis¬ 
tortion  constants  and  the  ground  state  rotational  constants  can  be  calculated. 
The  results  of  these  calculations  are  given  in  Table  11.  The  errors  quoted 
are  95%  confidence  intervals.  The  values  of  the  ground  state  rotational 
constants  (Ref.  62)  as  determined  from  the  vibrational- rotational  spectrum 
are  also  included.  The  agreement  is  excellent.  Also  included  are  calculated 
values  of  the  six  quartic  distortion  constants.  These  were  calculated  from 
relations  given  by  Kivelson  and  Wilson  (Ref.  63)  between  vibrational  force 
constants  and  the  distortion  constants.  The  force  constants  were  determined 
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Table  10.  Equilibrium  Rotational  Constants 
(cm"1) 


A  =27.110? 
e 

B  =  14,  5922 
e 

Ce  =  9.5120 


A 

B 

C 

Q1 

0. 8168 

0.  2075 

U, 1694 

q2 

-3.  4957 

-0.  1561 

0.  1387 

q3 

1.  1418 

0. 1096 

0.  1467 

from  the  harmonic  frequencies  of  H^Q  and  D^O  S*ven  by  Benedict,  et  al. 
(Ref.  64).  The  agreement  between  observed  and  calculated  values  of  the 
quartic  distortion  constants  is  quite  good.  It  should  be  noted  that  the  six 
values  are  not  linearly  independent.  They  are  linear  combinations  o£  four 
independent  values. 

5.  LINE  WIDTH 

a.  Direct  From  Interferogram 

As  was  mentioned  previously  (Sections  II  and  III-A)  the  fall-off  of 
intensity  in  an  interferogram  is  due  in  a  large  measure  to  the  width(s) 
of  the  lines  in  the  spectrum.  If  the  lines  of  the  spectrum  are  very  nearly 
evenly  spaced,  then  it  will  be  the  signatures  that  are  dampled  out  as 
-exp  -  2irx€  (Ref.  22). 

A  number  of  interferograms  (conditions  of  different  pressures)  were 
obtained  for  CO  in  a  15-cm  absorption  cell  (see  Fig.  19). 


Table  1 1 .  Ground  State  Rotational  Constants 
of  the  Water  Molecule 


This  Work 

Previous  Work  (Ref.  62) 

A  =  27.8761  ±0,0034 

0 

27. 877 

Bq  =  14.  5074  1  0,  0090 

14.  512 

C  =  9. 2877  £  0. 0021 
o 

9.  285 

This  Work 

Calculated 

T°  =-0.1084  10.0012 

aaaa 

-0. 1088 

Tbbbb  =  ~00083  ±0.0012 

-0.00672 

T°  =-0.00107  10.00027 

cccc 

-0.00093 

Taabb  =  +0- 01985 

+0.01846 

Ccc=-°- 00491 

-0.00515 

Tbbccs-°- 00108 

T°  =  -0.  01602 
aa 

Tbb=  '°  02368 

-0, 00077 

T°  =  ;0.  02584 
cc 

[aaaaaa]  =  (6.  89  1  0.  26)  X  10"5 
[bbbbbb]  »[ccccccj*  0 

NOTE:  Only  four  of  the  T”p  g  are  independent,  the  others  being  linear 
combinations  of  these. 


In  the  analysis  of  the  signa¬ 
ture*  ,  the  intensity  was  taken  as 
the  average  between  the  largest 
maxima  and  two  adjacent  minima . 

The  fir st through  the  ninth  signa¬ 
tures  were  used.  The  zeroth  sig¬ 
nature  was  not  used  because  any 
imbalance  between  the  reference 
beam  and  sample  beam  would  be 
most  evident  here  (i.  e.  ,  in  the 
neighborhood  of  x  s  .0,  zero  opti¬ 
cal  path  difference).  Signatures 
beyond  the  ninth  were  not  used 
for  reasons  discussed  later  in  this 
section.  The  resulting  intensities  were 
tensity  and  then  the  natural  logarithms 


Fig.  33.  The  natural  logarithm  of 
the  reduced  optical  path  difference 
for  CO  (1  =  15  cm,  p  =  521  Torr). 


normalized  to  the  first  signature  in- 
of  these  relative  intensities  were 
plotted  as  a  function  of  x.  (The 
x's  were  taken  where  the  largest 


Fig.  34.  Linewidth  parameter 
versus  pressure  for  CO. 


maximum  of  the  signature  occurred. ) 
The  resulting  graph  is  a  straight 
line,  within  the  experimental  error, 
and  its  slope  is  2 ire.  A  typical  re¬ 
sult  is  shown  in  Fig.  33. 

This  procedure  was  carried 
out  for  four  runs  at  pressures  of 
521,  440,  375.  and  254  Torr.  The 
€' s  obtained  were:  0.  065  £0.  004, 
0.057  1  0.003,  0.047  ±0.005,  and 
0.  033  £  0,  004  cm"  * ,  respectively. 
These  results  are  shown  in  Fig.  34, 
where  <  is  plotted  as  a  function  of 


pressure.  The  result  is  a  straight  line  whose  intercept  is  the  origin,  within 

-4  -1 

experimental  error.  The  slope  is  c/p  =  (1.2_  ±  0.  i£)  X  .  0  cm  /Torr. 

*  -4  •  1 

The  microwave  value  for  e/p  is  (1.17  4  0.  07)  X  10  cm  /Tor r  and  the 
microwave  data  (Ref.  65)  were  obtained  at  a  pressure  of  less  than  1  Torr. 
Although  the  two  results  are  in  agreement  within  the  experimental  error,  the 
discussion  to  follow  will  show  that  the  value  obtained  from  the  interferometry 
data  should  be  larger  than  the  microwave  value  and  indeed  represents  an 
average  upper  limit  for  this  parameter. 

In  general,  this  method  is  subject  to  the  following  systematic  errors: 

First,  the  signature  shape  is  distorted  because  CO  (or  any  other 
molecule  for  that  m.  tter)  is  not  a  rigid  rotor.  Thus,  its  pure  rotational 
spectrum  consists  of  lines  that  are  not  exa  ctly  equally  spaced  in  wavenumber. 
(If  the  rotational  lines  were  exactly  equally  spaced,  the  signatures  would  be 
symmetric  and  undistorted. )  The  lines  converge  slowly  due  to  centrifugal 
distortion  effects,  and  this  is  sufficient  to  distort  the  signatures  so  that  the 
fall-off  in  intensity  is  greater  than  is  expected  from  the  exp(-2irx«)  factor 
alone.  The  method  used  for  measuring  the  intensities  compensates  for  this 
distortion  to  some  extent;  ho  /ever,  comparison  between  calculated  signature 
shapes  including  and  neglecting  the  centrifugal  stretching  term  indicate  the 
intensities  of  the  higher  signatures  (Sth,  9th,  etc. )  were  underestimated  by 
'*5%  (Ref.  22).  The  distortion  of  the  signatures  becomes  more  pronounced 
as  one  goes  to  higher  signatures  at  larger  x. 

Second,  the  lines  are  as sumed  to  be  Lorentsian  with  a  Lambert- Beer's 
Law  for  absorption.  That  is  I  =  I  e"^v^  where  y{v)  is  the  Lorentx  line 
shape  function  and  L  is  the  absorption  path.  In  Ref.  22  it  is  shown  that, 
when  the  exponential  is  expanded  in  a  power  series,  the  effect  of  the  higher 
order  terms  (3rd  and  beyond)  yields  a  factor  that  will  decrease  the  intensity 
of  the  signatures  in  addition  to  the  exp(-2wx«)  factor.  Again  this  becomes 
more  important  as  one  goes  to  signatures  at  larger  x. 

Because  of  these  two  errors,  only  signatures  one  through  nine  were 
used  in  the  analysis,  that  is,  signatures  of  relatively  small  x. 
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Third,  in  general  the  lines  will  have  different  width*  because  of 
possible  rotational  resonances  (Ref.  66)  and  Doppler  broadening.  For  a 
temperature  of  295*K,  ^(J  =  0  —  i)  ~  0.  05  X  10~4  cm**  and  the  J  =  19  -•  2C 
transition,  it  will  be  ~1 .  0  X  10"4  cm"* .  Our  measurements  were  at  pressures 
corresponding  to  «  of  0.  03  to  0.  07  cm-*;  therefore,  the  Doppler  broadening 
is  negligible. 

Fourth,  there  is  the  problem  of  instrumental  broadening  (in  interferom¬ 
eters  due  to  finite  entrance  and  exit  apertures,  and  the  usual  optical  problems, 
parallelnes*  of  the  two  grating  faces,  etc. ).  The  corresponding  problem  for 
conventional  spectrometers  is  a  very  great  one  and  in  fact  extensive  "slit 
width"  corrections  must  be  made,  since  often  the  width  of  the  lines  under 
investigation  is  much  narrower  than  the  band  of  radiation  passed  by  the 
entrance  and  exit  slits.  From  the  data  presented  here  (i.  e. ,  the  fact  that 
the  c  vs  p  curve  extrapolates  through  the  origin  well  within  the  experimental 
error),  it  is  apparent  that  any  instrumental  broadening  is  within  the  experi¬ 
mental  error. 

The  result  obtained  may  therefore  be  considered  an  average  upper 
limit  for  the  line  width  parameters  of  CO,  as  all  the  major  systematic  errors 
tend  to  give  a  larger  «  value. 

b.  From  the  Spectrum 

The  preceding  treatment  is  valid  as  long  as  each  of  the  lines  have  the 
same  line  width  parameter  and  the  absorption  is  not  too  strong;  even  then  it 
is  applicable  only  to  diatomic  and  linear  molecules,  and  the  result  is  an 

Q 

average  value.  The  more  practical  and  interesting  cases  are  those  where 


3 - — 

Individual  lines  can  be  studied  by  using  a  numerical  filtering  technique  dis¬ 
cussed  by  Connes  (Ref.  67).  Essentially  this  involves  convolving  the  inter  - 
ferogram  (obtained  between,  say,  0  and  vc  cm"*)  with  the  mathematical 
function  that  is  the  Fourier  transform  of  an  ideal  filter  passing  all  the  radia¬ 
tion  between  vi  and  wg  (the  region  of  interest)  and  none  of  the  radiation  outside 
this  range.  The  result  of  this  convolution  is  the  interferogram  one  would 
have  obtained  if  the  input  radiation  had  been  physically  filtered  by  a  filter 
having  the  characteristics  of  the  ideal  filter.  In  this  way,  single  lines  may 
be  studied  and  the  necessary  correction  for  the  systematic  errors  may  be 
applied  more  easily. 
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the  maximum  absorption  is  strong  and  the  line  width*  do  vary.  The  follow¬ 
ing  is  a  discussion  of  some  preliminary  results  obtained  for  the  pure  rotational 
spectrum  of  CO. 

A  beam  of  radiation  traversing  a  length  I  of  an  absorbing  medium  will 
absorb  intensity  according  to  the  well-known  Lambert-Beer's  Law,  and  the 
intensity  passed  I(v)  can  be  written 


I(v)  =  I0(v)e^(v)i  (22) 

where  y  v)  i*  the  incident  intensity  and  yW  is  the  absorption  coefficient.  As 
discussed  in  Section  Ill,  a  two-beam  interferometer  operating  in  a  double¬ 
beam  differencing  modes  does  not  measure  l(v)  and  in  fact  gives  little 
information  about  Iq(v).  What  is  measured  is  y  v),  and  it  is 

yv)  =  Ic(v)  -  I(V)  =  Io(v)  [1  -  e’Y(v)/]  (23) 


Now,  aside  from  the  fact  that  little  is  known  about  IQ(  v)  from  a  double - 

beam  run,  the  measurement  of  intensities  in  the  FIR  is  fraught  with  other 

difficulties  (see  Section  V).  Thus,  we  shall  restrict  ourselves  to  a  discussion 

of  what  information  is  available  from  I^{v),  making  only  the  assumption  that 

o 

Io(v)  does  not  change  radically  over  small  intervals  of  v.  Differentiating 
IA(v)  with  respect  to  v  yields 


dIQ(v) 

dv 


[1 


e-\(v)/ 


<iyv) 

Tv~“ 


*  [-* 


m  this  respect,  it  has  been  our  experience  that  the  largest  variation  of 
Io(v)/cm~l  is  caused  by  the  channel  spectrum  from  the  lamp  envelope  of 
the  source  (see  Section  V). 
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The  values  of  v  where  this  equation  is  satisfied  are  where  dl^(v)/dv 
has  Its  maximum  and  minimum  values.  In  fact,  for  a  symmetrical  line  they 
are  equally  displaced  about  v  .  We  define 


for  those  two  values  of  v  which  satisfy  Eq.  (26).  The  quantity  Av^  is 

the  wavenumber  difference  between  v  and  where  dl.(v)/dv  has  it*  extreme 

4  o  A'  ' 

values.  Dividing  Eq.  (26)  by  (Av^)  and  letting  y  =  ( «/A  v^)  and  Q  =  At  /i, 
we  obtain 


y4  -  2y2(l  -  Q)  -  3  a  0 


Now  Q  is  the  maximum  absorption  coefficient  times  the  absorbing  path  length 
f  and,  in  general,  may  range  from  near  aero  to  very  large  numbers.  Since 
Q  cannot  be  negative  and  y  =  ‘/Av^  is  also  a  positive  quantity,  the  only  root 
of  Eq.  (27)  having  physical  significance  is 


1  /2, 1  fZ 


y  =  {1  -  Q  +  [(1  -  Q)  +  3J  '  } 


1/2  I/? 

For  Q  *  0,  the  weak  absorption  case,  y  =  3  or  <  =  3  A  v^,  which 

is  the  relation  usually  used  in  the  microwave  region  (Ref.  68).  For  Q  >  0^ 

.  ,1  /2  '•  ,  . 

•  •  *  av 

Now  Eq.  (26)can  be  solved  for  «  ifAv^,  A,  and  t  are  known,  -NowAv 
can  be  measured  (see  discussion  on.  the  CO  pure  rotational  spectrum  in  this 
section),  and  I  is  known  from  the  experiment.  The  parameter  A  is  a  function 
of  pressure,  pojxilation  of  the  states  involved  in  the  transition  (upper  state 
also  since  stimulated  emission  is  not  negligible  in  the  FIR,  a  thermodynamic 
equilibrium  is  assumed  however),  the  dipole  moment  of  the  molecule,  etc. 

An  explicit  formula  for  A  is  given  in  Townes  and  Schswlow  (Ref.  66).  Thus 


-92- 


if  this  information  is  known,  A  may  be  calculated.  If  this  is  not  the  case, 
then  it  is  still  possible  to  get  «  by  making  two  runs  at  different  values  of  i. 
Then  two  Eqs.  (26'  will  result  and  the  factor  involving  A  can  be  eliminated. 
The  resulting  equation  will  have  f  as  a  function  of  (Av^)^,  (Av and  the 
ratio  (I  j  /i 2'- 

In  practice  it  is  necessary  to  correct  the  measured  4v^,  since  in 
general  the  resolution  of  the  interferometer  is  finite  and,  thus,  the  line 
shape  observed  will  be  the  tru*  line  shape  convolved  with  the  scanning 
function  of  the  interferometer.  As  might  be  expected,  the  raw  4^ 

(L  finite)  is  larger  than  the  true  Av^(L  -*  *>,  Fox  AJ  f*  «i,  a  correction 
curve  of  the  form  (Av^)^/(Av vs  L t  can  be  plotted  and  is  shown  in 
Fig.  35.  Plotted  also  is  «^/«^  vs  L«  .  As  is  seen  this  quantity  changes 
very  slowly  with  L,  and  thus  a  reasonable  value  of  is  obtained  directly 
from  the  spectrum;  L  is  known  and,  therefore,  the  correction  factor  can 
be  determined. 

Such  a  correction  curve  is  inadequate  when  (AJ/*)  is  not  small 
compared  with  one.  We  have  calculated  correction  curves  for  various 
values  of  (AJ  /« ).  A  typical  case  is  shown  in  Fig.  36. 

Preliminary  results  for  the  pore  rotational  spectrum  of  CO  are  shown 
in  Fig.  37.  The  results  have  been  normalised  to  1  - xtm  pressure  (assuming 
*  varies  linearly  proportional  to  pressure).  The  solid  lines  are  curves 
from  a  theoretical  calculation  (Ref.  69)  of  self-broadened  line  widths  for 
linear  me le coles. 

Because  the  dipole  moment  of  CO  is  so  small  (0.  112  Debye),  it  is 
•the  qufcdrupole  moment  Q  that  is  responsible  (at  least  for  low  J)  for  the 
pr ensure  broadening.  Although  the  data  points  have  considerable  scatter, 
a  vt>lue  of  Q  somewhat  less  than  4X10  ^  esu  is  indicated. 

t>.  ABSORPTION  O F  GASES  IN  THF  MEDIUM  PRESSURE 

RANGES  (~l»i£  aim) 

The  preliminary  results  of  line  width  measurements  have  pointed  out 
the  desirability  of  making  measurements  at  higher  pressures.  In  the  low  to 
trtedium  pressure  range  (where  only  two  body  collisions  are  important) 
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slut  ion  for  weakly 
ariation  of  t  .  /«  . 


Correction  curve  for  (Avq)  due  to  finite  resolution  for 
medium  absorbing  lines  {Al.  A  =  3,2813). 


« Uw’/otm) 


studies  in  the  near-IR  (vibration-rotation  spectra)  and  microwave  (usually 
low  J  pure  rotational  spectra)  region,  it  has  been  found  that  the  absorption 
lines  are  broadened  linearly  with  pressure.  At  higher  pressures  (where 
three  or  more  body  collisions  become  important),  the  lines  are  broadened 
in  an  asymmetric  manner  (i.  e.,  the  line  shape  is  nc  longer  symmetric 
about  the  peak  absorption).  It  has  also  been  observed  that  the  lines  are 
shifted  in  wavenumber  because  of  the  increased  pressure.  Most  of  these 
studies  have  been  carried  out  in  the  near-IR  (Ref.  70). 

Before  the  development  of  two-beam  interferometers  for  the  FIR, 
such  investigations  were  not  feasible  since  the  S/N  was  so  poor  as  to  preclude 

using  narrow  enough  slits  in  conven¬ 
tional  spectrometers  to  make  studies 
on  line  shapes.  Usually  only  investi¬ 
gations  reporting  wavenumber  measure¬ 
ments  were  carried  out.  The  per¬ 
formance  of  the  Aerospace  FIR 
lamellar  grating  has  indicated  that 
some  of  these  studies  are  now 
feasible.  A  small  cell  has  been 
designed,  fabricated,  and  tested  in 
the  1  to  10  atm  range.  A  photograph 
of  this  cell  is  shown  in  Fig.  38.  It 
has  a  clear  circular  aperture  of  ~2.  5 
cm  diam  and  an  absorption  length  of 

oo  _*•  ,1  t  j.  ~6.  5  cm.  Some  preliminary  tests 

Fig.  38.  Absorption  cell  for  medium  r  7 

pressures  (~1  to  10  atm)  have  been  made  with  CO  in  the  absorb¬ 
ing  cell,  but  no  analysis  or  measure¬ 
ments  have  been  made. 

B.  OPTICAL  CONSTANTS  OF  SOLIDS 

The  index  of  refraction  of  optical  materials  is  commonly  measured 
by  a  variety  of  techniques  for  wavelengths  in  the  visible  and  near  visible 
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regions.  In  the  FIR  spectral  region  (80  to  500  p,  20  to  125  cm"*)  most  of 
these  methods  are  not  used  either  because  they  are  not  feasible  in  this 
spectral  region  or  because  they  require  special  equipment  that  has  not 
been  developed.  An  exception  is  the  channel  spectrum  technique  in  which 
the  radiation  is  multiply  reflected  between  the  faces  of  a  flat  sample  so 
that  interference  in  the  transmitted  beam  produces  a  spectrum  consisting 
of  a  series  of  maxima  nearly  equally  spaced  in  frequency  ('  channels"). 

The  maxima  occur  at  those  wavelengths  for  which  the  optical  thickness  of 
the  sample  is  an  integral  number  of  half  wavelengths;  hence,  from  wave¬ 
length  and  thickness  measurements,  the  index  of  refraction  can  be 
determined. 

Interferometers  are  admirably  suited  to  this  determination  since  they 
are  fundamentally  wavelength -measuring  devices.  In  the  FIR  region,  they 
are  superior  to  grating  type  instruments  because  of  their  high  sensitivity 
and  spectral  resolution.  In  recent  years  the  required  instrumentation  and 
techniques  have  been  developed  to  a  high  degree  (Refs.  8,  9,  71,  72).  The 
application  of  the  Aerospace  Corporation  lamellar  grating  interferometer 
to  the  measurement  of  refractive  indices  is  discussed  here. 

The  transmittance  of  a  material  is  determined  by  its  index  of  refraction. 
This  index,  in  general,  must  be  complex  to  account  for  absorption  within  the 
material.  The  channel  spectrum  technique  is  effective  for  determining  the 
index,  provided  the  absorption  is  not  too  large.  The  energy  transmitted 
through  the  sample  can  be  calculated  in  the  following  manner  for  the  geometry 
of  Fig.  39.  Let  n  =  n  +  ik  be  the  complex  index,  r  =  (n  -  l)/(n  +  1 )  be  the 
amplitude  reflectance  and  t  =  2 /(n  +  1)  be  the  transmittance  at  each  surface. 
Then 


F’g.  39.  Geometry  of  sample  for 

channel  spectrum  index  of 
refraction  measurements. 


is  the  vector  sum  of  the  amplitudes  of 

the  transmitted  components.  Here 

v  is  the  wavenumber 

h  is  the  sample  thickness 

P  is  the  angle  between  normal  to 
the  i  ample  surface  and  the 
ray  within  the  medium,  and 

6  is  the  phase  shift  produced  by 
absorption  within  the  medium. 

The  energy  transmittance 

T(v)  =  A  •  A*  (30) 


has  been  derived  in  many  forms  (Ref.  21).  Here,  where  we  are  interested 
in  the  general  solution  with  absorbing  materials,  and  where  the  measurement 
is  mad*  by  Fourier  transform  methods,  it  is  convenient  to  express  this 
transmittance  by  (Ref.  8) 


T(v)  =t(v) 


p  cos  16 


) 


where 


t(v) 


4n 


(n  +  IV 


“exp( -crh/cos  P) 


1 


Ifete.2  ±  fc2 lex£l-?.h Iss*  P ) 

{l(n  +  ir+kZ)2(l-p2)> 


(31) 


(32) 
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(33) 


_ Kn  -  U2  ±  isLl - 

{[(n  +  1)  +  kjexp(ah  cos  P)} 


6  =  4irnhv  cos  P  +  6  =  4irnhv  cos  P  +  Tan 


‘0  77~. , ) 


(34) 


and  a  =  4trkv  is  the  customary  absorption  constant.  The  more  familiar 

2 

form  of  Eq.  (31)in  terms  of  the  same  parameters  is  T(v)  =  t(v)  (1  -  p  )/ 

2 

(1  +  p  -  2p  cos  6).  p  is  significant  in  that  it  is  the  fraction  of  the  energy 
that  survives  two  reflections  and  two  transmittance s  through  the  sample. 
The  angle  of  the  beam  within  the  sample  p  is  related  to  the  external 
angle  -y  by  Snell's  law,  which  becomes,  when  solved  for  cos  P  for  a 
sample  of  complex  index 


* 

v2 

1/2 

2fl  1 

cos  P  =  y 

f)  ♦(?) 

n 

\ 

» 

t(v)  is  the  average  transmittance  that  would  be  observed  if  there  were 
no  interference,  i.  e. ,  it  is  the  scalar  sum  of  the  transmitted  energies. 

Equation  (33)  is  obviously  the  Fourier  cosine  series  expansion  of 
T(v).  i.e.  . 

co 

T( v)  =  A |  cos  18 

i=0 
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where 


Aq  =  t(v)  ,  Af*0  = 

Because  the  interferometer  used  for  these  measurements  generates  an 
"interferogram"  that  is  the  Fourier  transform  of  T(v)  (multiplied  by  the 
source  function),  Eq.  (33)  is  a  very  convenient  form  to  work  with.  The 
interferogram  function  (see  Section  II)  can  be  written  as 

00 

F(x)  =  J"  Tf  v)(l  +  2^  Dp*  cos  i0j  cos  2pvxdv 
o 

where  x  is  the  optical  path  difference  in  the  interferometer.  In  this  equation 
the  source  function  has  been  absorbed  into  t(v)  and  the  constant  term  in  the 
interferogram  has  been  dropped.  If  t(v)  and  p  were  constant  and  v  extended 
from  zero  to  infinity,  F(x)  would  simply  consist  of  a  series  of  6-functions 
with  coefficients  A|  located  at  x  =  18/2*.  For  a  finite  range  of  6,  these 
6  functions  are  broadened  into  sine  (=  sin  x/x)  functions,  which  are  further 
distorted  by  the  nonconstancy  of  t  and  p. 

These  distorted  sine  functions  have  been  given  the  name  of  "signatures". 
Since  they  contain  all  the  available  information  about  the  sample,  **  it  is 


Here  we  note  that  the  modular  resolution  definition  is  the  appropriate  one, 
since  no  "channels"  in  the  observed  spectrum  will  be  seen  unless  at  least 
the  first  signature  is  observed. 
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possible  to  replace  the  interferogram  between  signatures  by  a  constant,  thus 
eliminating  a  major  fraction  of  the  recorded  noise  without  significantly  dis¬ 
torting  the  spectrum  obtained  from  the  transformed  interferogram. 

From  Eq.  (31)  it  is  evident  T(v)  will  be  a  maximum  when  all  the  cosine 
terms  in  the  sum  are  in  phase;  i.  e.  ,  when  0  =  m  •  2ir,  where  m  is  an  integer, 
from  which 


2  co3  (3  +  ~  m  (36) 

This  leads  to  a  series  of  maxima  equally  spaced  in  v,  provided  n  is  constant. 
From  measured  values  of  these  vmax  and  the  other  parameters  in  Eq.  (36), 
it  is  possible  to  calculate  n(v).  The  phase  shift  6  can  be  neglected  for  the 
first  approximation. 

From  a  measured  mean  transmittance  and  the  above  index  value,  the 
absorption  coefficient  a  ~  4irvk  can  be  calculated  (see  Eq.  (32),  which  in  turn 
can  be  used  to  determine  6  (see  Eq.  (34)  so  that  a  corrected  n  can  be  found). 
The  process  can  be  repeated  if  necessary  until  n  and  k  achieve  stationary 
values. 

This  method  is  limited  by  the  requirement  that  the  sample  be  trans¬ 
parent  enough  that  a  significant  fraction  of  the  radiation  can  survive  at  least 
two  internal  reflections  and  three  passes  through  the  sample.  We  now  make 
some  semiquc'.ntitative  estimates  of  the  maximum  thickness  permitted  for  a 

given  n  and  absorption  coefficient.  For  a  difference  spectrum  ( see  Section  III-C) 

.3 

the  noise  can  be  as  low  as  10  E  ,  where  E  is  the  scarce  intensity.  For 

00 

n  <  4  the  amplitude  of  the  channels  may  be  approximated  by  2t(v)oEq,  the 

first  term  in  the  series  of  Eq.  (31),  which  is  the  contribution  of  the  first  pair 

of  internal  reflections.  (The  r(v)  term  is  removed  in  the  differencing 

technique.  )  Assuming  a  S/N  of  at  least  2  is  needed  to  make  meaninglul 

.  -j 

measurements,  then  2t(v)p  >  2  X  10  ",  For  the  absorption  constants  and 
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frequencies  involved,  k  «  n,  and  an  adequate  approximation  becomes 


e-20.h<  10-3jn_4  IjV  tl) 

n(n  -  1 } 

Figure  40  gives,  as  a  function  of  the  absorption  constant,  a,  the  maximum 
value  of  h  which  satisfies  this  relation,  with  the  refractive  index  as  a 
parameter. 

We  now  turn  to  an  analysis  of  the  accuracy  with  which  the  parameters 
needed  for  calculating  n  can  be  determined. 


Fig.  40.  Maximum  sample  thickness  usable  with 
present  techniques  for  channel  spectrum 
determination  of  refractive  index  in  the 
FIR,  for  values  of  the  real  index  n  and 
absorption  coefficient  a  =  4irkv. 
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Since  the  order  number,  m,  must  be  an  integer  there  is  no  difficulty 
determining  it  exactly,  provided  that  the  accuracy  is  sufficient  to  distinguish 
between  adjacent  integers.  Since  in  this  wavelength  region  m  is  never  laiger 
than  a  hundred  or  so,  accuracy  is  easily  achieved.  A  simple  procedure  for 
determination  of  m  is  the  following:  Two  values  of  v  are  determined 
from  the  spectrum  for  m  =  rrij  and  m  =  nij  +  p.  Substituting  twice  in 
Eq.  (36)  and  subtracting,  we  have  2nh  cos  0(v.  -  v,,  )  =  p,  from 

which  an  average  n  can  be  found.  Inserting  this  value  in  one  of  the  above 
equations  then  leads  to  a  value  of  m.  This  technique  will  be  adequate  only 
as  long  as  n  is  relatively  constant,  but  when  it  is  not,  there  will  be  strong 
absorptions  and,  as  shown  above,  the  method  becomes  unsatisfactory. 

The  determination  of  frequency  furnishes  an  ultimate  limit  to  the  final 
precision.  Measurements  with  our  interferometer  on  pure  rotational  spectra 

indicate  a  limiting  instrumental  wavenumber  accuracy  of  about  0.  003  cm  *, 

- 1  5 

or,  for  60  cm  radiation,  about  5  parts  in  10  . 

To  achieve  this  precision  requires,  however,  locating  the  peaks  in  the 
spectrum  with  a  similar  precision.  The  precision  with  which  the  experi¬ 
mental  peaks  can  be  located  is  in  turn  limited  by  the  inherent  sharpness  of 
the  peaks,  governed  by  how  many  terms  in  the  series  of  Eq.  (31)  contribute 
significantly  to  the  spectrum.  The  rate  at  which  the  terms  in  this  series 
die  off  is  in  turn  controlled  by  the  index  of  the  material.  The  precision  is 
also  limited  by  changes  in  the  base  line, which  may  make  the  lines  appear  to 
shift,and  by  noise.  Figure  41,  the  channel  spectrum  of  fused  quarts,  dhows 

two  problems:  First,  if  there  is  any 
change  in  absorption  there  will  be  a 
changing  base  line;  second,  the  channel 
spectrum  wc  are  interested  in  uses  only 
a  small  part  of  the  dynamic  range  of  the 
instrument.  Both  of  these  problems 
may  be  significantly  decreased  by  means 
of  the  facilities  in  the  Aerospace  inter¬ 
ferometer  for  double-beam  difference 

transmittance  change  caused 
by  changi^r*  absorption. 
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Fig.  41.  Transmittance  of  2.  14  mm 
thick  fused  quartz  T(v) 
showing  both  channels  and 
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interferometry.  We  can,  while  obtaining  the  interferogram,  subtract  the 
interferogram  of  a  similar  sample,  different  only  in  that  it  is  wedged  and 
therefore  does  not  exhibit  a  channel  spectrum.  The  resulting  interferogram 
and  its  Fourier  transform,  the  spectrum,  appear  in  the  upper  pa  *t  of  Fig.  42. 

The  changing  base  line  has  indeed  been 
removed,  but  apparently  at  the  expense 
of  adding  a  good  deal  of  noise.  However 
because  of  the  localized  nature  of  the 
information  ne'»d<-  lor  Hie  channel 
spectrum  in  the  interferogi  *.m  signa¬ 
ture,  ar.  discussed  s.bc> .  . ,  it  Is  possible 
to  artificial. y  remove  much  of  the  noise 
from  the  incerferogram  before  it  is 
transformed  and  still  not  harm  the 
information  content  of  the  spectrum 
far  as  determining  refractive  indices 
i-S  concerned.  The  results  of  this  ore 
.mown  in  the  lower  part  of  Fig.  42. 

With  all  these  procedures  for  improving 
the  measurement  of  the  individual  maxima,  it  is  possible  to  decrease  the 
scatter  in  points  on  the  index  curve  by  -is  much  as  a  factor  of  4  or  5  in  the 
case  of  ar.  absorbing  material  of  low  index  such  as  fu$&d  quarts. 

If  the  refractive  indices  are  to  be  determined  with  an  accuracy  limited 
only  by  the  frequency  measurement,  ihc  error  in  the  thickness  value  must  not 
be  limiting.  For  a  typical  sample  thickness  of  3  mm,  this  means  a  thickness 
determination  to  within  a  fraction  o'  a  fringe  of  visible  light  The  flatness  and 
parallelism  of  the  faces  must  be  held  within  similar  tolerances.  Sample 
preparation  to  these  tolerances  has  been  carried  out  in  our  optical  shop.  The 
measurements  of  thickness  to  within  ±0,  2  u  were  made  by  a  commercial 
metrology  laboratory.  If  measurements  to  this  precision  are  to  be  meaningful. 

1  ^Metrolonics  Inc.  ,  Burbank,  Calif 


mum 


Fig.  42.  Difference  interferogram 
and  spectrum  for  fused 
quarts,  without  editing 
(upper)  and  with  editing 
( lower). 


the. material  must  be  homogeneous  to  a  similar  extent.  Since  independent 

methods  of  analyzing  the  homogeneity  are  difficult,  we  have  started  with  high 
,13 

purity  material. 

Within  the  sample  all  libht  rays  should  see  the  same  thickness  of 
material.  This  effective  thickness  depends  however  on  the  angle  with  which 
the  light  is  incident  on  the  sample,  the  rays  at  large  angles  seeing  an 
apparently  foreshortened  thickness  as  indicated  by  the  cosine  (3  in  Eq.  (34). 
Ideally  then  the  measurements  should  be  made  in  parallel  light,  but  since  this 
is  not  possible  with  the  Aerospace  instrument,  the  observed  T(v)  is  a  sum 
over  the  possible  angles  weighted  according  to  the  solid  angle,  and  is,  except 
for  a  constant  multiplier. 


T(v)  =  t(  v ) 


OP 

1  +  2  / ain  ?  co>  c6<p)de 


where  ±s  the  maximum  internal  angle.  Carrying  out  the  integration. 


f  \ 

/COS0  +  1  '* 

1  sin  p  cos  1 6(P)dp  =  j 

cos  I  4*.nh|  — - — ? - j 

\  f 

ain  y 


where  v  =  2*lnhv{co8  P  -  1). 

m 

Since  the  position  of  the  channels  <s  determined  by  the  first  term  in  the 
summation,  the  sin(y)/y  term  in  the  above  equation  with  m  =  1  is  the  one  th.»r 
attracts  our  attention.  In  nur  instrument,  with  its  cutoff  frequency  of  125  cm 
the  mi-xirnum.  value  y  can  have  is  250  ehn(cos  0  -  1).  If  we  wish  to  have 

y  <  *,  leaning  our  channel  amplitude  never  goes  Through  zero,  then 


Germanium  and  silicon  c-me  from  Exotic  Materials.  Costa  Mesa,  Calif. 
Fused  quartz  was  Infra  sil  grade  low  H_>G  quartz 
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hn  <  [250(cos  -  1)]  For  an  '  Eternal  half  angle  oi'  ?,  5  deg,  this  leads 

to  a  maximum  thickness  for  samples  of  silicon  and  germanium,  i.  i  and  1. 8 
cm,  respectively.  If  this  criterion  ie  met,  the  additional  effects  introduced 
bv  tbs  convergence  of  the  beam  on  the  amplitude  of  the  modulation  is  less 
pronounced  than  the  change  in  source  intensity  with  energy.  The  only 
correction  needed  then  is  to  replace  the  thickness  h  with  [h{l  +  cos  {3  )]  /2. 

A  sufficiently  accurate  determination  of  n  to  correct  for  beam  convergence 
can  usually  be  obtained  from  the  positions  of  the  signatures  in  the 
interferogram. 

For  materials  sufficiently  transparent  frr  this  method  to  be  successful, 
the  imaginary  part  of  the  dielectric  constant  k,  which  gives  rise  to  both 
absorption  and  the  phae?  factor  6,  must  be  small.  An  estimate  of  sufficient 
precision  can  therefore  be  made  by  means  of  Eq.  (32)  from  a  spectrum  of  the 
average  transmittance  r(v),  using  preliminary  channel  measurements  to 
obtain  a  value  of  n.  The  average  transmittance  spectrum  is  also  easily 
obtained  from  a  single -beam  interferogram  by  using  only  that  portion  of  the 
interferogram  around  x  -  0  and  then  dividing  this  spectrum  by  a  spectrum 

obtained  under  similar  conditions, 
but  without  the  sample  in  the  beam. 

For  the  materials  reported  on  here, 
the  phase  factor  was  small  enough  to 
ignore  in  comparison  to  the  uncertainty 
in  locating  the  channel  spectrum  peaks. 

Figure  43  shows  the  results  of 
the  analysis  as  described  when  applied 
to  the  fused  quartz  spectrum  in  Figs. 

4i  and  42.  The  error  bars  on  the 
refractive  index  data  indicate  the 
statistical  scatter  between  nearby 
points  and  are  smaller  near  the  middle 
of  the  spectral  region  of  interest, 
since  this  is  the  region  where  our 


Fig.  43.  Absorption  coefficient 
o  =  4rrkv,  and  real  re¬ 
fractive  index  for  fused 
quartz. 
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instrument  has  the  best  S/N.  The  absorption  curve  was  calculated,  by  means 
of  £q.  *32),  from  an  average  transmittance  spectrum.  The  error  bars  on  the 
absorption  curve  indicate  the  uncertainty  of  each  point  on  the  curve.  This 
uncertainty  arises  almost  entirely  from  the  uncertainty  in  the  transmittance 
measurements. 

Because  silicon  has  a  higher  index  and  lower  absorption,  more  internal 
reflections  contribute  to  the  spectrum,  as  may  be  seen  from  the  presence  of 
3  signatures  in  Fig.  44  Figure  45  shows  the  refractive  index  and  absorption 


Fig.  45.  Absorption  coefficient  and 
real  refractive  index  for 
silicon. 


constant  in  the  FXH  for  silicon.  The 
error  bars  have  the  same  significance 
as  those  used  on  the  fused  quartz  data. 
The  higher  number  of  internal  reflec¬ 
tions  making  significant  contributions 
to  the  spectrum  makes  the  peaks  in  the 
spectrum  sharper,  leading  to  lower 
scatter  in  the  index  data  points.  The 
transmission  of  one  sample  of  silicon 
has  been  determined  independently  with 
a  commercial  grating  spectrometer 
(Perkin-Elmer  301)  over  a  wide  spectral 


range.  These  results  are  presented  in  Fig.  46,  together  with  the  measured 
transmittance  with  the  interferometer  in  the  region  where  the  instruments 
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overlap.  F i  om  these  measurements, 
estimates  of  the  uncertainty  in  trans¬ 
mittance  have  been  made. 

Figure  47  shows  the  refractive 
index  and  absorption  constant  for 
germanium  in  the  FIR.  Figure  48  is  a 
plot  of  transmittance,  again  taken  over 
a  wide  range  with  the  Perkin-Elmer  301, 
for  comparison.  Because  of  the  chemical 
similarities  between  silicon  and  ger¬ 
manium,  a  spectrum  similar  to  Fig.  44 
for  silicon  might  be  expected.  Germanium,  however,  exhibits  more  dispersion 
and  more  absorption.  This  may  be  caused  by  impurities  in  the  sample  but, 
we  think  we  are  observing  the  low-frequency  side  of  the  band  at  100  cm"* 
reported  by  Aronson  and  McLinden  (Ref.  73). 


Fig.  47.  Absorption  coefficient  and 
real  refractive  index  for 
germanium. 
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We  have  presented  a  method  for  determining  the  refractive  index  of 
materials  from  channel  spectrum  measurements.  The  equations  used  to 
obtain  values  make  no  assumptions  except  that  of  nearly  normal  incidence,  so 
the  polarization  of  the  radiation  is  unimportant.  This  method  is  particularly 
well  suited  for  interferometric  measurement  techniques  since,  as  is  shown, 
the  spectrum  maybe  written  as  a  Fourier  cosine  series.  Since  the  inter¬ 
ferometer  gives  us  the  cosine  transform  of  the  spectrum,  the  approximations 
which  are  justifiable  in  treating  the  data  are  easily  seen. 

The  method  is  good  for  materials  which  are  not  strongly  absorbing  in 
this  spectral  region.  Limiting  sample  thicknesses  for  varying  combinations 
of  real  index  and  absorption  have  been  estimated. 

These  methods  have  been  applied  to  a  number  of  materials,  with  data 
presented  for  three,  to  show  the  application  of  the  technique.  It  should  be 
noted  that  the  values  are  close  to  these  measured  in  the  radio  frequency  part 
of  the  spectrum  and  may  differ  widely  from  values  in  the  visible  region,  as 
the  fused  quartz  shows. 

C.  INDEX  OF  REFRACTION  OF  GASES 

The  results  just  described  have  encouraged  us  to  try  to  apply  the  same 
principles  to  the  measurement  of  the  indices  of  refraction  of  gases  in  the 
FIR.  The  high-pressure  absorption  cell  can  be  modified  into  an  "interference 
cell.  "  The  faces  against  which  the  cell  windows  are  mounted  have  been 
ground  flat  (~3  fringes  of  HG  5460  A)  and  parallel  (~8  arc  sec).  The  windows 
are  phano  concave  lenses  of  crystalline  quartz.  The  focal  lengths  of  the 
lenses  are  chosen  such  that  the  entering  beam  is  made  parallel  inside  the 
cell  and  then  restored  to  its  original  divergence  on  leaving  the  cell.  The 
interior  surfaces  of  the  windows  are  plain  and  provide  the  reflective  surfaces 
for  interference.  In  the  FIR  the  reflectivity  of  a  single  surface  of  quartz  is 
~  13%  (in  this  respect  germanium  would  be  a  better  window  material).  Only 
preliminary  tests  of  the  cell  as  an  interferometer  have  been  completed.  The 
first  channel  signature  (at  x  x  13-cm  optical  path  difference)  has  been  observed. 
The  observed  S/N  for  the  signature  is  not  as  large  as  hoped  for,  and  part  of 
the  reason  for  this  is  that  the  source  is  extended  in  nature.  This  results  in 
difficulty  in  locating  the  focus  in  the  sample  chamber  and  subsequent 
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placement  of  the  interference  cell.  A  short  arc  Hg  lamp  is  presently  being 
tested  (for  reasons  in  addition  to  its  more  point-like  feature)  and  evaluated 
for  use  as  a  source  in  the  FIR. 

D.  NARROW  BANDPASS  FILTERS 

Fabry-Perot  etalons  consisting  of  two  parallel  wire  meshes  have  been 
used  by  earlier  investigators  (Refs.  74-77)  as  efficient  narrow  band  filters 
for  the  FIR  and  microwave  regions  (wavelengths  greater  than  about  50  p). 
These  filters  have  been  further  developed  in  this  laboratory,  the  fabrication 
procedures  and  the  performance  obtained  being  the  subject  here.  This  work 
has  been  greatly  facilitated  by  the  availability  of  the  Aerospace  Corporation 
lamellar  grating  interferometer,  which  has  a  sensitivity  and  spectral  resolu¬ 
tion  unavailable  to  the  earlier  workers.  This  permitted  much  better  evalua¬ 
tion  of  the  improvements  in  the  filter  design. 

The  spectrum  transmitted  by  such  a  filter  is  described  in  texts  on 
physical  optics.  The  power  transmissivity  T  is  given  by 

t2 

T= - r— 5 - 5 -  (37) 

[  ( 1  -  r j  +  4r  sin  6/s] 

where  t  and  r  are  the  power  transmissivity  and  reflectivity,  respectively, 
of  the  individual  meshes,  and  6  is  the  phase  shift  of  the  radiation  passing 
twice  through  the  gap  between  the  meshes,  i.  e.  , 

6  =  •  2vbp  cos  9  +  20  (38) 

where  v  -  i/\  is  the  wavenumber  of  the  radiation,  b  is  the  mesh  separation, 

P  is  the  index  of  refraction  of  the  assumed  nonabsorbing  material  in  the  gap, 

9  is  the  angle  of  incidence,  and  0  is  the  phase  shift  caused  by  reflection  at  a 
mesh  surface. 
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All  the  filters  discussed  will  have  an  air  (or  vacuum)  gap  between  the 
meshes  and,  hence,  n  =  1.  If  the  meshes  are  made  of  a  metal  with  high 
conductivity  the  losses  are  small,  or,  t  as  i  -  r;  and  thus 


T  = 


t 


i 


t  +  4r  sin^6/2  1  +  (4r/t^)  sin  6/2 


(39) 


It  i3  evident  from  this  equation  that  T  =1  for  5  =  m  2tr,  where  m  is  an 

M  max 

integer,  i.  e.  ,  there  is  a  series  of  transmittance  peaks  equally  spaced  in 

wavenumber  for  v  =  v  ,  where 

max 


mir  -  0 

max  ~  2irb  cos  9 


(40) 


Also,  for  6  ={m  +  1/2)  2v,  T  =  (1 
’  '  '  mm  ' 


r/1  +  r)  ,  therefore  T  -*•  0  as 
'  mm 


r  -*  1.  The  conventional  measure  of  the  sharpness  of  the  fringes  is  the 
finespe  F,  defined  as  the  ratio  of  the  separation  between  adjacent  fringes 
to  the  width  at  half  maximum  (half  width)  of  a  given  fringe.  For  r  near  1 
this  is  given  very  closely  by 


,  ,1/2  ,  ,1/2 
F  =  HilL—  =  _ 

1  -  r  t 


(41) 


and  F  -*  ao  as  r  -*  1.  A  large  finesse,  of  course,  indicates  narrow  fringes. 

The  above  theory  assumes  an  ideal  filter,  i.  e.  ,  no  resistive  losses, 
perfectly  flat  and  parallel  meshes,  and  the  angle  of  incidence  6  is  either 
zero  or  a  constant.  In  order  to  achieve  a  finesse  F  for  nearly  normal 
incidence,  the  allowable  variation  in  b  is  6b»i/2vF,  thus,  for  F  =  100, 
v~l00  cm"*.  6b  =  5  X  10"^  cm,  or  approximately  one  visible  light  fringe. 
This  result  is  surprising  in  view  of  the  general  impression  that  the  construc¬ 
tion  tolerances  for  FIR  components  can  be  rather  sloppy  because  of  the 
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large  wavelengths  involved.  Attempts  have  been  made  by  other  workers  to 
construct  tunable  filters,  with  adjustable  mesh  spacing,  but  the  performance 
has  been  mediocre  because  the  mesh  parallelism  was  not  adequately  main¬ 
tained.  Also,  simple  rugged  filters  have  been  constructed  by  cementing 
meshes  onto  the  faces  of  a  wafer  of  fused  quartz.  Again  the  performance  has 
been  mediocre,  perhaps  due  to  inhomogeneities  or  lack  of  flatness  in  the 
quartz,  but  more  probably  due  to  lack  of  perfect  contact  between  the  quartz 
and  the  mesh  caused  by  inclusions  of  cement. 


i 


i 


Fig.  49.  Details  of  filter  construc¬ 
tion.  The  apace  is  cut 
from  precision  stainless 
steel  shim  stock  and  the 
burrs  carefully  removed. 
The  faces  of  the  end 
pieces  touching  the  grids 
are  machined  and  lapped 
flat. 
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The  available  meshes  are  from 
5  to  2000  mesh  (rulings  per  inch),  and 
are  about  0.  0002  in. «  5-p  thick.  For 
the  FIR  region  considered  here 
(v  =  10  to  100  cm-*),  250-  to  1000-mesh 
rulings  are  most  suitable.  These 
meshes  are  fragile  and  have  no  rigidity 
of  their  own  and  so  must  be  carefully 
stretched  over  a  supporting  ring  to 
remain  flat,  as  shown  in  Fig.  49.  The 
meshes  are  cemented  as  smoothly  as 
possible  onto  the  two  temporary  sup¬ 
ports  shown.  No  wrinkles,  as  seen  by 
reflected  light,  can  be  permitted  at 
this  stage.  These  are  then  hung  over 
the  other  members,  as  shown,  the 
weight  of  the  support  rings  stretching 
the  meshes  taut.  The  assembly  is 
then  clamped  between  the  end  pieces, 
the  excess  mesh  trimmed  off,  and  the 
joint  sealed  with  epoxy  cement. 

T3 — : - - — 

Available  from  Buckbee-Mear s  Corp.  , 

St.  Paul,  Minn. 


Since  the  wire  meshes  must  be  very  accurately  flat  and  parallel,  the 
most  critical  component  is  the  stainless-steel  shim  space  shown.  Ours  have 
been  cut  from  precision  shim  stock,  any  burrs  on  the  edge  being  carefully 
removed  by  lapping.  The  end  pieces  also  are  carefully  machine  J  and  lapped. 

A  completed  filter  with  a  i-in.  aperture  is  shown  in  Fig.  50. 

Calculations  of  the  transmittance  01  'eflectance 
of  a  metal  mesh  are  difficult  because  of  the  complex 
two-dimensional  geometry.  The  most  successful 
approach  involves  the  calculation  of  an  equivalent 
transmission  line  and  has  been  solved  for  the  related 
problems  of  one-dimensional  arrays  of  various 
shaped  wires  (Ref.  78-80)  and  the  problem  of  circular 
openings  in  a  metal  sheet  (Ref.  81).  We  find  however 
a  very  satisfactory  analysis  of  filter  performance 
can  be  based  on  the  empirically  determined  values 
of  power  reflectance  r  and  power  transmittance  t 
Fig.  50.  for  a  single  mesh.  Scaling  the  distance  d  between 

pleted  filter.  mesh  opening  by  the  wavenumber  v  allows  us  to 
to  plot  results  in  terms  of  the  dimensionless 
variable  vd,  making  it  possible  to  compare  data  taken  on  different  meshes  in 
different  frequency  ranges.  Figure  51  is  a  plot  of  the  transmission  of  four 
different  mesh  materials  as  a  function  of  this  reduced  variable  vd.  The 
transmission  curves  for  this  region  are  nearly  straight  lines  on  a  loglog 
plot,  indicating  a  simple  power  law  relation.  The  curves  for  two  of  the 
materials  coincide  as  anticipated.  The  250-mesh  material,  however, 
exhibits  a  higher  transmittance,  as  might  have  been  expected,  since  each 
cell  in  this  mesh  has  about  17%  more  open  area  than  the  cells  in  the  other 
meshes.  The  nickel  mesh  exhibits  a  higher  transmission,  perhaps  due  to 
its  higher  resistivity. 

For  vd  greater  than  that  shown  in  the  figure,  the  measured  trans¬ 
mittance  first  rises  to  a  maximum  of  about  90%  at  vds  0.  85  and  then 
decreases  and  oscillates  for  greater  vd.  For  vd  >  1,  a  mesh  begins  to  act 


-11  5- 


TRANSMISSION 


Fig.  51.  Transmission  of  a  single  wire  mesh  as 
a  function  of  reduced  fequency,  vd, 
where  d  is  the  spacing  between  openings; 
□  250-mesh  (d  =  102p)  copper (Buckbee- 
Meers  ruling  no.  25027);  A  500-mesh 
(d  =  204p)  nickel  (ruling  no.  unknown); 

0  500-mesh  (d  =  204p)  copper  (ruling  no. 
509);  •  750-mesh  (d  =  3Q6p.)  v.«pper 
(ruling  no.  761). 


as  a  two-dimensional  diffraction  grating;  hence,  varying  amounts  of  trans¬ 
mitted  power  ar^  diffracted  out  of  the  acceptance  angle  of  the  measuring 
instrument,  making  details  of  the  measured  transmittance  dependent  on  the 
particular  experimental  geometry  used  for  the  measurement.  In  view  of  the 
variability  of  mesh  fabrication  and  the  difficulties  of  measuring  absolute 
transmissions  in  this  spectral  region,  we  feel  the  present  results  are  in 
satisfactory  agreement  with  those  of  others  (Refs.  9,  74-76,  82).  Indeed 
the  data  in  Fig.  8  of  Russell  and  Bell  (Refs.  9,  83)  superimpose  on  our 
500-  and  750-mesh  data  for  values  of  vd  from  0.  25  to  0,  75. 

The  Aerospace  interferometer  cannot  measure  reflectance  at  normal 
incidence,  and  hence  we  cannot  show  that  r  -  t-t,  thus  demonstrating  the 
absorption  in  the  mesh  to  be  small.  However  the  data  of  Mitsuishi, 
et  al.  (Ref.  75)  indicate  that  losses  are  indeed  small,  at  least  for  vd  <  0.  85. 
Also  the  performance  of  the  filters  discussed  in  Section  V  could  not  have 
been  achieved  had  there  been  large  losses  in  the  meshes. 
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The  choice  of  mesh  material  for  a  particular  filter  application  depends 
.a  the  wavenumber  and  finesse  specified.  For  example,  from  Eq,  (40)  if 
F  =  100  is  required,  r  as  97%  or  t  ^  3%.  From  Fig.  51,  this  gives  vd  =  0.  19 
for  the  500-  or  750-mesh  grids,  and  v  x,  54  cm"*  or  36  cm"1,  respectively. 
Note  that  the  above  finesse  assumes  an  ideal  filter  and  makes  no  allowance 
for  fabrication  imperfections. 

Filters  have  been  successfully  made  using  a  number  of  meshes  and  shim 
spacers,  A  measured  spectral  response  curve  fer  one  of  the  better  samples 
(750  mesh,  0.  0807-cm  space)  is  given  in  Fig.  52.  This  is  the  actual  energy 


Fig.  52.  Energy  throughput  of  a  wire  mesh  filter  consisting  of 

two  750-mesh  screens  separated  by  a  0.  08C5  cm  spacer. 
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throughput,  and  it  is  multiplied  by  the  spectrum  of  the  source  and  the 

admittance  of  the  interferometer.  Transmission  peaks  appear  approxi¬ 
mately  every  6  cm  ,  the  three  at  6,  12,  and  18  cm  being  missing 
because  of  the  low  energy  available  in  this  interval.  The  widths  of  the  peak3 
increase  with  v  as  expected  from  the  mesh  transmittance  given  in  Fig.  51. 
From  the  noise  level  at  the  base  of  the  peaks,  it  is  evident  that  peak 
transmittance  values  and  off-peak  rejection  can  be  determined  at  best  to 
within  about  5%.  This  is  poorer  than  the  measurements  on  single  grids 
because  of  the  small  amount  of  energy  transmitted  with  a  narrow  trans¬ 
mittance  peak. 

Wavenumber  determinations  are  much  more  accurate  than  intensity 
measurements  as  they  are  calculated  directly  from  the  grating  displace¬ 
ment.  Transmittance  peaks  can  be  read  to  about  0.  003  cm  *  and  half-widths 
to  about  0.  005  cm 

The  haH -width  thus  de*-~rmined  is  not  the  true  half-width  of  the  filter, 

but  it  is  broadened  by  the  finite  resolution  of  the  interferometer  and  by 

the  convergence  of  the.  measuring  beam  of  radiation.  The  resolution  of  the 

interferometer  is  limited  by  the  maximum  displacement  permitted  the 

lamellar  grating.  However,  before  this  limit  is  reached  the  interior ogr am 

signal  tails  into  the  noise,  and  hence,  in  practice  the  resolution  is  actually 

limited  by  the  S/N  In  order  that  there  be  a  usable  signal  it  is  necessary 

that  a  conveigent  beam  be  transmitted  through  the  {liter,  8  varying  from 

0  to  8  and  v  varying  from  m*  -  A/Z^b  to  me  -  cos  a,  according 

to  Eq.  (40)  For  a  beam  of  half-angle  8  it  can  be  shown  that  this  broadens 
^  ^max 

the  peaks  by  approximately  5v^  8^  ^‘■v/2. 

In  Fig.  5-3  finesse  calculated  irom  the  measured  half- widths  is  plotted 
against  v  for  thre^  750-mesh  filters.  As  expected,  F  decreases  with  a.  The 
largest  values  (100)  were  obtained  in  a  filter  with  a  thin  spacer  and  corre¬ 
spondingly  wide,  more  easily  measures  peaks.  The  two  curves  tor  the  filters 

with  the  wide  space  were  obtained  with  8  i-  deg  and  8  4  deg 

max  max  ' 

and  clearly  show  the  broadening  due  to  beam  convergence.  The  finesse 


values  reported  here  have  no  corrections  applied,  and  yet  they  are  about 
three  times  larger  than  values  reported  by  earlier  workers  (Refs.  74,  76) 
for  which  a  spectrometer  slit-width  correction  was  applied. 

The  magnitude  of  the  corrections  which  could  be  applied  here  is  indi¬ 
cated  by  Fig.  54  for  the  750-mesh,  0.  0807-cm  shim  filter.  The  upper  curve 


Fig.  53.  Measured  finesse  for 
750-mesh  filters  of 
different  spacing. 


Fig.  54.  Reduction  in  finesse  for 

750-mesh,  0.  081-cm  spac¬ 
ing  filter  due  to  limitations 
of  measuring  technique. 


is  the  finesse  calculated  from  the  measured  mesh  transmittance  (Fig.  51). 
It  assumes  a  perfect  filter  fabrication  and  ideal  measuring  instrumentation. 
Introducing  the  4  deg  beam  angle  convergence  lowers  the  finesse  to  the 
second  curve.  Including,  in  addition,  the  finite  resolution  of  the  particular 
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i.uerferometei  measurement  lowers  the  finesse  to  the  bottom  curve.  It  is 
seen  that  this  curve  is  in  close  agreement  with  the  actual  measured  finesse. 
Becatvse  the  corrections  are  an  order  of  magnitude  larger  than  the  ideal 
half-width,  the  true  finesse  cannot  be  calculated  with  any  precision.  We  can 
conclude  however  that  it  is  significantly  larger  than  the  value  of  58  measured 
at  24  cm  .  The  ideal  finesse  of  1000  shown  in  the  figure  is  unreasonably 
large,  as  it  would  require  meshes  flat  and  parallel  to  within  about  500  A 
Further  evidence  that  the  measured  finesse  is  too  low  is  given  by 
Fig.  55,  which  is  a  graph  of  the  measured  transmittance  at  the  peaks  for 

the  same  750-mesh,  0.  0805-cm  space 
filter,  measured  again  with  4  and  8 
deg  beam  angles.  Increasing  the  beam 
angle  and  hence  widening  the  peak 
decreases  the  measured  transmittance 
as  would  be  expected  from  an  unresolved 
peak.  The  peak  transmittance  at  this 
wavenumber  for  the  750-mesh, 

0.  0142-cm  spacer  filter  is  abotit  40%, 
or  about  three  times  as  large.  Since 
the  meshes  and  wavenumbers  are 
identical,  the  transmittance  should  not 
be  greatly  different.  This  would  indi¬ 
cate  a  true  finesse  of  at  least  180. 

The  off-band  rejection  cannot  be 
adequately  measured  with  our  instru¬ 
ment- -or  any  existing  FIR  instrumenta¬ 
tion.  From  the  data  presented,  however, 
we  have  some  confidence  that  an  adequate  mathematical  description  of  the 
filter  can  be  made  using  Eq.  (37)  and  the  measured  transmittance  of  the  single 
meshes  (Fig.  51).  Filter  imperfections  will  broaden  the  peaks,  but  they 
should  have  little  effect  on  the  off-band  rejection.  Resistive  losses 
however  will  degrade  the  rejection.  The  minimum  transmittance 


Fig.  55.  Transmission  as  a  function 
of  frequency  for  a  750- 
mesh,  0.  081 -cm  spacing 
filter  for  two  beam  con¬ 
vergence  angles. 


halfway  between  peaks  calculated  from  this  equation  is  shown  in  Fig.  56.  The 

low  frequency  peaks  (Fig.  53)  do  indeed 

intersect  sharply  with  the  noise  base  line, 

. .  T . T-  _ _ T  ...  , T _ T _ 7 _ so  these  results  are  not  unreasonable.  The 

high  rejection  levels  shown  in  Fig.  56  are 

of  course  to  some  extent  meaningless  in  a 

!  /  _  room  temperature  environment,  since  the 

S  blackbody  radiat'on  from  the  filter  itself, 

/  or  reflected  from  other  portions  of  the 

T«»  /  instrument  by  the  filter,  can  easily  be 

/  much  larger  than  the  oil-band  energy 

/  transmitted  by  such  a  filter. 

/  We  feel  that  these  interference  filters 

w4  -  / 

f  can  now  be  made  to  order  with  finesse, 

f  peak  transmittance,  and  off-band  rejection 

better  than  can  be  measured  with  any 
®o  to  to  to  to  no  fife  mo  existing  instrumentation.  In  particular 

■r  (cai'1) 

'  applications,  additional  filtering  will  be 

needed  to  eliminate  the  unwanted  higher 
Fig.  56.  Calculated  transmis-  ,  ,  .  ...  .  .  . 

sion  halfway  between  hatm0nlC‘'  ™*  filtering  will  h.v.  to  be 

peaks  of  a  750-mesh  tailored  for  each  case.  A  study  of  suitable 
filter,  materials  and  techniques  is  in  progress. 


<o  to 

»(«■•■) 


Fig.  56.  Calculated  transmis¬ 
sion  halfway  between 
peaks  of  a  750-mesh 
filter. 
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V.  NOISE,  EXPERIMENTAL  ERRORS,  AND  PROBLEMS 


Fig.  57.  Appearcnce  of  a  Lorentz  line  when  an  error  haa  been 

made  in  the  location  of  the  origin;  curve  A,  Pv  =  0, 

curve  B,  Pv  =  1/8,  and  curve  C,  Pv  =1/4.  ° 

o  o 
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f.gurc,  the  line  shape  can  be  severely  distorted  and  significant  apparent  line 
shifts  are  possible.  The  amount  of  distortion  depends  on  P  and  v  .  When  an 
error  of  this  type  is  made  in  the  interferogram,  its  effect  in  the  spectrum  is 
usually  quite  striking.  As  an  example,  assume  ft  ~  50p  (a  large  error)  in  an 
interferogram  made  up  of  lines  spread  out  in  the  wavenumber  region  0  to 
100  cm  .  In  the  region  ~0  cm  ,  the  lines  would  appear  almost  as  they 
should;  near  v  -  50  cm  * ,  where  Pvq  ~  1/4,  the  lines  would  have  the  most 
asymmetric  appearance,  somewhat  like  their  first  derivative  shape;  near 
v  *  100  cm”  ,  the  lines  would  again  have  near  their  true  shape  but  would  be 
upside  down.  It  is  obvious  that  good  spectral  data  will  be  obtained  only  if 
{3vq  «  1/4. 

When  ft  is  sufficiently  small  so  that  Zvftv^  is  small  enough  for  the 
small  angle  approximation  to  hold  in  Eq.  (42),  then  the  shifts  and  distortion 
in  the  line  shape  can  be  easily  calculated.  These  are  summarized  in  Table  12. 


Table  12. 


Quantity 


I(v) 


I( 


v) 

max 


at 


Distortions  Caused  by  Mislocation  of  Origin  (x  =  0) 


ft  =  0 


M0  (Bmall  angle  approx.) 


A* 


,  ,  2 

(V  -  V  )  +  € 

O 


At  /,  V  '  *0  \ 

- -j— jll  -  — j—q) 

V  -  V  )  +  <  \  / 


o 


.  ta 

>  2 


I(v) 


max 


A 

( 


Zt 

q  =  2itPvo 


t(‘  -  ^*) 

2*  (‘ +  ^f-) 


Half-width 


B. 


PERIODIC  ERROR 


Another  possible  error  in  the  measurement  of  x  is  a  periodic  error, 
a  type  of  error  expected  if  a  screw  type  mechanism  (Ref.  6)  (and,  as  we  shall 
show,  a  linear  inductosyn  system)  is  used.  Let  us  assume  a  monochromatic 
line  input. 

Input  Spectrum  I(v)  60'q) 

Then 


F(x)  =  cuo  2»txv 


Now  let  us  assume  a  periodic  error  of  period  L  ,  and  amplitude  B;  thus, 


x  ,  =  x.  .  4  B  sin 

measured  (true) 


:irx 


where  0  it  a  phase  factor  The  ^aturec  interferogram  function  will  then  be 


and  now,  calling 


e  e 

FJx)  -  cos  (2rxv  4  2r Bv'  sin(2"xv  4  *}] 
r)  O  O  C 

=  cot(2*xv  }cos|  2*Bi  sin(2*xv  ♦  #)] 
o  o  e 

-  *in{2«xv  )sin(  ZeBv'  sin(2*xv  4  *)) 
o  o  e  T  * 


We  now  assume  2i?Bv  is  small  enough  that  the  small  angle 
approximation  is  valid;  therefore, 

Fp(x)  ^  cos(2iTxv'r)  -  2ttBv  sin( 2tt-x vQ)sin( 2ttx e  -!  <f>) 

=  Cos(2'rx,v  )  -  ttBv  cosf  2ttx(v  -  v'  )  +  o]  +  ttBv  cosl  2ttx(v  +  v  )  4  <ftl 
»o  o  oe  o  oe 


In  the  spectrum,  the  first,  term  gives  the  line  at  v  ;  the  second  and 
third  terms  yield  ghosts  at  v  -  ±  v  ,  where  shape  will  in  general  be  dis¬ 

torted  by  the  phase  factor  <p  in  much  the  same  way  as  discussed  in  Section  V-A. 
It  is  obvious  that  if  2'n'Bv^  is  not  small  enough  that  the  small  angle  approxi¬ 
mation  is  valid,  so  that  higher  order  terms  should  be  carried  in  the  expansion, 

then  these  higher  order  terms  would  contribute  ghosts  at  r  ~  v  :fc  nv  We 

o  e 

can  estimate  the  first  relative  ghost  intensity: 


ghost  intensity  %  ^ 
line  intensity  o 


for  <p  -  0 


We  have  examined  our  linear  inductosyn  scale  for  this  type  of  error. 
The  following  technique  was  employed.  The  signature  of  a  diatomic  molecule 
should  be  evenly  spaced  in  optical  path  difference.  The  positions  of  the 
maximumg  of  the  first  eight  signatures  of  HBr  were  measured  (the  isotopic 
splitting  is  negligible  in  this  molecule  for  signatures  in  the  neighborhood  of 
the  origin).  These  eight  measurements  were  then  fitted  to  an  equation  of  the 
form:  x( signature  maximum)  =  an  +  b,  where  n  is  the  signature  number  and 
b  is  included  to  account  for  any  origin  error.  The  best  values  of  a  and  b 
were  determined,  and  the  deviation  calculated  and  plotted  as  a  function  of  x 
The  results  are  shown  in  Fig.  58  with  a  sine  curve  of  amplitude  1  and 
period  2  mm  fitted  to  best  cover  the  points.  The  error  bars  represent  the 
uncertainty,  due  to  noise,  in  locating  the  signature  maxima.  There  is  little 
doubt  that  a  periodic  error  exists.  For  3  —  1 .  c;?i.  and  a  period  of  —2  mm,  the 
ghost  intensity  at  ~50  cm  *  line  should  appear  with  '-2%  of  the  line  intensity 
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Fig.  58.  Optical  path  difference  error  in  a  linear  inductosyn 
scale.  The  superposed  curve  is  a  sine  curve  of 
amplitude  1.  5u  and  period  2  mm. 

at  displacements  of  ~  *5  cm  1  away  from  the  line.  The  over-all  noise  level 
of  most  of  our  spectra  would  hide  these  ghosts;  however,  see  Fig  52  .  There 
are  features  at  ~  *5  cm  *  about  each  of  the  '  lines  that  appear  to  be  ghosts. 

Tiie  remedy  for  this  error  is  of  course  a  more  accurate  measurement 
of  x,  the  optical  path  difference.  The  possibility  of  using  a  single  mode  laser 
is  being  examined. 

r  CHOICE  OF  Ax 

As  mentioned  previously,  the  optical  path  difference  increment  Ax  is 

fixed  b>  sampling  theory,  i.  e.  .  Ax  5  fZv  )  where  is  the  cvTofl  wave*' 

r  c 

number  of  the  spectrum.  Three  points  ar*  discussed:  (1)  determination  of  v  , 

c 

(2)  effect  of  choosing  Ax  too  large,  and  (3)  efs  ct  of  chopping  Ax  smaller  than 
need  he. 

We  determine  Ax(or  v  }  in  the  following  manner  An  interfere#  ram  oi 
the  source  is  obtained  with  a  relatively  small  Ax,  corresponding  to  a  large 
v  .  This  interferogram  is  then  transformed  and  the  spectrum  of  the  temp  a* 
modified  by  the  filtering  is  obtained  Where  the  intensity  go*.*  to  rrro 
(within  the  noise  level)  is  the  true  v  and  thus  the  maximum  Ax  is  determined 
This  is  done  whenever  sources  or  filtering  as  angemen**  «rc  chin#*d 
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When  Ax  is  rh-iM'D  Too  ia  rg<  .  a  :  •  . :  j.  '.<*«  s  plat  ••  Tins  is  best 

understood  in  terms  ot  the  s<  aiming  fum  tt->n  The  breadth  of  the  scanning 
function  is  determined  by  L,  the  maximum  optical  path  difference  generated. 
The  range  over  which  the  scanning  function  stans'  is  determined  by  Ax. 
Actually,  except  as  Ax  -  0,  there  is  more  than  one  scanning  function.  For  a 
given  Ax  the  first  scanning  function  covers  the  spectral  region  0  -  v  .  the 
second,  Zv  *  v  ,  the  third,  2e  •»  3e  ;  etc .  Thus  when  Ax  is  taken  too  large 
or  there  is  intensity  in  the  spectrum  beyond  n  (2Ax)~  one  will  alias  the 

spectral  information  beyond  »•  back  into  the  region  <v  .  To  iilustraf". 

;  C  c  .1 
suppose  the  true  cutofl  wavenumber  of  a  spectrum  is  125  cm  (for  which 

Ax  -  1  / Z 5 0  cm  4 Op)  but  Ax  i-s  chosen  to  be  5 Op  (corresponding  to  a  v  of 

-1  c_  ! 

100  cm  ).  Under  these  conditions  the  spectral  data  between  0  -100  cm 

will  be  displayed  as  they  should  be,  but  in  addition  the  spectral  data  between 

100 -*125  cm  *  will  be  folded  over  on  top  of  the  data  between  100  cm  and 
1 

75  cm  .  This  is  somewhat  similar  to  the  overlapping  of  orders  in  a  grating 
spectrometer  and  is  usually  undesirable, 

There  is,  however,  one  case  where  the  multi-order  scanning  function 
is  an  advantage.  Suppose  there  is  no  intensity  in  the  regions  corresponding 
to  the  1st,  2nd,  -.  .  .  ,  (n  -  1  )th  order,  but  the t  e  in  intensity  in  the  nth  and 
none  beyond.  Then  Ax  can  he  chosen  as  [2(n  4  1)*^  -  2nv^  j"1  (2v  instead 

of  Ax  =  (2(n  +  l)v  j'1,  and  for  the  same  resolution  fewer  points  need  be 
measured.  This  haw  been  exploited  by  the  Cornua  (Refs.  6,7)  in  their  near- 
1ft  sthtdie  a. 

Making  Ax  smaller  than  need  be  is  equivalent  to  making  v  larger  than 
need  be,  and  obtaining  the  same  resolution  requires  more  points  to  be 
measured  in  theory  there  is  no  disadvantage,  but  in  practice  it  amounts  to 
A  waste  of  measuring  time  and  computer  time  This  is  expanded  on  in  the 
discussion  of  noise  in  interferometry 


1  5 


There  is  nsj  scanning  id  the  formal  sense. 


in  this  manner 


but  it  its  eotr-  cnierrt  To  discuss  rf 
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D. 


CHANNEL  SPECTRA 


One  of  the  most  annoying  phenomenon  encountered  in  FIR  interferometry 
is  channel  spectra.  In  this  frequency  region  (10-125  cm  1 ),  surfaces  not 
plane  parallel  tinder  the  usual  optical  criteria  (a  fraction  to  a  few  wavelengths 
of  visible  light)  may  be  plane  parallel  in  the  FIR  and  will  produce  channel 
spectra.  The  effect  is  quite  easy  to  see  in  IR  interferometry,  and  signature- 
type  features  will  be  easily  detected  in  interferograms  at  x  =  2nd  (Ref.  72), 
where  n  is  the  index  of  refraction  of  the  plane  parallel  slab  of  material,  and 
d  is  its  thickness.  Windows  can  be  suitably  wedged  to  prevent  this  phenomenon, 
and  all  the  windows  have  been  so  treated  in  our  interferometer.  However,  the 
thin  pieces  of  black  polyethylene  used  as  filters  cannot  be  wedged  and  are 
believed  to  contribute  somewhat  to  the  previously  discussed  residual  unbalance 
at  the  origin.  The  channel  signature  will  come  close  to  x  =  0  since  d  ~  lOOp 
and  n  »  1.  5;  of  course,  there  is  some  cancellation  because  of  double  beaming. 

More  troublesome  is  the  fused  quartz  envelope  of  the  high-pressure  Hg 
arc  used  as  the  source.  It  seems  to  have  local  areas  of  plane  parallelism 
and,  since  the  reference  beam  and  sample  beam  do  not  originate  from  exactly 
the  same  region  of  the  plasma  and  envelope,  the  effect  is  iar  from  cancelled 
by  double  beaming.  We  have  tried  dimpling  the  lamp  as  suggested  by 
Richards  (Ref.  18),  but  this  led  to  great  difficulties  in  making  an  equivalent 
match  between  the  two  beams.  We  have  therefore  left  our  lamps  undimpled. 

The  result  is  a  cosine  variation  of  the  source  spectrum  that  is  superimposed 
on  the  sample  spectrum  with  a  period  of  ~1.  3  cm  and  whose  amplitude 
dec rea ss s  at  higher  wavenumbers.  Fortunately  these  channels  are  quite  wide, 
and  their  maxima  and  minima  are  easily  determined;  thus,  they  are  not 
confused  with  the  sharp  lines  observed  in  the  high  resolution  spectra.  How¬ 
ever,  as  the  width  of  the  lines  approaches  the  period  of  variation  the  error 
in  the  frequency  determinations  is  greater.  It  is  estimated  that  the  amplitude 
of  oscillation  (in  intensity)  of  this  channel  spectrum  is  four  to  eight  times  the 
noise  level  in  the  low  wavenumber  region  of  the  spectra. 
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As  mentioned,  the  channels  are  much  stronger  in  the  low  wavenumber 
region  and  are  almost  nonexistent  above  55  cm  ,  From  tbie  we  infer  that 
most  of  the  energy  in  the  low  wavenumber  end  of  the  spectrum  (<  55  cm"*) 

comes  from  inside  the  lamp  (i.  e.  ,  the  plasma)  and  most  of  the  energy  above 

- 1  16 
55  cm  comes  from  the  hot  quartz  envelope. 

E.  LONG  TERM  SOURCE  FLUCTUATIONS  -  STEP  FUNCTIONS 

The  medium-high  pressure  Hg  arc  lamp  has  been  the  most 
frequently  used  source  for  the  far  infrared,  Although  no  comprehensive 
comparison  of  possible  FIR  sources  has  been  made,  it  is  rather  firmly 
established  that  the  medium-high  ores  sure  Hg  arc  is  a  brighter  FIR  source 
than  a  glowbar  (Ref.  83),  a  >.e  discharge  (Ref.  84),  or  a  high  pressure  Na 
lamp  (Ref.  85). 

It  has  been  our  experience  that  £he  stability  of  Hg  arcs  such  as  the 
G.  E.  AH3  and  AII4  is  *he  limiting  factor  (i.  e.  ,  largest  contributor  to  the 
noise)  in  cur  observed  3/N  (if  the  channel  spectra  are  not  considered  noise). 

At  first,  it  would  ae  m  *hat  monitoring  the  source  and  using  this  to  operate 
a  feedback  loop  so  t  s  io  compensate  for  these  fluctuations  would  eliminate 
this  difficulty.  Aside  irern  the  negative  resistance  characteristics  of  these 
lamps  such  a  schema  is  not  practical  since  the  Hg  arc  is  in  reality  two 
sources,  the  plasma,  and  the  hot  envelope  of  the  lamp.  Most,  if  not  all,  of 
the  fluctuations  are  in  the  plasma.  Their  onset  is  quite  rapid  and  may  have 
quite  a  long  duration  (~15  minutes).  We  have  seen,  quite  often,  a  step  function 
type  discontinuity  in  the  interferogram.  Conversely,  the  FIR  energy  output 
from  the  envelope  should  be  quite  steady,  and,  even  if  it  did  change,  such  a 
change  would  be  expected  to  be  slow  and  quite  small  (if  power  input  is 
regulated  to  some  degree).  The  arc  fluctuations  are  caused  by  two  inter¬ 
related  effects:  changes  in  the  voltage -current  characteristics  and  actual 
arc  wandering  over  the  electrodes.  The  phenomenon  is  very  noticeable  in 

°A t  room  temperature,  fuzed  quartz  does  show  a  cutoff  in  the  50-60  cm 
region:  however,  it  is  impossible  to  say  what  the  spectral  transmission 
properties  are  at  the  temperatures  the  envelope  reaches  when  the  lamp  is  on. 
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double-beam  operation  (~4%  of  F'(x)  average).  In  this  respect  we  have  found 
that  the  AH3  lamps  appear  somewhat  more  stable  than  the  AH4  lamps.  We 
surmise  that  this  is  due  to  the  capillary  restrictions  in  the  AH3  lamps  not 
present  in  the  AH4  lamps,  and  thus  the  actual  physical  wandering  of  the  arc 
in  the  AH 3  lamps  is  less  than  in  the  AH4  lamps. 

We  now  consider  what  effect  a  step  function  in  the  interferogram  has 
on  the  spectrum  obtained  by  Fourier  cosine  transformation.  Consider  a  step 
function  of  amplitude  A  starting  at  x^  and  ending  at  x^.  The  effect  of  the 
"box"  in  the  spectrum  is  ascertained  by  Fourier  cosine  transforming  the  "box" 
from  optical  path  difference  space  (x)  into  wavenumber  space  (v).  Thus 


uu  — 

=  J' A(xj,  x^cos  2rxvdx  =  h.  J 


cos  2rtxi'dx 


which  can  be  written 


!( v)box  =  tA(x2  “  xj)l  '  tc°8  nv(x2  +  xi)l  ‘  »8inc  Mx2  *  Xj  )1 


The  first  term,  the  amplitude,  is  proportional  to  a,  which  is  the 
amplitude  of  the  step  function  in  optical  path  Difference  space  and  to  (x^  -  Xj), 
the  duration  of  the  step  function  in  optical  path  difference  space.  The  second 
term  gives  the  rate  of  the  variations  in  wavenumber  space.  If  (x^  +  Xj)  is 
large,  the  fluctuations  will  be  rapid  and  will  appear  to  th?  observer  as  high* 
frequency  noise.  The  last  term  is  the  envelope  (ir.  general)  of  this  function, 
and  it  is  seen  that  tne  largest  fluctuations  cf  t_cu.  ir.  the  low  wave- 

number  region. 

It  ie  possible  to  "edit"  the  interferogram  by  e«timatlnr  A  from  the 
record,  then  subtracting  it  from  the  interferogram  function  from  points  j 
to  X2*  This  will  eliminate  the  gross  effects  of  the  extraneous  box,  but  in 
no  way  compensates  for  the  fact  that  there  has  been  a  real  spectral  charge 
of  intensity  (plasma  change)  or  an  apparent  change  (arc  winner). 
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We  are  testing  and  evaluating  some  "short  arc"  dc  operated  Hg  lamps 
in  the  hope  that  they  will  provide  a  solution  to  this  problem. 

F.  STOP-AND-GO  SCANNING  AND  TIME  CONSTANTS 

As  described,  digital  computation  is  the  most  flexible  and  accurate 
means  of  performing  the  Fourier  transforms  required  in  interferometric 
spectroscopy.  This  requires  that  the  continuous  interferogram  be  sampled, 
meeting  particular  requirements  of  frequency  and  total  number  of  samples 
as  discussed  in  this  report.  The  simplest  scheme  for  taking  those  data  is  an 
on-the-fly  technique  in  which  the  grating  moves  continuously.  Data  are  taken 
during  a  time  interval  sufficiently  short  that  the  grating  has  not  moved 
appreciably.  If  a  high  S/N  were  available  in  the  interferometer  output,  this 
scheme  would  allow  rapid  data  acquisition.  However,  because  of  the  low  S/N 
available  in  the  signals  being  detected,  obtaining  adequate  S/N  in  data  taken 
by  this  method  would  require  extremely  slow  scanning  rates  leading  to  long 
total  run  times  with  much  of  the  time  being  spent  unproductively  in  going 
between  points.  This  problem  can  be  eliminated  of  course  by  moving  rapidly 
between  positions  where  data  are  desired,  then  stopping  grating  motion  during 
the  data-taking  process.  This  stop-and-go  technique  is  the  one  used  in  the 
Aerospace  instrument. 

Proper  operation  of  such  a  system,  however,  requires  a  proper  regard 
for  the  time  constants  of  the  detecting  system.  There  are  three  principal 
time  constants  in  the  detecting  system  in  our  instrument:  (1)  The  time 
constant,  t,  of  the  synchronous  detector  in  the  19-cycle  amplifier;  (2)  the 
gate  time,  T,  used  by  the  preset  counter;  and  (3)  the  period  AT  *  1/3  seconds 
required  for  the  grating  to  move  to  a  new  position.  The  first  controls  the 
highest  frequency  noise  that  will  appear  in  the  output;  the  second  controls  the 
length  of  time  ("integration  time")  over  which  the  data  are  averaged;  the  third 
is  only  indirectly  under  the  control  of  the  experimenter. 

These  filtering  parameters  are  adjusted  for  any  particular  run  to 
optimise  the  conflicting  requirements  of  a  low  noise  level  and  a  reasonable 
total  run  time.  Failing  to  recognise  that  the  adjustments  are  interdependent, 
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however,  may  lead  to  systematic  errors  in  intensity  in  the  final  spectrum. 

If  the  ratio  of  the  integration  time  to  the  detector  time  constant  is  too  low  and 
if  the  signal  level  is  quite  different  from  the  preceding  point,  the  signal 
averaged  by  the  counter  will  be  changing  over  an  appreciable  portion  of  the 
period  of  averaging.  This  adds  to  the  interferogram  another  function 
proportional  to  the  first  derivative  of  the  interferogram. 

The  magnitude  of  this  contribution  can  be  estimated  with  the  aid  of 
Fig.  59,  a  schematic  representation  of  the  signal  changes  during  one  data, 
acquisition  cycle.  At  t  ,  data  acquisition  at  the  preceding  location  is 
completed;  between  t  and  the  grating  moves  to  a  new  location  and  stops. 
For  convenience  the  signal  (solid  line)  is  assumed  to  vary  linearly  during  this 
time.  At  tj,  data  taking  at  the  new  point  is  initiated  and  continued  throughout 
a  period  T  until  time  t^.  The  output  signal  from  the  19 -cycle  detector,  with 
a  poor  choice  of  time  constants,  is  shown  as  the  dashed  line.  The  fractional 
error  in  the  data  taken  then  will  be  the  ratio  of  the  hatched  area  to  the 
rectangle  of  area  TE^,.  This  fractional  error  is: 

><  [m  -  »-T/T)] 

The  first  set  of  parentheses  represents  the  amount  of  change  in  the  signal  and  thus 
the  largest  possible  error.  The  first  square  brackets  represent  the  decrease 
in  this  maximum  error  because  the  output  signal  has  some  time,  AT,  in  which 
to  start  to  follow  the  input  signal  change;  the  second  brackets  introduce  the 
error  caused  by  the  input  to  the  counter  changing  during  the  data  taking. 

Since  the  noise  is  expected  to  go  down  with  tha  squsre  root  of  the 
Integration  time,  increasing  T  above  about  10  seconds  increases  the  overall 
run  time  more  than  is  worthwhile  for  the  improvement  in  noise.  Since  a  pen 
recorder  runs  off  the  output  of  the  detector,  an  offhand  choice  for  t  would  be 
a  little  longer  than  the  i -second  response  time  of  the  recorder,  possibly  3 
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Fig.  59.  Schematic  representation  of  signals  into  and  out  of  the 
19-cycle  detector  during  one  data  acquisition  cycle. 
During  time  AT,  grating  is  moving  to  new  position. 
During  time  T,  the  "integration  time",  data  are  taken. 
Hatched  area  is  the  amount  by  which  the  final  inte¬ 
grated  information  will  be  in  error  due  to  the  constant 
of  the  detector. 
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seconds.  With  these  values  aim'-.-,?  -.o  correction  takes  place  during  the 
moving  time  (  [  ]  =  0.  95),  whereas,  the  error  during  the  integration  process 
is  still  about  30%.  The  solution  is  to  decrease  the  detector  time  constant  as 
much  as  possible  and  to  independently  increase  the  recorder  time  constant  if 
desired.  If  we  reduce  t  to  0.  5  second,  25%  of  the  change  can  occur  while 
the  grating  is  moving.  The  error  in  the  integration  is  about  b:nt so  the  total 
error  is  now  only  a  little  above  3%,  a  ten-fold  reduction.  The  effects  oi  this 
improvement  in  accuracy  in  intensity  measurements  on  the  transformed 
spectrum  may  be  seen  from  Fig.  60,  which  shows  two  spectra  of  a  mercury 
vapor  source  taken  under  identical  circumstances,  except  for  the  time 
constant  of  the  i9-cycle  detector.  From  this  figure,  it  is  apparent  that 
proper  choice  of  time  constants  is  imperative  for  a  system  having  several 
cascaded  systems,  each  with  a  comparable  individurl  time  constant. 


Fig.  60.  Effect  on  the  spectrum  of  a  large  detector  time  constant. 

Both  curves  are  spectra  of  a  medium-pressure  mercury 
vapor  source,  and  except  for  the  detector  time  constant, 
were  taken  under  identical  operating  conditions.  The  in¬ 
tegration  time  was  10  sec/point.  The  solid  line  corres¬ 
ponds  to  a  detector  time  constant  r  of  0. 50  sec;  the  dashed 
line  is  the  iptcv/um  obtained  with  a  r  of  3  sec. 
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G. 


BACKGROUND  RADIATION 


In  any  spectrometry,  and  particularly  in  the  FIR,  background  radiation 
competing  with  that  from  the  sample' can  introduce  serious  errors  in  measure¬ 
ments.  It  is  important  that  this  source  of  error  be  understood  and  that  all 
possible  means  to  minimize  it  be  used.  In  order  to  assess  this  error,  con¬ 
sider  first  the  PlnnrV  r3  diction  law  as  applied  to  these  long  wavelengths.  It 
can  be  written 


cA 

e<x>=tot 

e  L  -  1 


where  and  are  constants  and  X.  and  T  are  the  wavelength  and  the 
absolute  temperature  of  source.  For  KT  »  c by  expanding  the  exponential 
and  keeping  only  the  leading  term, we  obtain  the  Rayleigh -Jeans  approximation 


E(M 


where  we  see  that  E(X.)  is  simply  proportional  to  T,  not  rapidly  varying  as 
it  is  for  small  X.  In  a  room  temperature  (300*K)  environment  with  a  source 
at  3000*K  (very  high  for  available  sources),  the  radiance  (W/cm^  ster)  is 
only  10  times  greater  for  the  source.  The  source  in  our  case  fills  a  solid 
angle,  as  seen  by  the  detector  through  the  interferometer,  of  about  1/50, 
whereas,  the  background  fills  l*,  which  is  larger  by  about  300.  Hence, the 
integrated  energy  from  die  background  exceeds  that  from  the  source  by  a 
factor  larger  than  30  even  before  absorbing  samples  are  inserted  in  the  beam. 

Fortunately,  two  factors  act  to  greatly  reduce  die  effect  of  the  back¬ 
ground.  First,  moat  of  the  background  radiation  does  not  pass  through  the 
interferometer  and  it  is  not  modulated  by  the  interferometer  displacement; 
hence,  it  introduces  a  constant  aero  shift  in  the  interferograzn  that  is  ignored 
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in  the  computer  transformation.  Second,  an  even  larger  fraction  of  the 
background  is  not  chopped  and  the  detector  response  it  produces  does  not 
contain  the  19-cycle  modulation  frequency  to  which  the  amplifier  is  tuned. 
Radiation  must  pass  through  bcth  the  interferometer  and  the  chopper  to  be 
transformed  into  an  output  spectrum.  In  any  case,  however,  it  adds  to  the 
detector  noise.  Also  the  chopped -nonmodulated  component  uses  up  some  of 
the  available  dynamic  range  ot  .he  amplifier. 

In  order  to  minimize  the  effects  of  stray  chopped  radiation.it  is 
necessary  to  baffle  the  system  carefully.  Baffling  in  this  spectral  region 
is  more  difficult  than  it  is  for  the  shorter  wavelength  region  because  most 
of  the  usually  black  materials  are  not  black.  Even  black  flocking  reflects 
about  20%.  Also,  ordinarily  opaque  materials  such  as  heavy  cardboard 
transmit  about  10%.  Sheet  metal  appears  to  be  about  the  only  truly  opaque 
baffling  material,  and  we  have  not  found  any  really  good  nonreflecting  surfaces 
We  have  been  able  to  reduce  the  chopped -nonmodulated  component  to  less  than 
10%  when  highly  transmitting  samples  are  being  measured,  but  the  percentage 
increases  drastically  for  absorbing  materials. 

If  an  absorbing  sample  is  in  the  sample  holder,  the  chopped  and 
modulated  radiation  received  by  die  detector  when  the  sample  is  being  viewed 
consists  of: 

1.  Desired  radiation  from  the  source  transmitted  by  the  sample 

2.  Radiation  emitted  by  the  sample 

3.  Radiation  from  the  surroundings  that  is  reflected  by  the  sample 
into  tihe  beam 

4.  Radiation,  similar  to  items  2  and  3,  from  the  sample  holder  and 
other  surrounding  structures. 

There  are  similar  components  in  the  reference  beam  which,  of  course, 
subtract  from  the  signal  due  to  the  sample  beams. 

It  is  easy  to  estimate  the  effect  of  sample  emission.  Consider  a  source 
of  effective  temperature  radiating  energy  through  a  sample  ar  ambient 
temperature  T  .  The  sample  has  an  absoTbtivity,  a,  that  equals  its 
emissivity,  e. 


If  the  sample  were  in  equilibrium  with  its  surroundings  at  T  ,  it  would 
absorb  a  fraction,  a,  of  the  radiation  of  wavelength  X.  passing  through  it 
within  a  given  solid  angle,  but  it  would  also  emit  an  equal  amount  into  the 
same  solid  angle.  In  our  case,  the  source  radiation  enteri  ig  the  sample  is 
T  /T  times  as  large  because  of  the  higher  temperature  of  the  source.  Hence, 
the  radiation  observed  is  proportional  to  T  (1  -  a)  +  aT  =  T  -  a(T  -  T  ). 

8  &  0  5  a 

An  apparent  absorption,  a'  -  a[  i  -  (T  /T  )] ,  is  observed  that  can  never  be 

2L  0 

larger  than  1  t  /t  With  the  commonly  used  quartz  Hg  arc  source,  the 

*  •  _l 

hot  quartz  is  the  radiating  material  for  v  near  100  cm  If  T  is  then 
a^OO’K,  the  error  in  a  is  *30%.  The  observed  transmittance  will  then  never 
be  less  than  30%. 

The  reflection  by  the  sample  allows  the  detector  to  see  other  emitting 
surfaces  near  the  sample.  These  surfaces  can  have  any  spectral  radiation 
profiles  and  hence  the-distortion  introduced  into  the  measured  spectrum 
cannot  be  predicted.  It  will  not  be  a  simple  multiplier  as  in  the  case  of 
absorption  but  may  even  introduce  extraneous  features  into  the  spectrum. 

An  interesting  special  case  arises  when  the  sample  acts  as  a  partially 
reflecting  mirror  perpendicular  to  the  beam.  The  detector  then  "sees  itself,  " 
and  since  the  detector  is  cooled,  the  magnitude  of  the  error  is  reduced.  This 
effect  has  been  observed  with  metal  meehea  run  single~beam,  the  reference 
beam  being  blocked  off  by  a  blackened  metal  plate.  The  metal  meshes  are 
highly  reflective  and  the  detector  sees  itself  in  the  sample  beam,  but  it  sees 
the  300*K  surroundings  in  the  reference  beam.  The  measured  "transmittance" 
is  negative  as  a  result. 

It  should  be  pointed  out  dust  under  these  circumstances  the  detector  does 
not  see  just  itself  but  half  of  itself  and  half  of  the  surroundings.  The  inter¬ 
ferometer  puts  half  of  incident  energy  into  the  central  maximum  and  half  into 
the  higher  orders.  Looking  backward  through  it,  we  see  an  extended  field  of 
view  with  half  transmittance  from  the  position  of  the  central  maximum  and 
hall  from  the  positions  of  the  higher  orders.  Since  the  detector  is  very  cold, 
the  effective  temperature  that  the  detector  sees  under  these  conditions 
approaches  T^/2  as  extraneous  scattered  radiation  is  reduced. 
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H. 


DRIFT 


During  the  data  taking,  changes  in  signal  intensity  with  time  constants 
of  the  order  of  hours  are  observed.  In  a  particular  run  there  are  varying  contri¬ 
butions  of  two  types  to  these  changes:  (1)  an  effect  correlated  with  the  motion 
of  the  lamellar  grating  as  reflected  in  motion  of  the  interferometer  entrance 
aperture  image  in  the  exit  aperture;  {Z)  nonreproducible  drifts  of  detector, 
source,  and  critical  amplifiers.  Independent  of  the  cause,  however,  is  the 
eiiect  of  these  slow  changes  on  the  output  spectrum.  They  lead  to  an  increase 
in  low  wavenumber  noise  in  the  spectrum  and  may  have  some  effect  on  line 
sh«pes  in  higher  wavenumber  regions  of  the  spectrum. 

To  illustrate  thisywc  discuss  results  of  a  detailed  study  of  a  water 
vapor  interferogram  which  showed  a  particularly  simple  long-term  change, 

.1  murr-or-less  linear  decrease  with  increasing  path  difference.  This  ia 
shown  in  Fig.  61,  which  s  a  plot  of  averages  of  10  adjacent  points  taken  at 


OPTICAL  PATH  DIFFERENCE  (cm) 

Fig  61.  Averages  over  adjacent  10-point  regions  of  a  water  vapor 
interferogram,  showing  a  case  of  nearly  linear  drift. 
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100-point  intervals  in  the  inlerferogram.  A  first  approximation  to  the 

apparcnl  Lrm  of  the  iuterferogram,  F*  (x)  -  r'  (x  )  !  F  (x)  is 

a  am  a 

F' (x)  =  F'  (x)  (  1  -  ax) 
a  w 

where  F^'x)  =  F'Jxm)  +  FJx)  is  the  interferogram  due  to  the  water  alone, 

and  the  quantity  in  brackets  introduces  the  linear  fall-ofi.  The  spectrum 

then  is  obtained  by  Fourier  transforming  F  (x)  and  is 

a 


x 

m 

r 

E[v)  ~  i  F  (x)cos  Zxu'xdx 

/  a 


X 

m 

=  al  (x  )  /  (x  -  x)cos  2irvxdx 

w  m  J  m 

o 


(43) 


The  first  term  in  Eq.  (43)  makes  its  main  contribution  at  v  =  0  and  falls  off 
proportional  to  v  although  it  is  large  in  magnitude  at  v  ~  0  in  the  region  of 
the  spectrum  where  there  is  sufficient  energy  to  make  observations, 
v  k,  20  cm  ,  this  term  will  make  little  contribution. 

To  gain  an  idea  of  the  effect  of  the  additional  (1  -  ax)  factor  in  the 
integrand  of  the  second  term  of  Eq.  (43)tif  is  convenient  to  assume  we  are 
discussing  lines  with  a  lx>rentzian  shape,  E(v)  =  [(v  -  v  p  4  (Av)  }~  ,  of  half 
maximum  width  of  Av  centered  at  frequencies  v  .  It  has  bee:i  shown  thwt  the 
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L n  FREQUENCY  r(cnf') 


Fig.  62.  Difference  between  the  spectrum  near  v  =  0 
as  obtained  from  the  unr  jrrected  interfero- 
gram  .nd  that  obtained  from  the  same  inter- 
ferogram  alter  dividing  each  point  by(l  -  ax). 
The  straight  line  has  slope  -Z,  as  suggested 
by  the  v'^factor  in  the  first  term  of  Eq.  (43), 
The  x-marks  along  the  v  axis  are  zeros  of 


interferogram  expected  from  this  line  is  an  exponentially  decaying  cosine, 
cos  2ttvqx  exp(-2irAvx).  Then  the  apparent  spectrum  given  by  Ej(«/)  for  this 
line  will  be 


ft) 

Ej(v)  =  J cos{2Tri/ox)cos(2'tr\Jx)e”^1Tx^v(  1  -  ax)dx 
o 

If  a  is  not  too  large,  (1  -  ax)  may  be  thought  of  as  the  first  terms  in  the  expan¬ 
sion  of  exp(-ax).  The  form  of  the  function  to  be  transformed  then  is  the  Lorentss 
line  interferogram.  Thus  the  line  shape  must  also  be  Lorentrian  with  an 
apparent  line  width  Av&  =  Av  +  a/2ir  somewhat  wider  than  the  true  line  width. 

The  applicability  of  this  approach 


Fig.  63.  Lines  obtained  from  the 
corrected  interferogram  (dashed 
line)  and  from  the  uncorrected 
interferogram  (solid  line).  Note 
the  decrease  in  width  at  half¬ 
maximum  intensity  obtained  at  the 
expense  of  some  increase  in  the 
apparent  noise  in  the  base  line. 


is  demonstrated  by  dividing  each  point 
in  the  apparent  interferogram  by  (1-ax), 
thus  obtaining  the  actual  interferogram. 
The  difference  between  the  transforms 
ox  the  corrected  and  uncorrected  inter- 
ferograms  near  u  =  0  appears  as  a 
loglog  plot  in  Fig.  62.  The  decaying 
sine  curve  is  obvious.  The  straight 
line  has  the  expected  die-off  slope  of 
-2  and  is  seen  to  agree  well.  Figure  63 
shows  one  of  the  lines  in  the  water  vapor 
spectrum.  The  solid  line  is  the  trans¬ 
form  of  the  uncorrected  interferogram, 
whereas  the  dashed  line  obtained  from 
the  corrected  interferogram  shows  the 
narrowing  effect  of  eliminating  the 
drift  from  the  interferogram. 
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The  particular  analytic  expressions  developed  here  are  applicable 
only  to  interferograms  of  separated  lines  of  Lorentzian  shape  which  show 
a  drift  linear  with  increasing  position.  The  effect  of  other  long  period  drifts, 
however,  would  be  expected  to  be  similar,  introducing  increased  noise  at  low 
wavenumbers  in  the  spectrum  and  making  changes  in  the  line  shapes  at  higher 
wavenumbers. 

I.  NOISE 

Let  the  interferogram  function  FJx^)  be  measured  at  points  xq,  x^, 
x~,  •  x  ,  ■  •  • , xM{  =  L).  We  assume  the  separations  Ax(=x.  -  x.  , )  are 

Cm  n  IN  |  J  J  “  l 

all  equal  and  chosen  such  that  Ax  =  (2v^)  in  accordance  with  the  sampling 

theory  condition  (see  Section  II).  In  each  F(xfl)  there  will  be  a  contribution 

that  is  the  noise  and  is  called  A_,  A.,  A_,  •  •  •  ,  A  ,  •  •  •  ,  A„T.  The  A  are 

u  l  c  n  In  n 

assumed  to  be  independent  of  x  and  F(xn)  and  furthermore  their  mean  value 
is  zero.  We  then  take  as  a  measure  of  the  noise  in  the  interferogram  its 
standard  deviation  Nj  given  by 


The  spectrum  is  obtained  from  the  interferogram  function  by  Fourier  cosine 
transformation  and  thus  the  component  of  noise  in  the  spectrum,  Bm,  at  vm 
is  given  by  ^ 


N 


B  =  Ax 


m 


n=0 


cos  2irx  v 
n  m 


For  the  term  n  =  0,  we  should  use  A  /2  rather  than  AQ,  but  as  long  as 
N  is  reasonably  large,  the  error  is  small. 
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We  shall  calculate  the  B  only  at  those  v  that  are  independent. 

m  7  j  j  m 

The  independent  v  are  located  At>  =  (T<)~  cm"  apart,  and  taking  the  first 
-i  m  -i 

v  at  0  cm  we  have  v  =  mAv  =  m{  L)  ;  furthermore,  x  =  nAx,  and 

m  j  in  n 

Ax  =  (2vc)  .  Making  these  substitutions  we  obtain 


B'n  '  2vc  E  A”  C0SGcL1 


The  mean  value  of  the  B  's  is 

m 


iVi 


M  N 


^m  2(M  +  l)v 


m=u  n  = 


1  cosfe\ 

n  \voh ) 


Interchanging  the  sum  over  m  and  that  over  n,  and  multiplying  by 
Am(  =  l),  and  then  approximating  the  sum  over  m  by  its  limiting  integral  form, 
we  obtain 


N  vcL 


m  ~  2(M 


and  making  the  change  in  variable  y  =  (mm/v  L)  the  integral  reduces  to  the 

fit  c  _ _ 

form  J  cos  nydy,  which  is  aero.  Therefore  Bm  =  0.  We  proceed  to  calculate 
the  standard  deviation  of  the  noise  in  the  spectrum.  Squaring  Eq.  (44)  and 
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summing  over  m,  we  obtain 


M  N 


m=0 


Again,  by  interchanging  the  sum  over  m  with  the  other  sums  and 

approximating  this  sum  by  the  limiting  integral,  then  dividing  by  M(=N/2), 

2 

we  obtain  the  standard  deviation  squared  of  the  noise  in  the  spectrum  N  . 

8 


(45) 


which  is  equivalent  to  the  expression  given  by  Connes  (Ref.  6)  and  Bell(Ref.  8). 

We  have  carried  out  the  following  numerical  exercise  to  confirm  Eq.  (45). 

One  thousand  numbers  were  chosen  randomly  with  a  normal  (Gaussian) 

distribution  such  that  the  standard  deviation  was  1000  units.  These  one 

thousand  numbers  were  then  used  on  a  noise  interferogram  and  were  Fourier 

cosine  transformed  with  various  values  of  L  and  Ax.  Ihe  standard  deviation 

in  the  noise  spectra  were  then  computed,  and  Hie  results  are  summarised  in 

Table  1 3.  The  ranges  of  L  and  Ax  were  Chosen  to  be  representative  for  the 

FIR  region.  Equation  (45)  predicts  that,  if  Ax  is  decreased  (or  v  multiplied) 

c 


Table  13.  Calculated  (from  Eq.  45)  and  Observed  Standard  Deviations 


L,  cm 

Ax,  p 

NI 

Ng  (calc) 

Ng(obs) 

4.  00 

40 

100C 

80 

86 

3.  60 

40 

1000 

H5 

83 

3.  20 

40 

1000 

80 

79 

2.  *80 

40 

1000 

75 

74 

2.  40 

40 

1000 

69 

69 

2.  00 

40 

1000 

63 

64 

8.  00 

80 

1000 

179 

172 

4.  00 

80 

1000 

127 

126 

6.  00 

60 

1000 

134 

131 

3.  00 

60 

1000 

95 

97 

2,  00 

20 

1000 

45 

44 

1.  00 

20 

1000 

32 

31 

i  .  00 

10 

1000 

22 

21 

0.  50 

10 

1000 

16 

15 

by  a  factor  R,  the  noise  in  ihe  spectrum  will  be  reduced  by  a  factor  R*  ^ . 

The  number  of  points  will  of  course  be  increased  by  R  (for  constant  L).  Thus, 

if  the  same  amount  of  time  is  available  for  die  measurements,  NL  will  be 

1  /2  1 
Increased  by  R  resulting  in  no  change  for  N  .  In  practice,  however,  it  is 

still  advisable  to  select  Ax  as  large  as  possible  equal  to  (2v  )  to  minimise  the 

c 

"dead  time"  (the  time  it  takes  the  interferometer  to  go  from  point  to  point, 
during  which  no  data  are  taken),  since  in  general  the  smaller  Ax  becomes 
the  more  dead  time  there  will  he  per  unit  length  of  optical  path  difference. 


Making  Ax  as  large  a*  possible  also  minimizes  the  total  number  of  points  in 
the  interferogram  which  can  be  of  some  consequence  when  computer 
economics  are  considered. 

Connes  (Ref.  6}  has  considered  the  case  of  an  interferogram  that  is  a 

cosine  function.  This  corresponds  to  a  line  in  the  spectrum  having  zero  width. 

The  signal  (maximum  intensity  of  the  line)  is  proportional  to  L,  and  therefore 
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the  S/N  is  proportional  to  L  .  For  such  a  case  the  interferogram  should 

be  carried  out  as  far  in  optical  path  difference  as  possible  since  an  increase 

in  L  always  increases  the  S/N  in  the  spectrum.  This  is  a  limiting  case  and 

is  appropriate  when  the  resolution  attained  in  the  interferogram  (governed 

by  L)  is  much  less  than  the  width  of  the  line.  When  the  resolution  becomes 

comparable  to  the  line  width,  then  the  dependence  of  S/N  in  the  spectrum  will 
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be  different  than  the  L  dependence. 

Consider  a  spectrum  I(v)  made  up  of  a  single  line  having  a  Lorentz  shape. 
This  line  shape  was  chosen  since  it  represents  a  good  description  of  weakly 
absorbing  pressure  broadened  lines,  and  is  mathematically  tractable  for  the 
Fourier  cosine  transformation  we  wish  to  carry  out.  And  now  calling  the 
signal  in  the  spectrum,  Sg  =  Il(vq)  (see  Eq.  10), 

Ss  •Uvjl  -  e"i2*tLf)  (46) 

and  therefore  the  S/N  in  the  spectrum  is 


(47) 


where  Sj  =  F(0),  the  maximum  signal  in  the  interferogram. 

The  S/N  in  die  spectrum  has  a  maximum  when  L<  *  0.  200.  For 
L«  <  0.  2,  die  S/N  curve  rises  ^  as  is  to  be  expected  since  L  <  1  /«;  for 
L«  ~  0.  2  the  curve  flattens  out,  then  falls  gently  toward  the  asymptotic 
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dependence  of  L  .  If  x  is  set  equal  to  (5«)  in  the  interferogram  function 
(see  Section  II),  then  the  amplitude  of  the  interferogram  function  at  this  value 
of  x  will  be  down  to  -28.  5%  of  its  maximum  value  at  x  =  0.  Thus,  for 
maximum  S/N  in  the  spectrum,  the  interferogram  should  be  terminated 
when  it  has  been  damped  down  to  28.  5%  of  its  maximum  value. 

In  Fig.  64  is  a  plot  of  I^(v)  for  Lf  =  0.  200,  and  also  I(v).  The  intensity 
of  IL(v)  differs  markedly  from  I(v).  It  is  -28%  lower.  Also,  the  character¬ 
istic  sidelobes,  or  feet,  are  quite  prominent.  For  Lc  =  0.  200,  the  intensity 
of  the  first  sidelobe  (the  strongest)  is  (measured  from  the  first  minima  to  the 
second  maxima)  ~16%  of  the  central  maxima  (measured  from  the  first  maxima 
to  the  first  minima).  In  such  a  spectrum,  if  the  S/N  is  greater  than  -6,  these 
sidelobes  might  be  interpreted  as  real  lines  and,  therefore,  it  would  be 
desirable  to  suppress  them.  The  sidelobe  intensity  depends  on  L,  and  in 
general  the  sidelobe  intensity  will  be  decreased  by  making  L  larger. 

The  quantity  Ip  (discussed  in  Section  II)  represents  the  maximum  intensity 
of  the  first  (and  strongest)  sidelobes;  therefore,  1^(1^) /Ip  measur®8  the 
maximum  intensity  to  sidelobe  intensity.  There  is  little  to  be  gained  in 
carrying  the  interferogram  beyond  that  L  which  yields  I*  (vQ)/I p  much  larger 
than  the  S/N.  The  function  I^(v)/lpVS  2wL*  is  shown  as  curve  D  in  Fig.  65. 
Also  shown  are  three  Sg/Ng  curves  characterized  by  the  parameter 
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-  100  (curve  A)  ;  300  (curve  B)  ;  and  500  (curve  C) 

Consider  for  example  curve  B.  If  maximum  S/N  in  the  spectrum  is  required 
the  interferogram  should  be  cut  off  where  L<  ar  0.  20, as  discussed  previously, 
and  this  would  yield  a  S/N  in  the  spectrum  of  — 1 90/ 1 .  However,  if  we  require 
that  the  sidelobes  be  the  same  size  as  the  noise,  then  we  are  interested  in 
where  curve  B  intersects  curve  D,  which  happens  at  2«Lc  •*  5.  35.  The  S/N 
in  the  spectrum  will  then  be  -130/1.  In  this  case  die  interferogram  will  be 
carried  down  to  e  »  0.  005  of  its  maximum  value.  Then  the  difference 
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I  (v  ) 

(CURVES  A,  B,  AND  C  ),  ( CURVE  D) 


400 


Fig.  45,  Spectrum  S/N  curve#  <A,  B,  and  G),  and  maximum 
isitenaity  tn  aidelobe  inter***  ty  curve  (D). 
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I(v0)  -  Ij^(v0)  will  be  also  within  the  noise.  If  (Sv^/rtt)  >  1  (and  generally 
it  is),  then  the  interferogram  will  be  carried  out  beyond  where  the  signal  in 
the  interferogram  is  lost  in  the  noise  of  the  interferogram. 

Two-beam  interferometry  is  usually  applied  to  rather  broad  band 
absorption  investigations  (this  does  not  of  course  exclude  high  resolution). 
This  will  lead  in  general  to  two  complications  not  considered  in  the  previous 
discussion.  First,  the  quantity  ws  have  been  calling  the  signal  in  the  inter¬ 
ferogram  (which  can  be  thought  of  approximately  as  the  variation  of  the 
interferogram  about  its  average  value)  is  usually  only  a  small  percentage  of 
the  total  signal  detected  by  the  detector.  This  is  die  well-known  dynamic 
range  problem  of  two-beam  interferometry.  It  can  be  suppressed  by  the 
technique  of  double-beam  differencing  (see  Section  HI).  Second,  usually 
more  than  one  line  is  present  in  the  spectrum  being  examined  interfero- 
metrically,  and  the  line  width  parameter*  <  may  not  all  be  equal.  Thus, 
the  interferogram  will  not  have  a  simple  damped  cosine  form,  but  rather 
will  be  the  net  result  of  many  of  these  terms  (one  for  each  line).  In  the  case 
of  diatomic  and  linear  molecules  the  arguments  of  these  cosines  are  related 
nearly  harmonically  and  the  grosser  features  of  the  interferograms  are  its 
well-known  signatures  (Refs.  22.  23), and  the  damped  amplitude  of  these 
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signatures  are  then  a  measure  of  the  (total)  signal  in  the  interferogram. 

At  the  other  extreme  is  the  pure  rotational  interferogram  of  the  water 
molecule.  This  interferogram  is  very  irregular  (Ref.  15);  however,  the 
amplitude  does  decay  (on  the  average)  as  one  goes  to  larger  x  and  an  average 
assessment  of  the  damping  is  possible.  To  a  certain  extent  in  pure  rotational 
spectroscopy  the  line  width  parameter  is  somewhat  in  the  control  of  the 
experimenter  since  it  depends  on  the  pressure  of  the  absorbing  gas.  For 
best  S/N  in  die  spectrum,  die  pressure  in  the  absorption  cell  should  be  fixed 
so  that  die  average  of  the  signal  in  the  interferogram  in  the  vicinity  of 


A  word  of  caution,  however.  In  the  pure  rotational  interferogram  (Ret.  If) 
of  DCi,  the  signatures  die  out  and  then  come  back  at  larger  L  because  die 
spectrum  is  made  up  of  a  series  of  closely  spaced  doublets. 


x  =  L  is  -28 %  of  the  average  of  the  signal  in  the  vicinity  of  x  =  0  (the  region 
about  x  =  0  should  be  excluded  from  the  average).  This  will  provide  (on  the 
average)  the  best  compromise  between  S/N  in  the  spectrum  and  sharp  lines 
and  presumably  the  best  condition  for  wavenumber  measurements.  If  there 
are  overlapping  lines  the  presence  of  the  sidelobes  will,  in  general,  introduce 
an  error  in  the  wavenumber  measurements.  This  is  discussed  later. 

The  line  shape  assumed  is  a  good  approximation  for  weak  absorption. 

For  more  strongly  absorbing  lines  the  Lambert-Beers  law  is  more  appro- 

19 

priate,  and  in  general  the  line  shape  will  be  modified.  If  such  a  line 
(hereafter  called  Lorentz-Lambert-Beers  line  or,  for  short,  LLB  line)  is 
compared  with  a  Lorentz  line  normalized  such  that  both  have  the  same 
maximum  intensity  and  same  effective  half-width  at  half-maximum,  then 
the  L  line  will  be  "squarer"  in  the  sense  it  is  broader  (more  intense)  between 
v  -  t  ,,  and  v  +  «  ,,  and  narrower  (less  intense)  outside  this  range.  We 
have  generated  a  number  of  I  (v)  for  various  values  of  y{vQ)  /  >  1  and  Fourier 
cosine  transformed  these  functions,  the  result  being  the  interferogram 
function  they  would  present.  In  the  region  xc  <  1  these  interferogram  a  were 
damped  out  -exp  -  |2*«  „x|,  but  in  the  region  x(«)ef-£  ~  1  the  cosine  modula¬ 
tion  came  back  with  amplitude  small  but  not  negligible  in  comparison  with 
its  amplitude  near  x  =  0.  To  a  first  order  of  approximation  the  difference 
between  these  two  line  shapes  iooks  like  two  lines  separated  by  «  „  and  thus 
one  would  expect  subsidiary  maxima  and  minima  when  *(0^  =  l,  2,  *  *  •  . 

The  magnitude  of  the  amplitude  of  the  interferogram  at  x { c  )  ^  *  1,2.  3,  • 
is  of  course  determined  by  the  difference  between  the  LLB  line  shape  and  the 
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When  an  interferogram  is  measured  using  the  double-beam  differencing 
mode  of  operation  (Ref.  15)  the  result  (after  Fourier  cosine  transformation) 
is  not  I(v>  the  spectrum,  but  rather  v)  --  [jL(v)  -  Ifv)),  where  I^v)  is  the 
spectrum  that  would  be  obtained  with  no  sample  in  the  cell.  If  we  assume  a 
Lambert-Beers  law  for  absorption,  we  hava  ^j(v)  -  I^v#!  -  v'~Y(v')f),  where 
y(v)  is  the  absorption  coefficient  and  .!  the  absorbing  path.  For  y(v)f  <<  1, 
this  reduces  to  v)  at  ,  For  pressure  broadened  lines  y(  v)  has 

the  form  given  by  Eq  (2  i). 
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Lorentz  line  shape.  As  far  as  the  spectrum  is  concerned  the  contribution 
from  these  parts  of  the  interferogram  will  "square”  out  the  line. 

Wo  have  assumed  no  error  in  the  measurement  of  x,  the  optical  path 
difference.  ^  Such  an  error  can  be  quite  degrading  as  far  as  S/N  in  the 
spectrum  is  concerned.  Consider  for  example  an  interferogram  which  has 
large  amplitude  features,  such  as  signatures.  In  the  vicinity  of  a  signature, 
the  rate  of  ch  :nge  of  intensity  in  the  interferogram  is  very  large  with  respect 
to  a  change  in  x,  and  thus  a  small  error  in  x  could  lead  to  a  large  change 
in  intensity  of  the  interferogram.  This  is  very  undesirable  in  that  the  apparent 
noise  generated  by  such  an  error  could  be  very  much  larger  than  the  true 
ncise.  Furthermore,  these  larger  apparent  noise  excursions  occur  at 
exactly  those  places  in  the  interferogram  that  correspond  to  the  positions 
of  greatest  contribution  to  the  intensity  of  the  lines  in  the  wavenumber 
spectrum.  Thus  the  noise  level  in  the  wavenumber  region  of  the  lines  could 
be  greater  than  in  the  wavenumber  regions  where  there  are  no  lines. 

Finally.,  we  consider  the  question  of  resolution.  Shown  in  Fig.  66  are 
two  identical  linas  with  Lc  =  0.  200  (optimum  S/N  situation)  separated  by  4« . 
The  dip  in  intensity  between  the  lines  is  -l  i%  (somewhat  smaller  than  the 
Rayleigh  criterion;  if  the  separation  were~4.  the  dip  in  intensity  would  be 
~18%,  very  close  to  the  Rayleigh  criterion).  The  lines  are  separated  by  4« 
but  appear  (measurement  of  maximum  intensities)  to  be  separated  by  ~3.  8 c. 
Well  separated  lines  could  be  measured  more  accurately  (depending  of  course 
on  the  S/N).  Apodization  of  the  interferogram  would  suppress  the  sidelobes, 
but  little  or  no  gain  in  measurement  accuracy  would  be  expected  since  the 
lines  would  then  be  broader  and  their  central  intensity  region  would  overlap 
even  more. 


If  the  error  is  constant  fi.  e. ,  mislocation  of  the  origin,  x  =  0)  then  serious 
shifts  and  distortions  of  line  shape  may  result.  If  the  error  is  periodic, 
ghosts  will  be  seen  in  the  spectrum  (Ref.  6).  See  Section  V-A,  B. 


Fig.  66.  Two  identical  Lorentz  lines  separated  by  4c ,  and 
L  c  =  0.  200. 


VI.  COMPUTATIONAL  TECHNIQUES 


A.  INTRODUCTION 

Development  of  efficient  means  to  carry  out  the  calculations  indicated 
in  Eq.  (5)  has  enabled  Fourier  transform  spectroscopy  to  become  a  practical 
means  of  spectroscopy.  Because  the  computational  techniques  are  so 
necessary,  this  section  is  devoted  to  a  discussion  of  the  algorithms  used 
and  tests  that  may  be  applied  to  them,  and  it  closes  with  a  brief  discussion 
of  programming  the  algorithms  for  use  on  a  high  speed  computer,  including 
flow  diagrams.  With  detectors  presently  available,  the  only  meaningful 
energy  information  available  is  relative,  not  absolute.  Because  of  this, 
constant  factors  are  customarily  left  out  of  the  equations  of  Fourier  transform 
spectroscopy.  When  this  is  done,  Eq.  (5)  becomes 


E(v)  =  Ax 


N 

F(0)  +  2£ 


j=l 


F(jAx)  cos  2TrjvAx 


(48) 


Notice  the  finite  sampling  of  the  interferogram  forces  the  computed 
spectrum  E(v)  to  be  periodic  with  period  1  /Ax.  Furthermore,  since 
coa(y  +  mr)  =  cos(y  -  ntr),  E[  v  +  (l/(2Ax)]=  E[(l/2)Ax  -  v],  all  the  information 
must  be  contained  in  a  frequency  band  between  0  and  (2Ax)  *  [  or  in  general 
v  and  v  +  (2Ax)  *  ]. 

B.  COMPUTATIONAL  ALGORITHMS 

The  computational  algorithms  used  must  carry  out  the  operations 
indicated  in  Eq.  (48)  as  efficiently  as  possible.  Two  algorithms  have  been 
developed,  each  with  particular  advantages. 

i.  DIRECT  SUM  ALGORITHM 

The  conceptually  simplest  way  to  carry  out  Eq.  (48)  for  a  given  value 
of  v  is  to  take  a  value  of  F(nAx),  multiply  by  the  required  cosine,  and  add 
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the  result  to  the  sum  of  the  preceding  manipulations.  This  involves  N 
multiplications  and  finding  N  cosines  for  each  output  point.  If  we  define 
a  "manipulation11  as  one  addition,  one  multiplication,  and  finding  one  function, 
to  obtain  M  output  points  from  an  interferogram  of  N  points  will  require  NM 
manipulations.  By  far  the  slowest  of  these  processes  is  the  calculation  of 
the  cosine  by  the  computer.  This  can  be  speeded  up,  however,  because  all 
the  arguments  for  the  cosines  needed  to  calculate  the  transform  at  one 
frequency  are  integrally  related  to  each  other.  This  allows  the  use  of  the 
recursion  formula 

cos(nx)  =  2  cos(x)  cos[(n  -  i)x]  -  cos[(n  -  2)x]  ;  n  =  1,2,  .  .  . 

The  recursion  formula,  however,  introduces  accumulated  roundoff  error 
into  the  cosines.  In  practice,  a  compromise  is  sought  by  periodically 
computing  a  set  of  cosines  from  the  basic  definition  while  obtaining  most 
of  the  cosines  from  the  recursion  relation.  Even  with  this  feature  the 
computer  takes  about  40  psec  per  input-output  point,  meaning  a  spectrum 
of  4000  points  in  and  4000  out  would  take  about  11  minutes. 

2.  COOLEY -TUKEY  ALGORITHM 

It  was  noticed  that  if  certain  simplifying  restrictions  were  made  on  the 
relationship  between  input  point  spacing  Ax  and  output  point  spacing  Av,  a 
number  of  similarities  existed  between  the  computation  of  Fourier  series 
and  the  factorial  experimental  designs  familiar  to  statisticians.  From  this, 
Cooley  and  Tukey  (Ref.  86)  have  developed  algorithms  suitable  for  general 
complex  Fouriev  series  computation.  Forman  (Ref.  87)  was  the  first  to 
bring  these  algorithms  to  the  attention  of  spectroscopists.  The  algorithm 
computes  the  output  spectrum  of  M  points  by  means  of  several  intermediate 
arrays,  each  of  which  also  contain  M  points.  The  trigonometric  weighting 
is  not  calculated  all  at  once  for  each  interferogram  point  but  partially  applied 
each  time  a  new  array  is  computed. 
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The  most  convenient  way  to  accomplish  the  weighting  is  to  use  complex 
exponentials,  since  the  weighting  function  of  the  sum  of  two  arguments  may 
then  be  obtained  by  the  multiplication  of  the  two  weighting  functions.  It 
becomes  very  convenient  to  also  require  that  AvAx  =  1  /M,  the  number  of 
output  point?.  If  we  let  E(v)  =  E(kAv)  =  2Ax  Real  [  T(k)],S(j)  =  F(jAx), 
j  =  1,  2,  .  .N,  and  S(0)  =  F(0)/2,  Eq.  (48)  may  be  rewritten  compactly  as 


N 


T(k)  =  £S(j)W 
j=0 


jk 


(49) 


where  W  =  exp(i2ir/M)  and  N  =  the  number  of  input  points. 

All  the  information  in  the  interferogram  can  be  obtained  if  M  =  N; 
however  it  is  often  convenient  to  have  the  output  points  more  closely  spaced. 
Forman  has  pointed  oat  that  this  can  be  effectively  done  by  extending  the 
interferogram  function  with  zeroes,  until  it  has  M  points.  In  the  following 
proof  of  the  algorithm  it  is  shown  that  the  same  effect  can  be  obtained  with 
the  calculation  of  fewer  intermediate  arrays  if  calculation  of  the  first  array 
is  modified  slightly.  The  success  of  the  method  depends  on  N  and  M  being 
highly  composite  numbers,  the  most  efficient  base  being  3.  It  may  be  shown 
however  (Ref.  86)  that  either  base  2  or  4  is  only  about  6%  less  efficient,  and 
for  some  methods  of  programming  digital  computers  a  binary-based  system 
is  much  more  convenient;  therefore,  this  will  be  the  system  used  for  our 
demonstration  of  the  algorithm. 

The  interferogram  is  required  to  have  N  =  2n  points,  extending  the 
actual  data  with  zeroes  if  necessary.  The  spectrum  is  desired  at  M  =  2m 
output  points  with  m  =  n  +  p;  m,  n,  and  p  all  being  positive  integers.  It  is 
convenient  to  express  j  and  k  explicitly  in  binary  notation 


m-i 

k  =  k  2m~1  +  k  ,2m"2...  k  =V  k.21 
m-i  m-2  o  La  i 

1=0 
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j  *  in.f'  +  kn.22°"E 


n-1 


When  the  sum  (Eq.  49)  is  rewritten  in  the  following  way,  thi  final  result  can 
be  computed  in  n  passes  through  the  data,  each  pass  involving  the  calculation 
of  M  points,  corresponding  to  the  n  summations  over  the  binary  digits. 


tw.e  z  i 

jo=0  j,=0  jn2  =  0 
n-2 


l 

.  S  soa.,  •  •  •  vw 

LW 


j  4k?. 
Jn-1 


n-i 


X  W 


j  ,k2‘ 
Jn-2 


.  .  .  W 


jok 


The  only  dependence  on  binary  d*git  j  ^  is  enclosed  within  the  square  brsokotsi 
thus  we  can  eliminate  it  from  the  sum  by  carrying  out  the  summation  anc. 
define  the  result  as  Sj. 


w“'V-  •  •  ♦vp> 


>j(k  ,  •  >  •  i  kQ,  •  *  *  » ^  ^  ®0n.i»  •  •  •  » 

P  "  jn-l 


X  W 


.iM(Vi+ 


Since  W  raised  to  any  integer  power  is  unity,  the  right-hand  factor  of  this 

equation  is  always  1,  making  Sj  depend  only  on  kQ  .  .  .  k^,  jQ  .  .  .  ,  jn„2* 

as  we  have  indicated  in  the  argument  of  Sj.  This  process  is  repeated  on  the 

next  innermost  sum  to  obtain  S,  and  so  on  until  the  last  summation  over  j 

c  o 

yields  the  desired  transform.  The  relation  between  arrays  on  two  successive 
passes  may  be  summarised  by  the  recursion  formula 
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sf+l(kp+l’ 


kl*  ko’  jn-f-2'  *  • 


=  £  Sl(kp+I-1'  •  *  *  ’  ko*  4-1-1’  •  •  *  ’  V 

Jn-f-i 

X  jn  ,  t(k  2n_1_1  +  .  .  .  +  k  ..t*1) 
w  n-/  -1\  o  p+|  / 


If  we  define  So(0, ...  0,  jn_j*  jn_2»  •  •  •  jQ)  =  S^n-1’  *  •  •  »  jQ)  =  SW*  the 
computation  can  be  carried  out  by  applying  the  recursion  scheme  Eq.  (50) 
to  the  M  points  in  the  S0  array  to  obtain  the  array,  and  repeating  the 
operation  a  total  of  n  times.  The  calculation  of  one  point  in  one  array 
requires  two  multiplications,  two  additions,  and  the  finding  of  two  exponentials. 
In  the  entire  computation  then  there  are  nM  points  to  be  computed,  meaning 

the  total  number  of  manipulations  (as 


NUMBER  Of  OUTPUT  POINTS 


defined  earlier)  is  4  nM  =  4Mlog2N. 

The  extra  factor  of  2  is  introduced 
because  of  the  complex  operations, 
making  this  relation  directly  comparable 
with  the  NM  manipulations  required  in 
the  Direct  Sum  scheme.  The  saving  in 
time  for  long  runs  may  be  seen  in 
Fig.  67,  which  is  a  comparison  of 
IBM  7090  running  time  for  the  same 
data  processed  by  two  programs,  dif- 
fering  only  in  the  algorithm  used  for  the 
transform.  A  similar  curve  has  not 
been  obtained  for  the  CDC  6000  series 


Fig.  67.  Comparison  of  computing 
times  for  the  Direct  Sum 
and  Cooley-Tukey  algo¬ 
rithms  on  the  same  data 
run  on  an  IBM  7090 
computer. 


v--’.  **• 


computers,  but  the  scaling  factor 
appears  to  be  about  6,  since  data 
requiring  4096  input  and  output  points 
takes  about  26  seconds  of  6600  central 
processor  time  as  compared  with  the 
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1,  8  minutes  required  on  the  7090  system.  The  time  required  for  a  single 
manipulation  is  still  not  the  same  for  both  algorithms,  since  by  proper 
choice  of  a  programming  scheme  only  4(M-1)  sines  and  cosines  are  needed 
by  the  Cooley-Tukey  algorithm  as  compared  with  the  NM  required  in  the 
Direct  Sum  algorithm.  For  this  number,  it  is  practical  to  calculate  each 
trig  function  from  its  argument  rather  than  use  recursion  relations  with 
their  accumulating  error.  However,  since  all  advantages  are  not  on  the 
side  of  the  Cooley-Tukey  algorithm,  it  is  reasonable  to  turn  next  to  a 
comparison  of  the  two  methods. 

3.  COMPARISON  OF  ALGORITHMS 

The  Direct  Sum  method  gives  output  points  at  intervals  Av  input  to  the 
program.  It  uses  a  minimum  of  computer  memory,  needing  only  enough  to 
store  the  N  points  of  the  interferogram  plus  the  actual  number  of  output 
points  Nq  desired  for  a  total  requirement  of  less  than  2  NQ  locations.  These 
two  features  make  this  algorithm  particularly  attractive  in  studies  of  small 
regions  of  the  spectrum  where,  for  convenience,  a  great  number  of  points 
per  resolution  width  is  desired.  No  information  is  gained  in  addition  to  that 
possible  from  interpolation  with  an  appropriate  function  in  a  spectrum  with 
the  number  of  output  points  equal  to  the  number  of  input  points,  but  it  is 
simpler  to  not  involve  additional  computer  programs.  A  disadvantage  of 
this  method  is  the  large  amount  of  time  required. 

The  Cooley-Tukey  method  solves  the  time  problem,  but  it  has  other 
disadvantages.  It  is  required  that  Av  =  1  /MAx.  Because  Ax  is  specified 
by  the  experimental  conditions  of  the  interferometer,  the  only  possible 
variation  in  A  v  is  through  the  choice  of  M;  but  since  M  must  be  an  integral 
power  of  2,  the  choices  are  limited.  The  resulting  Av  is  seldom  a  simple 
number,  which  may  be  of  some  inconvenience  if  further  hand  computation  on 
the  spectrum  is  required.  A  more  serious  limitation  is  the  memory  required 
by  this  method.  The  algorithm  calculates  the  arrays  at  all  the  points 
between  0  and  1  /L  at  a  spacing  of  Av.  As  pointed  out  earlier,  the  only 
portion  of  the  transform  of  interest  is  from  0  to  1  /2  L;  thus  M  -  2  N0  or 
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twice  as  many  points  are  calculated  as  are  desired  in  the  final  spectrum, 
requiring  2  Nq  locations.  Points  are  needed  for  calculating  the  array 
in  a  different  sequence  than  they  occur  in  the  array,  so  some  provision 
must  be  made  to  save  at  all  times  at  least  part  of  the  preceding  array  (the 
minimum  amount  convenient  for  saving  is  half).  This  requires  N0  more 
locations.  Since  each  of  these  3  Nq  locations  must  store  a  complex  number, 
with  both  real  and  imaginary  parts,  the  memory  requirements  are  for  6  Nc 
locations  or  three  times  the  memory  requirements  of  the  Direct  Sum  method. 
Because  even  large  computers  typically  have  only  about  32,  000  cells  of  core 
memory,  NQ  is  limited  at  present  to  a  maximum  of  4096  points,  in  turn 
limiting  the  smallest  spacing  between  points,  Av  ^  to  a  value  somewhat 
greater  than  would  be  convenient  on  many  occasions.  We  are  at  present 
considering  the  possibility  of  writing  a  special  program,  taking  advantage 
of  th-.?  60-bit  word  size  in  the  CDC  6000  computers,  to  store  in  the  same 
computer  word  the  real  and  imaginary  part  of  each  element  in  the  S  array, 
allowing  oar  maximum  number  of  points  to  be  doubled  and  the  minimum 
spacing  between  output  points  halved. 

Each  algorithm  then  has  its  advantages.  The  Cooley-Tukey  algorithm 
is  particularly  suited  for  scans  of  wide -frequency  regions  at  high  resolution. 
The  Direct  Sum  approach  is  particularly  suitable  for  detailed  studies  of  a 
narrow -frequency  range  with  a  small  output  point  spacing  chosen  for  con¬ 
venience  In  further  analysis. 

4.  TEST  PROBLEM 

In  checking  any  computer  program  it  is  desirable  to  have  a  problem 
simple  enough  to  be  worked  by  hand,  and  yet  one  which  tests  the  capabilities 
of  the  program.  The  function  we  have  chosen  to  check  the  Fourier  cosine 
transform  algorithms  is  a  straight  line  Tlx.)  E  (i  -  x/x  T  for  0  s  x  5  x  , 
and  F(x)  =  0  for  x  >  x  •  The  interferogram,  as  plotted  by  the  computer,  is 
shown  in  Fig.  68.  The  transform  as  obtained  by  the  computer  is  listed  in 
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IHPUT  INTERFEROGRflfl  POINT  NO. 


Fig.  68.  Test  interferogram  as  plotted  by  the  computer. 

Table  14  and  shown  as  plotted  by  the  computer  in  Fig.  69  for  Eq  =  8000, 

Ax  =  0. 1,  which  leads  to  x__  =  1.6  and  v  =  5.  The  computations  have 

me 

intentionally  been  extended  to  2v^,  so  the  periodicity  of  the  function  discussed 
earlier  could  be  demonstrated.  The  two  series  algorithms  agree  to  within 
*1  in  the  eighth  significant  figure,  and  they  agree  with  test  hand  calculations 
made  at  a  few  output  points  to  within  the  error  of  tho  cosine  tables  available 
for  the  hand  calculations. 

The  Integral  transform  for  this  simple  interferogram  can  of  course  also 
be  calculated  in  closed  form  aud  the  result  is 

E(v)  =  Eo(l  -  cos  2lr*mv>/(4,r2*mv2)  (51) 
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Fig.  69.  Transform  of  test  interferogranv  as  computed  and  plotted 
by  the  computer. 

This  has  been  evaluated  for  the  values  of  v  where  sum  results  are  available, 

21 

and  those  results  are  also  listed  in  Table  14.  From  the  above  dist  ussion 
of  computational  precision  it  is  apparent  the  differences  between  the  evaluation 
of  £q»  (501  and  the  svun  evaluation  cannot  be  due  to  computational  errors. 

The  differences  arise  rather  from  tb«  finite  sampling  of  the  interferogram, 
as  may  be  seen  from  Fig.  70,  v/hare  the.; difference*  between  the  integral 
transform  and  sum*  arc  plotted  for  4  Ax  of  d,  05  and  0.  1,  which  leads  to 

To  facilitate  comparison  bet^en  the  integral  results  and  computer  results, 
the  factor  of  2  dropped  in  obtaining  Eq  (48)  is  reinserted  as  a  multiple  of 
this  integral 


Table  14.  Comparison  of  Integral  Transform  and  Sum  Transform 
of  Straight  Line  Interferog  ram 


2E(v) 

v  .  ~“2  ?  ^  ~  cos  ^Wxm^  (computer  sum  with 

(cm-1)  W  xmV  Ax  =  0.  1  sampling) 


0 

12800.  000 

12800.000 

0.  15625 

10375.  289 

10383.  627 

0.  3125 

5187.  6446 

5204. 3434 

0. 46875 

1152.  8099 

1161.1795 

0.  625 

0 

~io"8 

0. 78125 

415.  01157 

423.44627 

0. 9375 

576. 40496 

593. 36480 

1. 09375 

211. 74060 

220.27446 

1.25 

0 

~10‘8 

1. 40625 

128.  08999 

136. 75898 

1. 5625 

207.  50578 

225.00743 

1. 71875 

85.746192 

94. 588898 

1. 875 

0 

~10~8 

2. 03125 

61. 392244 

70.450783 

2. 1875 

105.  87030 

124.23754 

2.  34375 

46. 112396 

55.433426 

2.  5 

0 

~10‘9 

2.  65625 

35.  900655 

45.  536630 

2.  8125 

64. 044995 

83.  675684 

2. 96875 

28.  740413 

38.  751 1 38 

3.  125 

0 

~10"8 

3.28125 

23.  526733 

33.  981318 

Tr  ble  14.  continued 


V 

(cm-1) 

E 

— ”3 - k  (1  -  cos  2-rrx  v) 

2  2  '  m  ' 

2ir  x  v 
m 

2E(u) 

(computer  sum  with 
Ax  =  0.1  sampling) 

3.4375 

42. 873096 

64.285108 

3.  59375 

19.  61302  3 

30.  592392 

3.  75 

0 

~10”8 

3.  90625 

16.  600463 

28.  200623 

4.  0625 

30.  696122 

54. 600961 

4.  21875 

14. 23222i 

26.  568573 

4.  375 

0 

~10"8 

4.  53125 

12.  336848 

25. 550089 

4.  6875 

23.056198 

50.485028 

4.  8437 

10. 796347 

25.060336 

5.  0 

0 

in'11 
~1 0 

POINT  NUMBER 

Fi«  70.  Difference  between  sum  and  integral  for  different  sampling 

frequencies  of  the  interferogram.  (1)  Exact  transform 

minus  sum  with  N  *  16  samples,  Ax  =  0. 1.  (2)  Exact  trans 

form  minus  sum  with  N  =  32  samples,  Ax  *  0.  05, 


interferograms  of  32  and  16  points,  respectively.  Equation  (51)  is  also 
plotted  for  comparison.  In  terms  familiar  to  spectroacopists,  these  dif¬ 
ferences  arise  because  the  filtering  in  the  "instrument"  which  determined 
the  original  interferogram,  the  straight  line,  did  not  cut  off  the  radiation 
above  =  1  /2  Ax.  This  radiation  is  being  aliared  back  at  lower  wavenumbers 
and  makes  the  energy  appear  higher  at  v  <  v  .  As  Ax  is  decreased,  in  the 
examples  shown,  the  cutoff  frequency  is  increased,  making  less  energy 
available  above  the  cutoff  frequency  for  aliasing  into  the  region  of  interest. 

C.  COMPUTER  PROGRAM 

The  computer  program  includes  facilities  for  plotting  both  the 
interferogram  and  computed  spectrum,  facilities  to  subtract  the  asymptotic 
value  from  the  raw  interferogram,  and  also  facilities  to  modify  input  points. 
The  latter  feature  makes  it  convenient  to  apply  experimental  corrections  to 
the  interferogram.  The  heart  of  the  program,  of  course,  is  the  sections 
which  perform  the  summing  of  the  series.  In  the  present  program  these  are 
subroutines  written  in  FORTRAN.  Descriptions,  flow  charts,  and  listings  of 
these  subroutines  follow. 

Identification 

FCTRAN  -  Fourier  Cosine  Transform  by  direct  sum 
6400/6600  FORTRAN 

Purpose 

To  provide  the  Fourier  Cosine  transform  of  set  of  data,  by  a  direct 
summation  of  the  series 


E(kAv)  =  Ax 


N 

F(0)4  F(jAx)  cos  iwjkAxAv 

j=i 


Method 

For  details  of  the  method,  see  Fig.  71. 
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Fig,  71.  Flow  Chart  of  Direct  Sum  Method, 
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Usage 


CALL  FCTRAN  (BEGS  DELX,  DNU,  NX,  FNUMX,  PF,  NNU,  LAV) 


where 


BEG  is  the  first  value  of  v  for  which  the  transform  is  desired. 

DELX  is  the  spacing  between  points  in  the  function  to  be  trans¬ 

formed.  The  units  of  BEG,  DNU,  and  FNUMX  are 
reciprocal  to  the  units  of  DELX. 

DNU  is  the  spacing  between  desired  output  points. 

NX  is  the  number  of  values  in  the  function  to  be  transformed. 

FNUMX  is  the  largest  value  of  v  for  which  the  transform  is  desired. 

PF  is  a  constant  phase  factor  shift  which  may  be  put  in, 

usually  zero. 

NNU  is  the  number  of  values. 

IAV  is  not  used. 


The  input  data  are  assumed  to  be  in  the  first  NX  locations  of  a  common 
array,  F.  The  output  transform  is  placed  into  the  first  NNU  locations  in  a 
common  array,  T.  The  content  of  F  at  the  end  of  the  execution  is  unchanged. 

Restrictions 

The  calling  program  must  specify  as  COMMON  three  arrays:  NOUT, 

F,  T,  with  the  dimensions  of  10,  8192,  16384,  respectively. 

Memory  Requirement 

256g  exclusive  of  COMMON. 
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SUSROUTINE  FCTRAN  (SCO.DELTtONU.NX.FMONT.FF.NNU.lAV) 
FCTRAN  DIRECT  SUM  ALOOBITMM  FOB  FOURIER  COSINE  TRANSFORM 
0ECLAR1TIVE*  FORMAT  AND  INITIALIZATION  STATENENTS 


COMMON  NOUT #F«T 

DIMENSION  FU09ft»2l»  T(il4FoZ)»  NOUTtlOl 

MINT  « 

TFlM.*§3i0931 

XNURaSEO 

INITILIZE  COSINE  TABLE  CALCULATION 

00  ?  KK»1«NNU 

MMNUMOCL* 

ABO«AHOn<FR.1.0) 

CSTI«Z . 0*COS ( AB#«TF I) 

ap9«aro*pf 

T(lt21*C0S(TBI*AR0> 

ABMAMOO ( *  ,  0*FR  1 1 .  •  >  ♦FF 
T«2t2)»Crt§<TPI«ARA) 

L»l 

CALCULATE  COSINE  TAtl.E 

00  S  K«1tNX 

IF  <M-lOfl«L>  1*3*2 

TtK»*>»CST»«T<K*l**>-T»K-f«n 

tO  TO  S 

IN0IC«100«L*1 

L*L*1 

60  TO  4 

IN0IC«1«ML 

AOt*AMOO  (FLOAT  (iNnlC)  »«i  I»t>*M 
T(INOKt?>«COS(TB!«ABO> 

CONTINUE 

*C«ICS  SUM  CALCULATION  ANO  CALCULATION  OF  NE*  FtfOUCNCT  VALUE 

T(RR*l)«n*t 
00  t  U*2«N* 

TIREtl»«F(  J)*T  <J»I*2I  ♦  T « HM •  \ » 


T(RR.li»n«L**iFU>*2.o4TiRR*tn 
I  nur*onu<*float  (RR  >  ♦«* 

T  CONTINUE 

C 

c  RETURNING  statements 

C 

ItEOMNU*! 

00  t  t»I**Ct»143t» 
t  TIIIM.4 

BfTUMN 

f  format  I4R.ITM0I8FCI  SUM  MET**) 

END 
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Identification 


FCTRAN  -  Fourier  Cosine  transform  by  Cooley-Tukey  algorithm. 
6400/6600  FORTRAN 


Purpose 

To  provide  the  Fourier  cosine  transform  of  a  set  of  data  by  the  Cooley- 
Tukey  algorithm  applied  to  the  series 


E(kAv)  =  Ax 


N 

F(0)+2^  F(jAx)  cos  2irjkAxAv 

j=i 


Method 

For  details  of  the  program,  see  the  flow  chart  (Fig.  72).  The  input 
routines  calculate  the  following  quantities:  The  number  of  input,  N  =  2n,  and 
output  points,  M  =  2m,  within  the  limitations  of  the  memory  of  the  computer. 
The  interferogram  is  filled  out  with  zeroes  if  necessary.  LD  =  2  *  M/N. 

KB  is  the  beginning  of  temporary  storage  in  the  F  array.  For  all  but  the 
largest  number  of  output  points  (M  =  8192  in  the  present  program),  this  is 
set  to  protect  the  input  interferogram.  J  =  M/4,  IBEG  =  M/2,  NBEG  =  0. 

The  calculation  of  the  array  resulting  from  one  pass  involves  boxes  i  through 
ii.  Initializing  another  pass  involves  boxes  iii  and  iv.  The  output  routines 
rearrange  the  spectrum  so  that  the  portion  of  interest  fills  in  the  lowest 
NNU  positions  in  T.  These  routines  also  multiply  by  Ax  and  insert  the  phase 
factor  if  nonzero.  The  naming  convention  of  FORTRAN  to  distinguish 
between  fixed  and  floating  point  numbers  is  observed.  Divisions  of  fixed 
point  numbers  is  assumed  truncated.  Operations  indicated  with  t  are 
complex. 

For  details  of  the  mathematics,  see  Section  B. 

the  calculated  points  when  leaving. 

NX  is  the  number  of  values  in  the  function  to  be  transformed. 
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Fig.  72.  Flow  Chart  of  Cooley-Tukey  Method. 


Usage 

CALL  FCTRAN  (BEG,  DELX,  DNU,  NX,  FNUMX,  PF,  NNU,  IAV) 

where 

BEG  is  the  first  value  of  v  for  which  the  transform  is  desired 
when  calling  this  routine.  When  returning,  it  will  be  the 
value  of  v  corresponding  to  the  value  of  E(kAv)  in  the  first 
location  of  T.  It  will  be  an  integral  multiple  of  the  value 
of  DNU  returned. 

DELX  is  the  spacing  between  points  in  the  function  to  be  trans¬ 
formed.  The  units  of  BEG,  DNU,  and  FNUMX  are 
reciprocal  to  the  units  of  DELX. 

DNU  is  the  desired  spacing  between  output  points  when  entering 
the  routine.  It  is  the  actual  spacing  between  the  calculated 
points  when  leaving. 

NX  is  the  number  of  values  in  the  function  to  be  transformed. 

FNUMX  is  the  largest  value  of  v  for  which  the  transform  is  desired. 

It  is  unchanged  by  the  subroutines  operations. 

PF  is  a  constant  phase  factor  which  may  be  inserted.  It  is 

usually  zero. 

NNU  is  the  number  of  v  values  desired  when  entering  the  program. 
It  is  the  number  of  values  beyond  BEG  for  which  the  trans¬ 
form  is  calculated  when  leaving  the  program. 

LAV  not  used. 

The  input  data  are  assumed  to  be  in  the  first  NX  locations  of  an  array 
in  COMMON,  F.  The  output  transform  is  placed  into  the  first  NNU  locations 
in  an  array  in  COMMON,  T. 
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Restrictions 


The  calling  program  must  specify  as  COMMON  three  arrays:  NOUT, 
F,  T,  with  the  dimensions  of  10,  8192,  16384,  respectively. 

The  contents  of  the  array  F  will  be  modified  by  this  subroutine  if 
effective  number  of  output  points  is  8192. 

Memory  Requirement 

552g  exclusive  of  COMMON  storage. 
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APPENDIX 


INTERFEROMETRY  AND  CONVENTIONAL  SPECTROSCOPY 


In  theory,  two-beam  interferometry  has  a  decided  S/N  advantage  over 
conventional  single«slit  diapersive  spectroscopy  when  detector  noise  is 
the  limiting  factor  (see  Section  1,  Introduction).  In  a  paper  on  FIR 
instrumentation,  Richards(Ref.  14)  describes  an  experimental  comparison 
of  a  Michelson  and  a  lamellar  grating  interferometer  with  a  conventional 
grating  spectrometer.  With  respect  to  the  comparative  performance  of 
these  two  types  of  instruments  in  the  FIR,  Richards  said, 

"We  experimentally  verify  the  Fellget  advantage  of  the  inter¬ 
ferometer  to  within  about  a  factor  of  3  in  time.  This  is  of  the 
order  of  the  uncertainties  in  our  estimates  of  noise  levels  and 
grating  and  filter  efficiencies.  " 

This  statement  was  disputed  by  Kneiibuhl,  Moser,  and  Steffen  (Ref.  55)  in 
a  paper  describing  the  construction  and  performance  of  a  conventional  FIR 
grating  spectrometer.  They  conclude, 

"From  these  spectra,  we  conclude  that  the  performance  of 
Fourier-transform  spectrometers  may  be  slightly  better  than 
that  of  grating  spectrometers,  but  the  difference  is  certainly 
smaller  than  that  stated  by  Richards.  The  spectral  resolution 
yet  obtained  by  bath  types  of  spectrometers  lies  between  0.  i  and 
0.  05  cm”*  in  the  spectral  region  between  30  and  60  cm"**.  " 

This  statement  was  challenged  by  Dowling  and  Hall  (Ref.  881,  who  said, 
"While  achieving  remarkable  performance  with  their  spectrometer 
wc  feel  that  some  comments  and  conclusions  arrived  at  by  KMS 
with  respect  to  the  comcarative  performance  of  these  two  types 
of  instruments  are  not  supported  by  the  evidence  presented  and, 
furthermore,  are  in  disagreement  with  our  experience  with  a  far- 
infrared,  lamellar  grating  interferometer.  " 


And  alto  with  respect  to  the  resolution  vet  obtained  by  both  types  of 
instruments  they  said: 

"To  the  best  of  our  knowledge  this  resolution  has  not  yet  been 
demonstrated  for  conventional  grating  spectrometers  and  has 
only  recently  been  demonstrated  for  interferometric  instruments.  " 
More  details  are  given  in  these  and  other  papers  (Refs.  15,  89,  90). 

With  the  exception  of  the  detector  employed,  the  opucal  parameters 
of  the  conventional  grating  spectrometer  (Ref.  55)  and  the  Aerospace 
lamellar  grating  interferometer  (Ref.  15)  are  quite  similar.  These  are 
given  in  Table  A-l . 

In  the  previous  comparison  (Refs.  55,  88,  90),  there  was  quite  a 
large  difference  in  experimental  conditions  (pressure  and  absorbing  path 
length)  in  the  respective  water  vapor  spectra.  Presented  here  are  data 
obtained  with  the  lamellar  grating  interferometer  under  experimental  condi¬ 
tions  closer  to  those  used  by  KMS  (Refs.  55,  90). 

22 

We  obtained  two  single-beam  runs  for  I^O  with  the  lamellar  grating 

interferometer  with  an  absorbing  path  length  of  ~9.  2  m  and  a  pressure  of 

"1-1/2  Torr.  The  first  of  these  was  a  reasonably  high  resolution  run 

(Avm  -  0.  10  cm-*)-  The  transformed  interferogram  (i.e. ,  the  spectra) 

showed  the  triplet  at -38  cm-*  quite  clearly  resolved  (Ref.  88),  indicating 

the  run  had  good  resolution  and  reasonable  sensitivity  (the  line  at  38.  642 

cm-*  is  quite  weak  and  close  to  the  strong  line  at  38.792  cm-*).  However, 

the  "new"  lines  seen  by  KMS  (Ref.  55)  at  46.3,  49.5,  and  50.25  cm'*  were 

23 

not  observed.  Although  the  apparent  S/N  for  this  run  was  somewhat  poorer 


In  order  to  obtain  a  long  absorbing  path  (to  correspond  to  the  long  absorb¬ 
ing  path  used  by  KMS)  the  entire  optical  path  (from  detector  to  source)  was 
used  by  admitting  HgO  to  the  whole  system.  This,  of  course,  precluded 
double-beam  ope  ration,  and  thus  we  encountered  the  dynamic  range  problem. 

23 

We  believe  this  is  due  in  part  to  the  dynamic  range  problem. 


Table  A"i.  Comparison  of  a  Conventional  Grating  Spectrometer' 
and  a  Lamellar  Grating  Interferometer*5 


Conventional 

Spectrometer 

Interferometer 

Grating  area 

1050  cm^ 

930  cm^ 

f- number 

-f/4 

-f/4 

Slit  or  entrance 
aperture  area 

0.  8  cm'’ 

0.  7  cm^ 

Detector  and  NEP 

Golay 

4  4  ^ 

(8  X  10  ) 

Ge  Bolometer  T  =  4.  2°K 

r* 

(5  X  10‘1?*) 

*Ref.  55 
bRef.  15 

cThis  is  from  advertising  literature  from  Eppley  Laboratories  (Golay)  and 
Texas  Instruments  (Ge  Bolometer).  It  has  been  our  experience  in  comparing 
our  bolometer  and  Golay  cells  that  the  bolometer  is  ~ 6  to  9  times  better 
NEP-wise  than  the  factor  of  16  deduced  from  the  advertising  literature. 


than  in  double -beam  runs  (Ref.  88),  we  feel  it  was  sufficiently  high  to  detect 
these  lines  if  the  relative  intensities  as  shown  in  Fig.  3  of  Ref.  55  are  any¬ 
where  near  correct.  Their  absence  confirms  our  contention  (Ref.  88)  that 
these  lines  are  second-order  water  vapor  lines  (and  they  do  indeed  correspond 
exactly  within  the  accuracy  of  the  wavelength  quoted  by  KMS)  appearing  in 
the  first-order  spectrum  because  of  incomplete  filtering. 

A  second  run  was  made  in  a  much  shorter  time.  The  optical  path 
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difference  attained  was  -5.0  cm  (Av.,  ~  0.  20  cm  )  in  increments  of  30p. 
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of  optical  path  difference  and  each  point  was  sampled  for  0.75  sec.  This 
interferogram  function  was  FGT  twice,  once  with  L  =  5.  0  cm  (total  time  of 
run,  including  dead  time  =41.8  min)  and  also  with  L  =  3.  5  cm  (by  deleting 
the  last  1,  5  cm  of  the  interferogram,  which  corresponds  to  a  running  time 
of  29.3  min).  The  55-56  cm  *  regions  of  these  spectra  are  compared  with 
the  corresponding  region  as  given  by  KMS  (Fig.  4  of  Ref.  55)  in  Fig.  73. 

The  experimental  conditions  are  listed  in  Table  A-2.  The  results  of  KMS 
are  "normalized  transmittance.  "  We  have  "normalized"  our  curves  by 
assuming  that  the  greatest  intensity  in  this  region  corresponds  to  100% 
transmission  (and  this  appears  to  be  what  KMS  have  done).  We  make  the 
following  observations: 


1.  In  Fig.A-1,  the  depth-of-absorption  curves  B  and  C  are  much 
greater  than  in  A. 

2.  The  B  and  C  lines  appear  much  more  symmetrical  than  in  A, 
especially  the  high  wavenumber  component. 


3.  In  A  the  low  wavenumber  component  appears  slightly  stronger  than 
the  high  wavenumber  component.  In  B  -ad  C  the  opposite  is  true. 
Theoretical  relative  intensities  predict  the  higher  wavenumber 
component  to  be  stronger  (Ref.  56). 


We  point  out  that  Ax  =  30  p  corresponds  to  v  = 
chosen  since  we  had  changed  our  filtering  slightly, 


166  cm"*.  This  Ax  was 
and  to  be  sure  of  avoiding 


aliasing  we  chose  Ax  =  30  p.  Actually  we  could  have  used  Ax  =  40  p  as  it 
wag  apparent  from  the  spectrum  Avc  ^  125  cm’*,  This  would  have  resulted 
in  a  -33%  reduction  in  time  with  small  loss  in  the  S/N  in  the  spectrum. 
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The  dead  time  (i.e.  the  time  it  takes  to  move  from  point  to  point)  during 
which  no  data  are  taken,  tor  30  p  optical  path  difference  steps,  is  also  0.  75 
sec.  Thus  the  efficiency  of  this  run,  measuring  time  per  total  time,  was  -  50%. 
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Table  A-2.  Comparison  of  Experimental  Parameters 
Illustrated  in  Figure  71 


A  (KMS) 

B 

C 

Absorbing  oath 

7.4  m 

9.2  m 

9.2  m 

Pressure 

2  Torr 

1-1/2  Torr 

1-1/2  Terr 

Computer  time 

0 

~  44  sec3 

-44  sec3 

Scan  time  for 

55-56  cm  *  region 

17.3  min 

0.  42  min^ 

0.  29  min^ 

Computer  is  a  CDC  6600.  The  actual  FCT  time  is  small  compared  with 
the  read-in,  read-out  time,  etc. 

^Bascd  on  the  total  time  it  took  to  get  the  inte  rferogram  and  the  fact  that 
100  cm  *  was  covered.  See  text  for  a  discussion  on  KMS's  objection 
to  this. 
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4.  Th  °  apparent  noise  level  increases  from  B  to  C,  whereas  one 

would  expect  it  to  decrease  (see  Section  V).  We  note  in  particu¬ 
lar  that  maximum  intensity  filiation  in  the  55.  8— *-56.  0  cm"! 
more  than  doubles  in  going  from  B  to  C  (no  H^O  absorption  line 
is  expected  heie).  We  interpret  these  variations  (and  they  are 
seen  in  proximity  to  the  other  strong  and  medium  lines)  as  being 
the  aidelobes  discussed  in  Section  II.  We  thus  infer  that  these 
lin  ;s  are  narrower  than  the  scanning  functions  of  these  two  runs 
and  the  variations  observed  are  not  true  noise.  We  make  no 
comment  on  the  noise  level  for  the  spectra  shown  in  Fig.  71 -A 
since  KMS  made  no  estimate  of  their  S/N 

Kneubiihl,  Moser,  and  Steffen  (Ref.  89)  have  objected  to  the  "equivalent 
scan  time"  used  by  Dowling  and  Hall  (Ref.  88)  in  comparing  results  obtained 
with  the  two  instruments,  on  the  basis, 

"The  equivalent  scan  time  introduced  by  Dowling  and  Hall  does 
not  seem  to  be  a  great  value.  Theoretically,  high  resolution 
of  a  srriaii  spectral  region  can  be  achieved  by  an  aperiodic 
interferometer  in  a  short  time  (e.g,,  i  min).  But  this  requires 
a  narrow-  band  filter  with  a  bandwidth  of  the  order  of  magnitude 
of  1  cm"  ,  e.g.,  a  diffraction-grating  spectrometer." 

Although  a  diffraction-grating  spectrometer  would  serve  admirably  as  a 
narrow  bandpass  filter  (as  long  as  higher  order  radiation  is  filtered  out),  such 
a  complicated  instrument  may  not  be  necessary.  The  progress  in  fabricating 
narrow  bandpass  filters  (Ref.  9)  for  the  FIF  is  quite  encouraging  in  this 
respect. 

In  his  well-known  treatment,  Fellgett  (Refs,  3,  4)  has  shown  that  the 
S/N  gain  of  a  multiplex  instrument  over  conventional  scanning  spectrometers 
is  v*N/2.  where  N  is  the  number  of  spectral  elements  (resolution  widths)  to 
be  studied.  The  advantage  is  1  wher,  N  =  4.  That  is.  all  other  quantities 
being  equal,  the  interferometer  Has  the  advantage  if  N  >  4,  the  conventional 
spectrometer  when  N  *  4.  Usually  many  more  than  four  spectral  elements 
are  desired. 

At  the  present  time  the  conventional  spcct remete  r  offers  the  advantage 
of  short  wavenumber  scans ,  but  with  the  development  of  narrow  bandpass 
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filters  (Ref.  91)  and  the  more  efficient  use  of  the  time  available  for 

measurement  (i.e.,  minimization  of  the  dead  time),  we  believe  this  advantage 

26 

is  a  temporary  one. 
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An  opinion  expressed  in  a  letter  to  J.  F.  Moser. 
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